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OPERA TION FARSIDE | 


From ‘off the shelf 
_ thousands of miles 


Balloon-carried to 100,000 feet and launched up 
through the big bubble, the four-stage Operation 
Farside missile achieved the greatest height @ver at- 
tained by a man-made object—an estimated 4,000 
miles to the edge of outer space. 

A cluster of four Thiokol solid propellant Recruit 
engines supplied the critical first stage take-off thrust, 
followed by a single Recruit as the second Stage. 

Farside again emphasized the flight reliability of 
Thiokol solid propellant rocket engines, whether 
ground, sea or air-launched. And the five Recruit 
engines were standard “off-the-shelf” models. The 
Recruit has been proven time after time for power 
and performance. Since the engine is a regular pro- 
duction item, precious time and dollars were saved 

_ for the Farside project. 


Thiokol: 


CHEMICAL CORPORATION 


TRENTON. N.J.- ELKTON. MD 
HUNTSVILLE, ALA.- MARSHALL. TEXAS 
MOSS POINT. MISS.+ BRIGHAM CITY. UTAH 


Chemical Corporation for its liquid polymers, rocket propellants, plasticizers and other chemical products. _ 


> 
"Registered trademark of the Thiokol | 
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GUIDANCE and CONTROL 


in SPACE TECHNOLOGY 


It is becoming increasingly apparent that many of the 
techniques and analyses, and much of the equipment, 
developed for the present Air Force ICBM-IRBM pro- 
grams will have a wide future application in space 
technology. For instance, many of the guidance and con- 
trol techniques for ICBM’s are applicable to the space 
vehicles of the near future. 

An important element of these applications is preci- 
sion. The pre¢ision required of the guidance and control 
system for vehicles aimed at the moon or one of the 
planets is not substantially greater than that required for 
the Air Force ICBM-IRBM programs. And, the precision 
needed to guide a vehicle into a near-circular orbit of 
Earth is even less than that required for ICBM’s. 

The problem of communication with lunar and plane- 
tary vehicles is, of course, made more difficult by the 
much greatet distances involved. This, however, is not 
an insurmountable difficulty if today’s trends continue 
in the use of higher transmitted power, narrower com- 
munication bandwidths and amplifiers with very low 
noise-figures. 

The problems of operating electronic equipment in the 
space beyond our atmosphere are already encountered 
on present ballistic missile trajectories. The principal dif- 
ference in the case of space vehicle applications is the 


requirement for longer equipment lifetimes. Electronic 
equipment and power supplies will have to last for several 
hours or days or weeks, instead of a few minutes, under 
conditions of vacuum pressure, zero “g” fields, and bom- 
bardment by micrometeorites, “ect i particles, and 
radiation. 

The preceding examples serve to illustrate some of the 
ways in which the ICBM-IRBM programs are advancing 
the basic techniques of space technology. 

Since 1954, Space Technology Laboratories has been 
providing over-all systems engineering for these pro- 
grams. Both in support of this responsibility and in antici- 
pation of future system requirements, the Laboratories 
are presently engaged in a wide variety of advanced 
analytical and experimental work directed toward the 
exploration of new approaches in space vehicle elec- 
tronics, propulsion, and structures. 


The scope of STL’s work requires a staff of unusual 
technical breadth and competence. Engineers and scien- 
tists who are interested in advanced experimental devel- 
opment projects (as distinct from development for 
manufacturing, in which STL is not engaged) are invited 
to investigate the many opportunities on the Laboratories’ 
Technical Staff. 


ae Space TECHNOLOGY LABORATORIES 


A Division of The Ramo-Wooldridge Corporation 


5730 ARBOR VITAE STREET 


LOS ANGELES 45, CALIFORNIA 


a 
FEBRUARY 1958 73 


| 4 publication of the 


—-AMERICAN ROCKET 


Research and Development 


IRWIN HERSEY—DIRECTOR OF PUBLICATIONS 


EDITOR MARTIN. SUMMERFIELD 


ASSISTANT EDITOR LARKIN: JOYNER 


ASSOCIATE EDITORS 
ALI BULENT CAMBEL, Northwestern University 
IRVIN GLASSMAN, Princeton University 
M. H. SMITH, Princeton University 
CONTRIBUTORS 

MARSHAL FISHER, Princeton University — 
GEORGE F. McLAUGHLIN 
ADVERTISING PRODUCTION MANAGER 
WALTER BRUNKE 

ADVERTISING & PROMOTION MANAGER | 
WILLIAM CHENOWETH 


ADVERTISING REPRESENTATIVES 


James & Co. 
6535 Wilshire Blvd., Los Angeles, Calif. 
Telephone: Olive 3-3223 


6B C. Emery & Associates 
—155 East 42 St., New York, N. Y. 
Telephone: Yukon 6-6855 


& Summers & Associates R. F. and Larry Pickrell 
35 E. Wacker Dr., Chicago, Ill. 4 3 18 Stephenson Bldg., Detroit, Mich. 
Telephone: Andover 3-1154 Telephone: Trinity 1-0790 


Louis J. Bresnick 
304 Washington Ave., Chelsea 50, Mass., Telephone: Chelsea 3-3335 


AMERICAN ROCKET SOCIETY : 
Founded 1930 


OFFICERS 


George P. Sutton 
John P. Stapp 
James J. Harford 

A. C. Slade 

Robert M. Lawrence 
Andrew G. Haley 


President 

Vice-President 

Executive Secretary 

Secretary 

Treasurer 

General Counsel 


ier BOARD OF DIRECTORS 
Terms expiring on dates indicated 


Krafft Ehricke, 1959 H. S. Seifert, 1958 
S. K. Hoffman, 1958 K. R. Stehling, 1958 
Simon Ramo, 1960 Martin Summerfield, 1959 
H. W. Ritchey, 1959 Wernher von Braun, 1960 
Maurice J. Zucrow, 1960 he ty 


ae TECHNICAL DIVISION CHAIRMEN 


John F. Tormey, Propellants and 
Combustion 


‘David G. Simons, Human Factors 
Lawrence S. Brown, Instrumenta- 
tion and Guidance Brooks T. Morris, Ramjet 
Edward N. Hall, Liquid Rocket William L. Rogers, Solid Rocket 
Krafft A. Ehricke, Space Flight 


Scope of JET PROPULSION 


This Journal is a publication of the Ameri- 
can Rocket Society devoted to the advance. 
ment of the field of jet propulsion through 
the dissemination of original papers disclosing 
new knowledge or new developments. As 
used herein, the term ‘‘jet propulsion” 
embraces all engines that develop thrust by 
rearward discharge of a jet through a nozzle 
or duct; and thus it includes air-consuming 
engines "and underwater systems as well as 
rockets. JET PROPULSION is open to con- 
tributions dealing not only with propulsion 
but with other —* of jet-propelled flight, 
such as flight mechanics, guidance, teleme'ter- 
ing, and research instrumentation. I[n- 
creasing emphasis will be given to the scientific 
problems of extraterrestrial flight id 


Information for Autl ors 


Manuscripts must be as brief as the proper 
presentation of the ideas will allow. hx- 
clusion of dispensable material and concise- 
ness of expression will influence the Editors’ 
acceptance of a manuscript. In terms of 
standard-size double-spaced typed pages, a 
typical maximum length is 22 pages of text 
(including equations), 1 page of references, 
1 page of abstract, and 12 illustrations, 
Fewer illustrations permit more text, and vice 
versa. Greater length will be acceptable 
only in exceptional cases. 

Short manuscripts, not more than one 
quarter of the maximum length stated for full 
articles, may qualify for publication as 
Technical Notes. They may be devoted 
either to new developments requiring prompt 
disclosure or to comments on previously pub- 
lished papers. Such manuscripts are usually 
published within two months of the date of 
receipt. 

Sponsored manuscripts are _ published 
occasionally as an ARS service to the indus- 
try. A manuscript that does not qualify for 
publication according to the above-stated 
requirements as to subject scope or length 
but which nevertheless deserves widespread 
distribution among jet propulsion engineers, 
may be printed as an extra part of the Journal 
or as a special supplement, if the author or 
his sponsor will reimburse the Society for 
actual publication costs. Estimates are 
available on request. Acknowledgment of 
such financial sponsorship appears as a 
footnote on the first page of the article 
Publication is prompt since such papers are 
not in the ordinary backlog. 

Manuscripts must be double spaced on one 
side of paper only with wide margins to allow 
for instructions to printer. Include a 100 to 

word abstract. State the authors’ 
positions and affiliations in a footnote on the 
first page. Do not type equations; write 
them in ink. Identify unusual symbols or 
Greek letters for the printer. References are 
to be grouped at the end of the manuscript 
and are to be given as follows: for journal 
articles: authors first, then title, journal, 
volume, year, page numbers; for books: 
authors first, then title, publisher, city, edi- 
tion, and page or chapter numbers. Line 
drawings must be clear and sharp to make 
clear engravings. Use black ink on white 
paper or tracing cloth. Lettering should be 
large enough to be legible after reduction. 
Photographs should be glossy prints, not 
matte orsemi-matte. Each illustration must 
have a legend; legends should be listed in 
order on a separate sheet. 

Manuscripts must be accomplished by 
written assurance as to security clearance in 
the event the subject matter lies in a classified 
area or if the paper originates under govern- 
ment sponsorship. Full responsibility rests 
with the author. 

Submit manuscripts in duplicate (original 
plus first carbon, with two sets of illustra- 
tions) to the Editor, Martin Summerfield, 
Professor of Aeronautical Engi i | Prince- 
ton University, Princeton J. Preprints 
of papers presented at ARS national meetings 
are automatically considered for publication. 


JET PROPULSION is published monthly by 
the American Rocket Society, Inc., and the 
American Interplanetary Society at 
Northampton Sts., Easton, Pa., U. 
Editorial offices: 500 Fifth "Ave., New Mork 
36, N. Y. Price: $12.50 per year, $2.00 
per single copy. Second-class mail privi- 
leges authorized at Easton, Pa. Notice of 
change of address should be sent to the 
Secretary, ARS, at least 30 days prior to 
publication. Opinions expressed herein are 
the authors’ and do not necessarily reflect the 
views of the Editors or of the Society. 
Copyright 1958 by the American Rocket 
iety. 
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Frederick F. Liu is presently executive 
vice-president of Dresser Dynamics, Inc., 
a subsidiary of Dresser Industries. Dr. 
Liu received his B.S.E.E. summa cum 
laude, China; Dipl. Ing., Germany; 
B.S.M.E. from CIT; and his doctorate 
from Princeton University. Entering the 
rocket research field at the Forrestal 
Research Center, he designed electronic 
and rocket measurement systems and co- 
developed the Li-Liu pickup there. While 
a research engineering specialist of 
Rocketdyne, he directed the technical 
activities of propulsion physics laboratory 
and made notable contributions in high 
speed rocket control and measuring de- 
vices, and shock wave and hypersonic 
studies. He serves as a consultant to 
AEC’s KEWB program. 
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UCLA, was research engineer at Collins 
Radio, and later a senior research engi- 
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neer at Rocketdyne. He has several 
patents in the electronic field and is cur- 
rently the head of the electronic division 
of Dresser Dynamics, Inc. 


Introduction 


HE development of methods for 

measuring dynamic pressures has 
been spurred in recent years largely as a 
result of the increasing requirements of 
combustion and other propulsion kinet- 
ics research efforts (1).! Prior to this, the 
main efforts in the development and use 
of such equipment were in the internal 
combustion engine research fields where 
significant progress was made. At the 
present time, however, even these earlier 
requirements are being exceeded by the 
demands of high velocity shock wave, 
detonation and other aerothermody- 
namic studies (2). _These demands have 
the same basic objective; namely, the de- 
termination of the true magnitude of 
pressure relating to the physical phenom- 
ena under study and the recording of 
its pattern of change with time. They 
differ, however, in speed requirements. 
Phenomenological changes can take 
place within a period of microseconds or 
shorter: Their pressure-time wave forms 
often have to be determined with a maxi- 
mum of accuracy and precision; in other 
words, with a minimum distortion in am- 
plitude, frequency and phase response. 

In this paper a review of some recent 
advances in pressure measurement as it 
pertains to rocketry, gas dynamics and 
nuclear kinetics research are presented 
together with a critical survey of a num- 
ber of recent developments. A_ brief 
description of several techniques de- 
veloped by the authors for measuring 
high speed pressure transients is also 
included. 

The problems involved in dynamic 
pressure measurement are invariably 
more complex than those associated 
with static pressure measurement. In 
speaking of these problems, it is no 
longer possible to confine our attention 
solely to gages and pickups. 

Presented at the ARS Semi-Annual Meet- 
ing, San Francisco, Calif., June 10-13, 1957. 

1 Numbers in 


parentheses indicate 


References at end of paper. 


Recent Advances 
In Dynamic Pressure Measurement Techniques 


It is well understood that whenever 
there is pressure propagation through a 
confined channel, there will be associated | 
pressure attenuation due to frictional 
losses and accompanying modification — 
of the time constant of the original pres- _ 
sure waveform. Moreover, since pres- 
sure is related to the rate of deforma- 


duces another time constant and corre- _ 
sponding attenuation factor. If 
length of tubing is introduced in the © 
measuring system, the acceleration of 
the gas as well as the reflection at the © 
end of the tubing will further distort | 
the waveform. These phenomena are 
well known to fluid dynamicists and 
recently have been studied rigorously by 
a number of workers including Iberall 
(3). From the applied standpoint, 
rocket engineers have since found out 
that the thrust chamber pressure pat- 
tern as measured through a length of 
tubing and as measured by a pickup 
mounted flush in a chamber wall may — 
appear altogether different. The flush 
mounted pickup is, of course, capable — 
of giving the truer pressure pattern. 
The majority of combustion phenom- 
ena pertaining to inflammable mis- 
tures may be considered either as de- 
flagration (subsonic) waves, or as det- 
onation (supersonic) waves involving 
shock. Shock fronts have pressure rise 
times estimated in the order of milli- 
microseconds; high frequency oscilla- 
tions in a small rocket chamber can be 
6000 eps or more. When these pressure 
changes are sensed by systems without 
the requisite transient and frequency 
responses, the measurements will be 
true neither in magnitude nor in “pres- 
sure vs. time’’ waveform. In the space 
domain a shock front which propagates 
at a velocity many times the speed of 
sound has a thickness in terms of mean 
free paths (20.5 X 10~* em for N: at 11 
atm at Mach 1.128) (4). A 6000 cps 
pressure osvillation propagating at a 
velocity of 3600 fps and at a tempera- 
ture of between 4000 and 5000 F has a 
wave length of 7.2in. Due to the tran- 
sit time effect, a pressure pickup which 
has a circular diaphragm diameter of, 
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for example,? 7 in. and which is 

mounted perpendicular to the propaga- 

tion direction, will not yield the correct 

_ transient waveform since at any time 

_ the indicated pressure to a step function 
pis given by 


A 
Pind = KF = xf pdA 
0 
“ay (x/c) =t 
A= 2f (r? — c%?)'/2 edt = 


(x/c) =t 
f [(r?/e?) — dt 


—r/c 


where pina and p are the indicated and 
actual pressure; A is the diaphragm 
area with radius r; c is the propagation 
_ velocity of the phenomenon; F is the 
pressure-summed force exerted upon the 
_ sensing element; ¢ denotes the time; 
_ and K isa constant of proportionality. 
Thus, in the case of measuring 6000 
_ eps pressure oscillations, the diaphragm 
covers a length of only about 7¢ of the 
~ wave length or angularly 22 deg. Its 
indicated value at any time, then, is the 
- mean squared value over this 22 deg are. 
Rocket Engine Dynamic 
Pressure Pickups 


Certain rocket propellant combina- 
tions have already attained combustion 
temperatures above 5000 F, 
higher combustion temperatures 
insight. Thus, the development of 
a dynamic pickup which can operate 
—ailemadans at such high temperatures 
for a reasonably prolonged period and 
withstand acceleration as high as 1000 g 
without sacrificing appreciably its meas- 
urement stability is of practical im- 
portance to rocket research. One of the 
major contributions in this direction, 
attributed to the cooperative effort be- 
tween Princeton University and MIT, is 
the development of the flush-mounted 
water-cooled Li-Liu pickup. By means 
of liquid cooling, errors due to thermo- 
elastic sources are kept to a minimum 
yet the additional coolant mass does not 
contribute excessive acceleration effects. 
During the last few years this pickup has 
proved itself to be a useful tool among 
rocket research institutions (5). In its 
earliest stage of development, barring 
any negligence in the cooling system, the 
Li-Liu pickup had successfully with- 
stood many rocket firings with esti- 
mated temperatures ranging from 4280 
to 6200 F under both stable and scream- 
ing conditions. The heat transfer rate 
across the diaphragm as determined by 
Tischler and his colleagues at NACA 
Lewis Laboratory was 3.92 Btu/ 
in.*-sec for a screaming run. During 
* A few pressure pickups, for instance 
the sensors produced by Omega Instru- 
ment Company, Pasadena, have dia- 
phragm diameters down to 3? in. 
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this firing the estimated temperature 
was 6200 F; the water temperature at 
the entrance of the pickup’s cooling 
system was 64 F; and at the outlet 
87.3 F; the coolant flow rate was 0.028 
Ib/sec.* With later adjustments, and 
with calculations based on the area of 
the diaphragm alone, a heat transfer 
rate of 7.6 Btu/in.?-sec was attained. 
Newer models of super Li-Liu pickups 
‘an now achieve even higher heat trans- 
fer rates. 

The use of water-cooling over thin 
diaphragms has certainly resolved many 
previous difficulties. Without efficient 
cooling a diaphragm made of materials 
now available cannot endure the flame 
temperature even for a brief period; the 
diaphragm would be warped (outward) 
and the sensing element be damaged so 
severely as to render dependable 
measurement virtually impossible. High 
temperature also causes drastic changes 
in the elastic moduli of the sensing ele- 
ment. The resultant signal output then 
becomes erroneous by the standard 
of room-temperature calibration. These 
advantages being apparent, it is never- 
theless necessary to make an analysis as 
to whether the coolant also may exert 
any adverse influence. 

Understandably, the coolant be- 
tween the diaphragms would constitute 
an additional mass to the moving parts. 
Consequently, the effects are an increase 
in the inertial term resulting in greater 
overshoot and, concurrently, a decrease 
in natural frequency. The liquid cool- 
ant also exerts some damping effect, 
causing the damping ratio to increase. 
These effects are evidenced in Figs. 1 and 
2 which show the response of the Li-Liu 
pickup and Photocon pickup to identi- 
cal step functions from a pneumatic 
shock. These oscillograms were taken 
both with and without the cooling 
water present. (Note the sweep rate of 
the trace is 20 microsec per centimeter 
for the Photocon pickup, and 10 micro- 
sec per centimeter for the Li-Liu pick- 
up.) The cooling water causes some 
additional overshoot in both the pick- 
ups; the rise time has been increased 
from approximately 6 microsec to 10 
microsec in the case of the Li-Liu pick- 
up, causing its natural frequency to 
drop from 31.2 ke to approximately 18.7 
ke. For the Photocon pickup the rise 
time increased from approximately 16 
to 20 microsec resulting in a decrease of 
natural frequency from 25 ke to ap- 
proximately 20 ke. These tests, how- 
ever, were conducted at extreme con- 
ditions with very fast driving functions. 
During lower frequency operation the 
cooling water apparently has little or 
no detrimental effects. 

In addition to the Li-Liu pickup and 


3 Communication from A. 0. Tischler 
to Forrestal Research Center, Princeton 
University. 


Photocon pressure indicator, several new 
water-cooled pickups show promise for 
rocket engine uses. The E. Broza 
water-cooled Drucksgeber, made in Tett- 
nang, West Germany, is a quartz-type 
pressure transducer; the Omega-Liu 
pickup utilizes a different cooling tech- 
nique, lateral cooling of small-diameter 
diaphragms. 


Developments in Pressure 
Transducer Design 


Sensing Elements and Methods of 
Transduction 


Numerous new types of transducing 
methods for pressure pickups have ap- 
peared in recent years (6). In rocket 
testing the resistive strain gage and ca- 
pacitance types of pickups are still used 
most extensively for high frequency 
combustion chamber pressure measure- 
ments. The Li-Liu pickups use the 
former principle and the Photocon is an 
example of the capacitance type. The 
bonded strain gage, particularly when 
employing mutually crossed type grid 
windings on a strain tube, has excellent 
point-by-point temperature compensa- 
tion and is characterized by low thermal 
drift over a wide temperature range. 
First order thermal effects can be re- 
duced to virtually an insignificant mag- 
nitude by such compensation tech- 
niques. The wire gages add very little 
mass to the moving elements. They 
contribute a negligible amount of ac- 
celeration sensitivity and do not mate- 
rially worsen the dynamic performance. 

Good linearity is another advantage 
of the well designed strain gage sensing 
element if, during the dynamic measure- 
ment, the strain Al/! of the elastic ele- 
ment is kept below 0.1 per cent. When 
coupled to its associated electronic sys- 
tem, the low output impedance of resis- 
tive strain gages is a desirable feature 
because a typical strain gage with resist- 
ance on the order of 200 ohm has a 
capacitance of about 10 mmf and an in- 
ductance of 0.015 microhenry. The 
time constant of the gage itself is thus 
only about 0.002 microsec. To this, 
however, must be added the capaci- 
tances of the excitation sources (on the 
order of 100 mmf) and the capacitance 
of the lead cable (which depends upon 
the length and type of cable used). 
These reactances and the stray capaci- 
tances should contribute relatively mi- 
nor factors to the frequency and transient 
responses of a shielded, d-c excited 
bridge using resistors of low residual 
reactances. The upper frequency limit 
(wim) of a strain gage bridge, namely, 
the 3db point, can be estimated from 
the relation: wim = 1/RC. Assuming a 
200-ohm strain gage pickup having a 
total active and stray capacitance value 
of 500 mmf, then, electrically, the strain 
gage itself should be able to reach an 
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upper frequency limit of about 1 mega- 
cycle. This, nevertheless, is not prac- 
tically obtainable due to the limiting 
dynamic response of the mechanical 
part of the sensing element. The latter 
limits most of the commercially avail- 
able strain-gage pickups to below 100 ke 
natural frequency. One problem of the 
bonded-type strain gage is its instability 
under irradiation by strong neutron 
flux for prolonged periods. There have 
been reports of erratic changes of strain 
gage output up to 10 per cent. The 
precise nature of such radiation effect, 
however, has yet to be studied (7). 
Many of the advantages of the bonded 
strain gages are shared by the unbonded 
strain gages barring such exceptions as 
the greater residual reactance values 
and the greater sensitivity to accelera- 
tion of the latter (8). The use of un- 
bonded strain gage sensing elements is 
predominant, however, for low pressure 
measurements. disadvantage of 
both types of strain gages is the rela- 
tively low signal output, generally in 
the millivolt range even when used in a 
bridge circuit with 6 to 12 volt excita- 
tion. Several promising approaches to 
improve the strain gage output have 
been made available recently. Statham 
Laboratories has developed the Zero- 
length unbonded strain gage which 
yields an output of up to 3 volt with 
12-volt excitation (9). Recently in 
co-operation with the senior author and 
under the auspices of a program con- 
ducted by the AEC and Atomics In- 
ternational, Statham Laboratories and 
Dynamic Instrument Company have 
successfully developed newer types of 
unbonded gage transducers capable of 
higher transient response and reasonably 
suitable for reactor use. 

Electronic approaches to this problem 
include the use of pulse excitation meth- 
ods, such as those reported by Yates et 
al. in England and by Sander and 
Brodie of NACA (10). With such 
techniques, a short pulse of excitation 
voltage up to 200 volts can be applied 
to the bridge momentarily to obtain a 
high output. In addition, several tech- 
niques producing frequency modulation 
by means of resistance change have been 
developed, including one by Berwin. 

The condenser type of pressure sens- 
ing element has large signal output 
made possible by the frequency modula- 
tion or phase-modulation techniques 
used in the Photocon and Rutishauser 
systems.* Using high carrier frequen- 
cies (800 ke for Photocon and 25 mega- 
cycles for Rutishauser) it is possible to 
achieve high sensitivity with minute 
changes in capacity and also to attain a 
wide range of frequency response (11). 


‘Newer models now manufactured by 
Omega Instrument Co., Pasadena; in the 
strictest sense of the word, Photocon 
pickups use a quasi-F'M electronic system. 
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Fig. 1 Response of Photocon pickup to 
shock wave using nitrogen 


With an FM system the “number of 
counts per time period”’ can be adapted 
for pressure integration purposes. Al- 
though the capacitance principle is not, 
by itself, strongly susceptible to tem- 
perature effects, some capacitance gage 
systems do suffer seriously from thermal 
drift. This may be due to the thermo- 
elastic effects on the diaphragm or to 
changes in the dielectric constant at 
elevated temperatures. In many atomic 
energy applications where neutron 
irradiation effects are anticipated, 
chanzes in the dielectric constant due 
to nuclear activity in a number of 
capacitative and piezo-electric types of 
pickups can result in static readings 
which have no besring on _ pressure, 
sometimes creating errors of from 50 
to several hundred per cent (12). The 
transient response, however, is not likely 
to be seriously affected over a very short 
exposure time period. 

The capacitance type of pickup has 
been used rather extensively in Europe. 
Diehl and Visser at Delft, Holland have 
developed several such pickups of unique 
designs (13). The British have a num- 
ber of outstanding pickups, including 
those made by Southern Instrument 
Co. These pickups are used primarily 
with the FM carrier system. 

Electronic pulse techniques have also 
been used in sensing capacitance 
changes. Linlor, Kerns and Mark at 
Berkeley have reported a novel method 
which they have applied successfully in 
measuring liquid hydrogen pressure in 
the presence of pulsed magnetic fields 
(14). Pulses from a common source feed 
two cables, one of which is terminated by 
the pressure measuring capacitor in the 
pickup, and the other by an adjustable 
reference capacitor. These pulses are 
reflected from the terminations in mag- 
nitudes corresponding to their respective 


| 
| 


mite 


WITH H,O 


SWEEP: 10 CM 


Fig. 2. Response of Li-Liu pickup to 
shock wave using nitrogen 


impedance values, and their difference 
is amplified and displayed as voltage 

output by suitable circuitry. The ad- 

vantage of this method is that distrib- 

uted capacitance of the cables does not 

affect the sensitivity of the system. 

Cables of up to several thousand feet in 

length can, therefore, be used. 

Another new and significant develop- 
ment pertaining to capacitance-type 
pressure measurement devices is the 
T-42 ionization transducer credited to 
Lion of MIT and Switzerland. The 
system makes use of the gaseous elec- 
tronic principle of migration of space 
charges. It consists of a gas ionization 
tube exposed to a radio frequency field. 
The electrodes of the ionization tubes 
are connected to both a reference capaci- 
tance Cres and the capacitance C; to be | 
measured. When there is change in 
C,, a relatively large d-c signal is de-— 
veloped so that 


C set 


where Cre:, C; and K have compatible | 
wr and K isa constant. Since C,; = | 
(A/d), the output can be made essen- 
tially a linear function of the change of 
distance d (for small changes of C;). 
Piezo-electrie gages have been used 
for dynamic pressure measurements for 
many years. Some years ago the 
“Princeton Unit’ was developed for 
shock wave study and cther research _ 
physics. Bleakney and Arons have given — 
an authoritative account of this type of ; 
pickups in their papers (15). In Ger- _ 
many many developments in piezo-elec- 
tric pressure transducers were achieved — 
in earlier years. Due to their low 
curie points and pyroelectric sensitivity, 
the piezo-electric pickups are not suit- 


(Continued on page 128) 
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oh the Validity of Continuum Theory fo Satellite and 
Flight Problems at High Altitudes 


A study is wanilie of the validity of continuum flow theory 
applied to hypersonic flight in a rarefied gas. Order of 
magnitude estimates are made to establish the flight con- 
ditions for which a continuum analysis should be applica- 
ble in the stagnation region of a blunt body. Differential 
equations are presented for the flow in the stagnation re- 
gion for the flight regime where neither boundary layer 
theory nor free molecule theory is valid. The results of 
the study are applied to the problem of re-entry of a satel- 
lite and it is concluded that a continuum analysis, with no 


slip at the body surface, is valid for the flight conditions 


where heating is important. a7" 


SIEN (1)* and others (see, e.g., Siegel (2)) have proposed 

a division of fluid mechanics into various regimes ac- 
cording to the degree of rarefaction as measured by the cor- 
responding range of values of the Knudsen number, defined 
as the ratio of the mean free path \ to a characteristic flow 
length. Tsien utilized as a flow length either a characteristic 
body dimension, which we will denote by R,, or the boundary 
layer thickness. Of the various realms into which he and 
others have categorized fluid mechanics only two of these are 
reasonably well defined: They are the limiting cases of or- 
dinary continuum gas dynamics where the density is high 
enough so that intermolecular collisions dominate (A/R, < 1), 
and free molecule flow where the gas is sufficiently rarefied so 
that collisions with the boundaries dominate (A/R, > 1). In 
between these two limiting regimes there is of course a wide 
class of flows, although it is generally accepted (see, e.g., p. 381 
of (3)) that for \/R, > 1 a gas flow may properly be called 
rarefied. 

In the literature on low density flows (see, e.g., Schaaf and 
Chambré (4)) it is usually stated that continuum theory 
breaks down for a rarefied gas flow and that flow problems 
in this regime require treatment from the point of view of 
kinetic theory. However, Liepmann and Roshko (3) have 
pointed out that such a statement is much too strong and very 
often results from a comparison of a rarefied gas flow experi- 
ment with a solution of the Navier-Stokes equations which 
does not apply to the problem. For example, we would not 
expect the results of ordinary boundary layer theory, which are 
valid for high Reynolds numbers, to apply to a rarefied gas 
flow which is usually a low Reynolds number flow. It is 
for this reason that the use of the ratio of the mean free 
path to the boundary Jayer thickness is not completely 
satisfactory for defining the regimes of high Mach number 
rarefied gas flows. Furthermore, regardless of the character- 
istic length utilized it has been usual to compare this length 
with the mean free path in the free stream, \,,. However, 
this too is not satisfactory since for hypersonic flight the com- 
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pression of the gas across the shock wave and at the body 
increases the local density appreciably, so that the mean free 
path accordingly decreases in comparison with its free-stream 
value. 


Analysis of Flow Parameters 


In order to clarify more specifically the differences in mean 
free path we note, since \,, ~ v,,/@q (vq is the kinematic vis- 
cosity and a. the speed of sound), that the Knudsen number 
can be expressed in terms of the free-stream Mach and 
Reynolds number as 


R, Re, 


On the other hand, the mean free path behind a shock w ave is 


giv en by 
8 


But a ~ 7”, and at high temperatures the viscosity p ~ T'”; 
hence A, ~ €A.,, Where € = p,,/p, is the density ratio across the 
shock wave. Therefore behind a shock wave 


Now for air behaving as an ideal gas at hypersonic speeds 
e — i, while for actual atmospheric hypersonic flight be- 
cause of the activation of more degrees of freedom (vibration, 
dissociation) € ~ ;4 to 3. Of course for a weak shock wave 
e— 1, and Equation [2] reduces to the usual comparison given 
by Equation [1]. 

Carrying our arguments one step further, let us determine 
the mean free path at the body for the prac tical hypersonic 
flight case where the surface temperature is of the order of the 
free-stream temperature.4 Under this condition the mean 
free path at the body A, ~ (p../p~,)Aw- If, for simplicity, we 
restrict our considerations to the stagnation point region 
then p,/p.. ~ p.U?/p.. ~ M,,?, and 


M1 
R, M.Req 


We may now inquire as to what our previous results imply 
in terms of the actual flow field, say in the neighborhood of the 
stagnation region of a convex blunt body at hypersonic speeds 
(see Fig. 1). Assuming a continuum flow, from elementary 
continuity considerations we can determine the thickness A 
of the shock layer at the stagnation point, that is, the dis- 
tance between the body and the outer edge of the shock wave. 
The result may be expressed as A/R, ~ p,./p where R, refers 


*The surface temperature will of course be greater than the 
free-stream temperature. However, unless the surface tempera- 
ture is considerably in excess of 2000 F, the effect on mean free 
path at the body surface will not be too significant. For ex- 
ample, at a surface temperature of 2000 F, , will be about five 
times larger than given by Equation [3]. 
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specifically to the body radius of curvature at the stagnation 
point and where f denotes a mean density which is equal in 
magnitude to the density corresponding to sonic conditions. 
Since p*/po ~ 0(1) (the asterisk denotes sonic conditions and 
the subscript 0 stagnation conditions) it follows that 


On the other hand, at hypersonic speeds without considering 
relaxation effects, the shock thickness is given approximately 
by 6, ~ u*/a*p*. Hence for up ~ T'”* we have 


(Fe) (ANE) 


But (po/p*) 0(1) and € =~ Pa / Po; hence 


Since the radius of curvature of the body is of the order of 
the radius of curvature of the shock wave and since « < 1, 
then one thing we see from Equation [5] is that the shock 
thickness will be small compared to its radius of curvature so 
long as A, /R, < 0(1). Now still assuming the continuum 
estimates given by Equations [4, 5] to be valid, under what 
conditions will the shock layer thickness become of the same 
order as the shock wave thickness (see Fig. 1)? By equating 
Equations [4] and [5] this will occur for 


M 
~ —=.~ 1 (shock thickness ~ shock layer thickness). . [6] 


R, Re, 


In other words, for the condition of Equation [6] the shock 
wave and shock layer will have merged into one thin layer. 
But are continuum concepts valid for this condition? 

According to our definition of a rarefied gas we can say that 
behind the shock wave whenever \,/R, < 1 the flow behaves 
as a continuum.’ However, from Equation [2] this is true 
whenever 


1 
<— (continuum behind shock). . . . [7] 
€ 


Thus it is clear from the above condition and Equation [6] that 
under these circumstances the flow within the merged layer 
behaves as a continuum. 

Mention should also be made of the ‘ordinary boundary 
layer thickness” for the flight conditions we have been con- 


that at hypersonic speeds for a highly cooled boundary 
layer the boundary layer thickness in the neighborhood of the 
stagnation point is given roughly by 


| 


If we make the same reductions as previously, assuming the 
body temperature to be of the order of the free-stream tem- 
perature, we find that 


\? 
( Ry ) ( R, ) [8] 


The ratio of the usual boundary layer to shock layer thick- 
ness then becomes 


° The shock wave in this case is not defined sharply; however 
\, is taken to be of the order of magnitude of the mean free path 
corresponding to sonic conditions. 


sidering. For example, from the work in (5) it can be shown | 
> 
so that for the mixed flow condition \,/R, ~ 1/M,?. 
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Fig. 1 Schematic diagram of stagnation region for rarefied 
hypersonic flow 


and the inviscid flow region disappears when 


We may now interpret this result in terms of an approach to 
rarefied gas conditions from the usual continuum gas dy- 
namic regime, That is, considering the direction of increasing 
mean free path in the free stream we find that the ordinary 
“inviscid” shock layer has merged with the usual viscous 
boundary layer to give a fully “viscous” shock layer, before 
the shock wave begins to thicken appreciably. Thus, we still 
have a detached shock wave whose thickness is very small 
compared to the body radius in the region where we can no 
longer distinguish between an ordinary inviscid and boundary — 
layer behind the shock. 

In a similar manner we can examine the familiar continuum | 
gas dynamic “no slip” and “no temperature jump” boundary © 
conditions at the wall. From the usual estimates (see, e.g., 
(1, 3 or 4)) these effects will be negligible whenever A,/R, <1... 
But from Equation [3] we have 


Hence for M,, > 1 (say M,, > 10) the condition of ‘no slip— 
complete accommodation” will certainly be satisfied. We | 
further note that slip will not become an important factor — 
until] \,,/R, is greater than M 2/10, (that is, when \,/R, 
0.1). 

While plausible arguments have been made thus far regard- 
ing the validity of a continuum analysis, some remarks should 
be made about the breakdown of this analysis. First of all we — 
should discuss the boundary condition at infinity (ie., the 
outer edge of the shock wave). In private discussions with 
M. Krook, Lecturer on Astronomy, Harvard University, it 
has been pointed out that in this region the Navier-Stokes 
equations should be valid, regardless of the mean free path, 
because the gradients of fluid properties approach zero. Now 
a breakdown in the Navier-Stokes equations occurs when 
gradients exist such that any fluid property changes by an ap- 
preciable fraction in a distance equal to the local mean free | 
path. A region between the outer boundary of the shock and 
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the body surface will undoubtedly exist where the continuum shock wave and shock layer are indistinguishable in the 
equations will be invalid. Without a solution to the problem, neighborhood of a stagnation point on a smooth blunt convey 


one cannot say how large this region is. It is reasonable to body. The magnitude of the terms in the differential equa- 

state however that the errors involved in using a continuum tions can be related to the parameter A/R,, i.e., the ratio of 

analysis should be comparable to those in the calculation of a the shock layer thickness to body radius. This ratio is taken 

shock wave structure. Further it has been shown (7) that the to be small, as has been discussed, and in the present approxi- 

Navier-Stokes equations give a good approximation to a shock mation terms of order A/R, are neglected in comparison with 

wave structure at low Mach numbers; however no experi- those of order unity. 

mental results exist for hypersonic conditions. It is reasonable One other condition is imposed before arriving at the ap- 

to assume that the departures from a continuum flow will not proximate equations, namely that a boundary layer type phe- 

lead to first order effects up to the point where the shock and nomenon occurs in a thin shear layer near the body surface. In 

shock layer have merged. other words, in this shear layer the flow must be governed by 
Summing up our previous arguments we can conclude that the ordinary boundary layer equations. 

our continuum equations are still applicable at least up to, and The coordinate system is shown on Fig.1. Standard nota- 

possibly somewhat beyond, the point where the shock wave tion is used and for the sake of brevity the symbols are not 

and shock layer have all merged into one. Thus the Navier- defined. 

Stokes equations (not the usual boundary layer or shock wave After considerable reduction of the Navier-Stokes equa- 

equations) with the usual free-stream boundary conditions tions, the retention of the dominant terms leads to the follow- 

and “no slip-complete accommodation” surface conditions can ing equations 3 

be applied up to the region where the flow begins to behave 

like a free molecule flow. It follows then that the Navier- z-momentum ee = _ 


Stokes equations should be valid for the “transition” region 
between ordinary continuum hypersonic aerodynamics and puv 
free molecule flow, and the solution of these equations or ap- _ puu, + pouy + 10} 
propriately simplified forms may give results which go over Oo 
smoothly to the free molecule answer for slip and accommoda- 
tion coefficients of one. Such a conclusion also implies that, 
in contrast to previous thought, the accommodation co- 4 
efficient for energy transfer in the transition region is not an pu, = —p, + 3 
important consideration, since from our crude analysis it baad ~~ 
must be of the order of unity. 


y-momentum 


continuity = 
Simplified Form of Navier-Stokes Equations in 
Stagnation Point Region for Transition Regime Ff : 
An order of magnitude analysis has been carried out on . j = 0 plane flow 
the complete Navier-Stokes equations for the case where the ae => j = 1 axisymmetric flow 
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Flow regimes for the stagnation region of a blunt body at the satellite velocity 
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The boundary conditions consistent with the order of these 
equations, are 


[13] 


@ 


h = wall enthalpy = h,. 


Solutions of this system of equations should not only cover 
the condition where the boundary layer and shock wave are 
mixed but also the high Reynolds number situation where the 
shock, inviscid flow, and boundary layer are distinct regions. 

It can be shown, as expected, that these continuum equa- 
tions will give a surface heat transfer less than the “free 
molecule” result. Thus the surface heat transfer calculated 
from these equations and boundary conditions should pass 


continuously, with increasing mean free path, from the 


boundary layer theory toward the free molecule result. 


Transition Concepts Applied to Satellite Re-entry 


One of the more important problems of low density-high 
speed flight is associated with the re-entry of a satellite. It is 
therefore of interest to examine the flight regimes for a body 
traveling at the satellite velocity and indicate the significance 
of the results which have been presented. 

In order to furnish reference curves for the discussion to 
follow, we will consider the heat transfer rate to the stagnation 
point of a sphere traveling at the satellite velocity. Two 
limiting heat transfer curves are shown on Fig. 2 as a function 
of the ratio of the body radius to mean free path in the free 
stream. One curve gives the results of boundary layer theory 
(6) while the other represents the free molecule heat transfer 
for an accommodation coefficient of unity. In (6) a simple 
formula is given for the stagnation point heat transfer rate 
through a boundary layer. Fora body with a surface enthalpy 
small compared to the stagnation enthalpy the result is 


R, U 
qv.1.R, = 8.15 ) ) Btu/ft-sec.... . [16] 
sat 


where Usa: denotes the satellite velocity. The free molecule 
heat transfer rate for an accommodation coefficient of unity is 


R, \f U \* 
Ry = 572 ( Btu/ft-sec...... [17] 
No 


The intersection of the extrapolation of these two heat trans- 
fer curves occurs at a value of R,/A,, ~ 2. 

It is instructive now to indicate the various flow regimes, in 
terms of the ratio of the body radius to mean free path R,/d.,. 
It should be emphasized that these estimates are crude and 
represent orders of magnitude only. In carrying out these 
estimates we have taken M ,, equal to 25 and the density ratio 
equal to 0.1. The regimes indicated on Fig. 2 are defined as 
follows: 

1 The boundary layer thickness becomes equal to the 


given by 
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shock detachment distance eliminating the inviscid flow 
region. From Equations [4, 8] this occurs when 
~- 
ho WE 
2 The shock and boundary layer merge and become in- 
distinguishable. From Equation [6] 


3 The limit of continuum theory according to Equation 
[7] gives 
Rs 
20 


4 A surface “slip” condition is not expected to become 
important until the ratio of surface mean free path to body 
radius is the order of ;'5 (see, e.g., (4)). According to Equa- 
tion [3] for a surface temperature of the order of the free- 
stream temperature this condition occurs when s 


R, 10 l 
60 


For a surface temperature of 2000 F then 


Finally the regions of important re-entry heating for the two 


satellites (USSR Sputnik and U. 8. Vanguard) are indicated 
on Fig. 2. These regions were obtained from the results of 
(6), combined with additional calculations performed by 
James Hayes of the Avco Research Laboratory. It should be 
noted that for both satellites the critical heating occurs in a 
region where the continuum theory should apply and no 
“slip” occurs at the body surface. Of course noncontinuum 
theory will play an important role in determining the lifetime 
of a satellite in its orbit before re-entry. 
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The eccentricity and characteristic radius of a close 
satellite orbit decrease with time because of atmospheric 
drag. Variation of parameters and the Krylov-Bogoliuboff 
- approximation method are used to find simple approxi- 
mate expressions for the decay of eccentricity with radius, 
- the decay of radius with true anomaly and the growth of 
. = true anomaly with time. The special case of spiral orbits 
and the effect of asphericity and rotation of the earth’s 
=< ee: shell are discussed. 


Introduction 


PROBLEM of major concern in on engineering is 
the duration of the satellite on its orbit. Given an initial 
altitude and orbital eccentricity, how long a time or how 
- many orbital revolutions does it take the satellite to descend 
to any other assigned altitude? A subsidiary question of im- 
- portance to tracking, communication and various other 
- operational problems is the way in which eccentricity decays 
_ with altitude, time or revolution number. 
There are several reasons why the effect of air drag on satel- 
_lite orbits is rather difficult to analyze. One is the present 
gross uncertainty as to atmospheric density at extreme alti- 
tudes. Another is the uncertainty about a correct drag law 
and drag coefficient for a body passing through a rarefied gas 
_ with extremely high velocity. Not least is the difficulty of 
solving the nonlinear differential equations of motion, a 
- formidable computation problem even for a high speed digital 
computer if the descent is to be followed accurately through a 
Jarge number of revolutions. 
It is because of just these facts that an approximate 
analytic solution is useful if it is easy to apply and gives fair 
estimates of the true behavior. Parameter values can be 
changed easily as drag coefficients are better determined, 
changes in the density function can be encompassed quite 
readily, and the analytical form of the approximation provides 
considerable insight into the way in which air drag causes 
changes in the elements of the orbit. 

The effect of air drag on satellite orbits has been discussed 
by Petersen (4)? and others.* Reference (4) summarizes 
some of the geometric and kinematic properties of circular 
orbits, surveys the literature on the density function for the 
upper atmosphere,‘ discusses the value of the drag coefficient 
most appropriate in the domain of free molecular flow, and 
derives the equations of motion of a satellite subject to air 
drag, interpreting the acceleration terms which enter these 
equations. 

Three approximate radial decay results are given in (4). 
The first is based on energy considerations and has been used 
previously by Sanger and Ehricke, as Petersen observes. The 
second is a spiral approximate solution to the equations of 
motion, which is correctly recognized to correspond only to 
cases of extremely small eccentricity. The third is “a very 
approximate estimate” obtained by ignoring air drag near 
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apogee and making a rough estimate of the energy extracted 
by air drag acting near perigee. 

This paper is based on the simple drag model and equations 
of motion presented in (4). Those equations are solved by a 
method of variation of parameters, approximations being 
made in the solution of the variation equations. The «p- 
proximate results are valid for orbital eccentricity which is 
small, but much larger than that permitted by the spiral solu- 
tion of that reference. It is shown here that the decay of 
radius is the same as that given by the approximate treatment 
from energy considerations.* Moreover, there is an additional 
result on decay of eccentricity which does not seem to have 
been presented before. Numerical results based on these 
approximate solutions are included only to compare them 
with “exact” digital computer sclutions for a few numerical 
cases, thereby to show that the approximate method gives 
re¢ »sonably good answers. 

Certain additional topics are discussed briefly in the last 
section of this paper: The class of spiral orbits found by Peter- 
sen (4) and the effect of oblateness and rotation of the earth 
on the orbits. However, the treatment is descriptive rather 
than analytical. 


Approximate Solution 


Basic Equations 


With (4), it is assumed that the air drag force is always 
directed along the satellite negative velocity vector with a 
magnitude 


1 
9 


Here Cp is the so-called “drag coefficient,” A is the area of the 
projection of the satellite on the plane normal to its velocity 
vector, p is the ambient atmospheric density as a function of 
the altitude h above the surface of the earth, and V is the 
satellite speed. Actually, Cp generally is a function of density 
and velocity and A depends upon vehicle attitude. However, 
the assumption that they are both constant over most of the 
descent path probably is reasonable in this first analysis. It 
also is assumed that no forces act on the satellite other than 
those of air drag and the attraction of a spherical earth, and 
that the atmosphere has a standard stable density which is a 
function only of the altitude above the surface of the earth. 
Under these conditions, the satellite remains in a space- 
fixed plane, its position described by its radial distance r from 
the center of the earth and its angular advance 6 from an 
arbitrary reference line in the plane of the orbit. It is trivial 
to show that the equations of motion are 
4 K pA 


— rp? = — a V [2] 


5 The model is capable of considerable generalization without 
voiding the approach used here, although the details may become 
more tedious. 

* This paper utilizes a formal approach through the equations 
of motion which superficially is quite different from the approacl 

through conservation of energy. The reason why they give the 


- same approximate result will become clear later. 
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tional constant. It is convenient to introduce auxiliary di- 


mensionless variables 
—&=a/r = Ka/(r8)?...... 


where a@ is the equatorial radius of the earth, and to use B as 
the independent variable. (Primes denote differentiation with 
respect to this variable.) Also, define the parameter v by 


For the spherical earth, h = r — a and 
[6] 


The basic equations, [2, 3], beeome 


if there were no air drag, the satellite would follow the 
ellipse described by 


[1 + ecos (8 — Bp)]............. [9] 
d 
a 


where # is a characteristic radius of the orbit related to the 
radius of the perigee by R = Rp (1 + 6), and Bp is the value 
of 8 at the perigee. 

The method of solution is based on the requirement that the 
form of the solution in the presence of air drag be the same as 
that of Equations [9-11] when R, e« and Bp are considered to 
be slowly varying functions of 8. There is no loss of generality 
in choosing Bp = 0 when 6 = 0, the “initial instant.’’ Denote 
initial values of R and ¢ by R; and ¢;, respectively. 

Bearing in mind that R, « and Bp are to be treated as 
variables, differentiate Equation [11] and use the result to- 
gether with Equations [9, 10] in Equation [8]. The result is 


where 
G(R, 6, 8 — Bp) = 
V1 + + 2€ cos (8B — Bp) Rk 13 
[1 + ecos (B — \it+ecos(B— Br) “| 


Additional differential equations connect R, «, and Bp be- 
cause there must be consistency among Equations [9-11]. 
That is, the &’ found by differentiating Equation [9] must 
equal that given by Equation [10]; also, because of Equation 
[7], the €” found by differentiating Equation [10] must equal 
the difference » — & formed in using Equations [9, 11]. These 
equations of consistency are 


e’ cos (8 — Bp) + €Bp’ sin (8 — Bp) = 4 [1 + €cos (8 — Bp)] 


FEBRUARY 1958 


Here K is the product of the mass of the earth by the gravita-_ 


Fig. 1 §8-8* phase plane > 


from which 


le + cos (8 — Bp)]G(R, — Bp) ... [16] 


=— sin (8 — Bp)G(R, B — Bp)...... {17] 


Equations [12, 16, 17] form a basic set for the variation of 
orbital elements. To provide a guide to the approximations it 
is reasonable to make in solving these equations, it is useful 
to consider the topological nature of the solution trajectories 
in a certain phase plane. 

An important preliminary question concerns the magnitude 
of vG/e. For any altitude where the lifetime of the satellite 
is more than a few days, and for any conceivably “peor 
satellite area-mass ratio, the product vG is very small, prob- 
ably one or more orders of magnitude smaller than 10~. On 
the other hand, « can be expected to exceed 10~% or 1074. | 
Therefore, the magnitude of »G/e ordinarily is small enough 
that its square is negligible for the purposes of this approxi- 
mate analysis. This chapter is limited to the case for which — 
this is true. 

Equation [17] can be rewritten as 


* 
= =1+ G (R, B*) sin B* 
€a 


[18] 


where 8* = 8 — Bp. The solutions of this equation are to be _ 
examined as a family of trajectories in a 68*-plane, in which : 

8* and B are the variables along the axes of ordinates and © 
abscissas, respectively. The particular solution of interest is 
the one passing through (0, 0) in this phase plane. The 
several parts of Fig. 1 show the details of several portions of 
the phase plane which are discussed below in turn. 

Because it is supposed that vRG/ae < 1 initially, Equation 
[18] shows that the slope of the phase trajectory is initially 
unity in the 68*-plane. Because it passes through the origin, 
the sclution appears as shown in part A of Fig. 1. It follows 
closely along 8 = 6* for many periods, as long as vRG/ae is 
quite small, without significant wiggles. This actually should 
be a good description of the 8, 8* relationship over most of — 
the satellite lifetime, in view of the very small eccentricity _ 
which must be reached before it becomes comparable to vG. 
Physically, this result means that @p remains zero and that 
the perigee does not begin to rotate as long as the eccentricity 
is considerably greater than vG. (Note, however, that 
perigee rotation from the oblateness of the earth occurs in- 
dependently of any effect from air drag.) a 

It is shown later that vRG/ae grows with 8. When this 
parameter has ceased to be small, 6* does not increase as — 
rapidly as 8, as shown in part B of Fig. 1. It is this range in 
which an approximate solution is especially difficult to obtain. 
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Physically, the perigee begins to rotate positively with a 
gradually increasing speed. 

As vRG/ae approaches and exceeds unity, the situation be- 
comes simpler. Fig. 1,C, shows the behavior for this range 
of the parameter. The scale of B* issuch that (2n — (1/2))7is 
_ the last position where sin 6* = —1 before the occurrence of 
_ vRG/ae = 1. As this parameter approaches unity from be- 
low, the trajectories cross the sin 8* = —1 lines at successively 
lower angles, so that the trajectory covers a large 6 range for 


- each 27 increment of 8*. Finally, when vRG/ae = 1, two 
_ lines whose equations are 
gin Go* = [19] 


_ spring into existence at (2n + (3/2))m and gradually diverge 
as the parameter grows (Fig. 1,D). These are lines of zero 
“a rate of change of 8* with respect to 8. Because they diverge 
slowly with 8, they are essentially horizontal themselves in 
; b any small 8 neighborhood, and act as a horizontal barrier to 
the phase trajectory. 

_ It can be shown that this qualitative description is correct 
_ by considering solutions of Equations [12, 16, 17] in the 
_ neighborhood of these singular lines. However, detailed ex- 
amination of the case where vRG/ae is of the order of unity or 
___ greater is not justified except perhaps for very special applica- 
tions. The radial variation implied by the € corresponding to 
_ this range of values would be very small—the order of a few 
meters for h = 1000 km. Such a variation is far over- 
shadowed by eccentricity effects arising from the oblateness of 
earth. 

This fortifies the restriction imposed previously that only 
_ vG@/e considerably less than unity is to be considered. The 
fact that the behavior of the solution in the 6, 3* phase plane 
does not change much until the eccentricity drops far below 
the range of practical interest, assures that this limitation does 
not exclude any cases of real importance. 


Approximation of G(R, ¢«, 8*) 


: With d8*/dB8 = 1 holding very closely over the 8 range of 
interest, Equations [12, 16] can be rewritten 


R? 
[20] 


dB* a 
d 
(€ + cos B*)G(R, €, B*)........ [21] 


Successful approximate solution of these equations depends 
upon the linearization of the nonlinear function G(R, «, *). 
In simplifying the form of G, one is tempted to expand the 
density function in Equation [13] in a power series in €, drop- 
ping all terms in ¢ of order higher than the first in the resulting 
expression for G. It is true that the argument of the p func- 
tion can be approximated very closely when second order 
terms in € are dropped. However, dropping similar terms 
when p is expanded in a Taylor series is equivalent to repre- 
senting the density at all points of the orbit by the density at 
the perigee, plus a term for the radial rate of change of 
density. This procedure is not very accurate for much of the 
eccentricity range of interest because of the distinctly non- 
linear behavior of density with altitude. Alternatively, one 
can take the point of view that eR, half of the maximum total 
variation of altitude as the orbit is traversed, is not very small 
relative to the altitude h = R — a, so that it should not be 
used as a small expansion parameter. For example, an ec- 
centricity of « = 10~ is not large, yet it represents an 
altitude variation of almost +70 km ath = 500 km. To fail 
to handle the expansion of p with care may be to sacrifice 
immediately the possibility of using the results to obtain 
good numerical information. 

A key step is now taken. It is assumed that in a restricted 
neighborhood of any particular value of &, the density function 
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can be approximated by 


where, in general, a and 6 vary slowly with the & (altitude) 
level. It can be found from Reference (1) that Equation 
[22] is an excellent representation, with a and 6 essentially 
constant above R — a = 300 km, approximately. This 
analysis treats a and 6 as constant over the altitude ranges of 
interest. 

Under these circumstances, one can show that 


R 
at = p(h) exp {xe cos B/(1 + € cos 


where 


G(R, €, B*) = p(h)eX* 8*[1 — € cos B* + O(e*)].... [26] 


and to within terms of the order of e?, Equations [20, 21] take 
the forms 


2 
p(h)eX* 8* (1 — € cos B*)...... [27] 
a 


de vR 
dB* a 


Because the density is very small and dR/d8* and de/d8* 
are correspondingly small, it suggests the applicability of the 
averaging method of Krylov and Bogolinboff [2]. Approxi- 
mate differential equations for the behavior of the elements 
are obtained by averaging the original differential equations 
over a full period in 6*, removing all of the trigonometric 
terms in the right-hand sides of Equations [27, 28]. We use 
the fact (5) that 


p(h)e** 8* [cos B* + €(1 — cos B*)].. . [28] 


where J,,(z) is the usual notation for the Bessel function of the 
first kind for imaginary argument. 
Upon averaging, Equations [27, 28] become 


pth) [To(ke) — el, (ke) + O(e*)]...... [30] 


pth) {Ts (Ke) + €[To(Ke) — Ii(ke)] + O(e*)}... 

dp a 

Note that 6 has been used interchangeably with 6* in view 
of previous results. The importance of this representation 
becomes evident when it is recognized that x can be of the 
order of 100 or larger. In this case, for eccentricities of € = 
0.01 or larger, the power series expressions corresponding to 
the functions Jy and J; could not be truncated with ¢? or é 
terms without disaster. Yet that is precisely what would 
have been done in expanding G completely into a power series 
ine. 

It is important in the sequel to note that the dominant term 
in the right-hand side of Equation [30] is Io(xe), while the 
dominant terms in Equation [81] are I;(xe) + elo(xe). The 
other terms are ignored in our approximate solution. 


Solutions 


Next, consider the solutions of Equations [30, 31]. Rather 
than trying to find e directly as a function of 8 through 
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Equation [31], it is convenient to divide this equation 

by Equation [30] to obtain 
de _€ 
dR 


RI(ke) 


which, in view of the definition of x, can be rearranged as 
d\n (eV p/R) F (ke) _k [33] 
This can be rewritten as an integral equation = 
p(Ri—a) Le Tome) 2) 
which can be used as the basis for an iterative solution. 
Note first that for xe << 1, J;(xe) ~ xe/2 and Jo(xe) ~ 1, so 
that the integral on the right-hand side of Equation [34] 


vanishes to within second order terms in small quantities. 
For this case, then 


eR _|p(R;i — a) 
R) =- 
Ri \ 


— a) 
which may be expected to be very good for ¢ of the order of 
0.001 or smaller, with « of the order of 100. 

Now let the result of Equation [35] be used in the right- 
hand side of Equation [34] as the first step of an iteration 
process. The product xe is of the form 


= aae:V pi 1 
Ke = — 


A 


eR; 
R 


[36] 
Ri Vp 
so that 
1 d(ke) kK 
[37] 


xe 


and by introducing xe as the variable of integration into the 
integral in Equation [34] this equation can be written 


eR; p(k a) e(xiei) 
= exp [@(kie:) — O(ke)} = [38] 
Here 
Ke 9 1 
0 x I(x) 
and 


R; — a)]* dR 

The functions ¢ and e® are tabulated in Table 1 for a range 
of xe which should span all values of interest. Actually, even 
the very great eccentricity mentioned for Vanguard (apogee 
altitude 1000 miles, perigee altitude 100 miles) gives a value of 
xe no greater than about 12. 

If desired, of course, the new value of ¢ obtained from Equa- 
tion [88] can be used in the right-hand side of Equation [34] 
to carry the iteration a step further. However, the inherent 
limitations on the method and uncertainties as to parameter 
values probably are such that a higher order iteration cannot 
be justified. 

In summary, e(R) is calculated by using Equation 40 to find 
xe(R); then Table 1 to find e%*©/e%*e), From this result, 
Equation [38] permits a direct determination of «(R). It is 
noteworthy that the result does not depend on »y, which is to 
say it is independent of the satellite area/mass ratio and the 
drag coefficient (subject only to the constancy of these quan- 
tities assumed previously). Thus it can be found once and for 
all from a knowledge of the atmospheric density function 
alone, given the initial altitude and eccentricity. 
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Table 1 Functions ¢ and e? 

Ke (Ke) eh (xe) Ke Ke) 

0 0 1.000 2.0 0.1926 1.212 

0.1 0.0010 1.001 3.0 0.3446 1.411 

0.2 0.0032 1.003 4.0 0.4921 1.636 

0.3 0.0063 1.006 5.0 0.6276 1.873 

0.4 0.0107 1.011 

0.5 0.0162 1.016 6.0 0.750! otis 

0.6 0.0227 1.023 7.0 0.8606 2.365 ; 
0.7 0.0303 1.031 8.0 0.9611 2.615 . 
0.8 0.0389 1.040 9.0 1.0531 2.867 | 
0.9 0.0483 1.049 10.0 1.1376 3.119 : 

1.0 0.0585 1.060 11.0 1.2156 3.372 | 
1.1 0.0695 1.072 12.0 1.2881 3.626 | 
1.2 0.0813 1.085 13.0 1.3561 3.881 
1.3. 0.0937 1.098 14.0 1.4201 4.1388" | 
1.4 0.1067 1.113 15.0 1.4801 4.393 
1.5 0.1201 1.128 16.0 1.5366 4.649 

1.6 0.1340 1.1438 

ef 0. 1482 1.160 Beyond xe = 16, the asymp- 

1.8 0. 1628 ce ig totic relation ¢ ~ In re + 
1.9 0 177 1.194 const is closely valid 


We return now to Equation [380] for R as a function of 8B, r 
keeping only the dominant term on the right-hand side. Be- 
cause K = xe is the function of R already found above, a — 
direct integration gives 


Ri 
vB = of 
R 


This integral is easy to evaluate numerically. 
If xe << 1,1)(K) ~ 1 and the relative change in altitude per 
revolution becomes closely equal to 


Ah 2nvR%p(h) 
— }/revolution = —— 
h ah 


dR 
R%p(R — 


which has the general form of Petersen’s Equation [40]. It — 
should be expected that this be the same result as obtained | 
from a conservation of energy argument, for it is well known | 
((3), p. 238) that the Krylov-Bogolinboff averaging method is — 
equivalent to a conservation of energy approach. However, 
the present result is a generalization which encompasses the 
effect of « on the radial decay. 

After R(8) [strictly, 8(R)] is found from Equation [41], | 
the time behavior of the orbit is obtained by returning to the 
definition of 7 in Equation [4] and its behavior in terms of R } 
in Equation [11]. From these 


From Equations [38, 41, 43] all the essential behavoir of an 
orbit having eccentricity v»RGa~! < <1 can be inferred 
semiquantitatively. Note that in both Equations [88, 42] 
it is natural to use R as the independent variable. 


Evaluation 


It is imports int to try to discover how close the ap Pee roximate | 
results given by Equations [38, 41] and d8*/d6 = 1 are rela- 
tive to the exact results that would be obtained by solving — 
Equations [12, 16, 17]. A comparison of results for two cases 
have been made which indicates that the approximate rela- 


Wes 
[23] 
9 
] take Pee 
27 
— d8 [43] 
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‘y) Equations [12, 16, 17] were solved as a coupled set numeri- 
ot cally® to give an “exact” solution for a few parameter values. 

The value® of » was taken as 5 X 108 cm?/gm and p(h) was 
taken as the upper bound of the estimates of (1). Initial 
conditions 6; = 0 and R; — a = 400 km, the latter chosen be- 
7 cause it is low enough that changes can soon be seen in the de- 

‘ _ pendent variables, yet high enough that the total lifetime 
occupies a reasonably large number of revolutions (15,000 to 
E 20,000 radians), thereby setting a rather stiff task for the 
approximate method. The two values = 0.001 and = 
0.01 were considered. 

First, it was found that 6* remained equal to 8 to within 
_ Jess than one part in one million over the 6-range considered, 
thereby fortifying the assumption 6* = £. 
Fig. 2 shows the results for ¢; = 0.001, both «/e; and 8 
plotted as functions of R. The solid curve is the exact solu- 

tion for B(R), the dashed curve is the exact solution for e(R)/e;, 

and the small open circles are the values obtained by the ap- 

proximate method of this paper. As might be expected for 
such a small ¢;, the approximate method gives excellent re- 
sults, departing only a few per cent from the true values for 
the smaller altitudes where the satellite is about to fall to 
earth within a few revolutions. 
Fig. 3 is the same as Fig. 2 except that the initial eccentricity 
oe is much larger, e¢; = 0.01. The agreement between approxi- 
foe & — mate and exact solutions for e(#) is quite good, as shown by 
the open circles in relation to the dashed curve. The crosses 
near these dots show the approximate result that would have 
_ been obtained ignoring the departure from unity of the right- 
hand side of Equation [38]. 

Unfortunately, the agreement for 6(R) is not as good. 
_ However, the approximate method is not so far in error as to 
be completely useless. It is of considerable interest to note 

the importance of retaining the Jo(xe) term in Equation [41]. 
Without it, the approximate 8(R) ee would be that 
shown by the open — in Fig. 3 


Other Topics 


Spiral Orbits 


The analytical development above is predicated on the 
existence of an elliptical first approximation to the motion. 
_ It ean happen that such an elliptical approximation does not 
exist, as can be seen physically as follows. 

Suppose that at some instant the satellite has a velocity 
vector with a small component directed toward the earth. 
Normally, this condition represents the return from an 

- apogee position. The satellite gathers speed as it continues 
7 - inward toward the earth and eventually reaches perigee with 
a higher total velocity than it had at the prev ious apogee. 
Se erent if the inward velocity component is very small, 
ee orresponding to an orbit which has very little eccentricity, 
_ the velocity gain between apogee and perigee is very small. 
_ The satellite falls such a short distance toward the earth in 
_ reaching perigee that its velocity increases very slowly. In 
the presence of air drag, the velocity loss from drag may ex- 
- ceed the normal gain from the fall toward the earth, and in 
_ this case the perigee is never reached. The path is a spiral 
rather than an ellipse with decreasing size. 

To find the change in radius for a spiral orbit, it is most 

convenient to return to Equations [7, 8] and seek their solu- 


: 7It is emphasized that these calculations do not purport to 
_ describe the actual behavior of an orbit because of the uncertainty 
as to a correct value of » and a correct density function. How- 
ever, using approximate values in these cases does permit results 
— to be compared which are calculated in two different ways. 

8’ By IBM-701 computer. 
* This large value of » should give no cause for alarm, for vp in 
_ the neighborhood of 400 km is of the order of 5 X 10-7. This » 


value is conservatively large, and more realistic, smaller values 
will further improve the approximation method. 


tions may be quite adequate for engineering purposes.’ 3 : tion by iteration. As initial values, use £; = a/R, and n; = 


a/R;, but leave open for the moment the choice of &;’.. Denote 
the order of the iteration by superscripts in parentheses and 
begin with 


= &; &'(B) = 0; = H....... [44] 
From Equation [8] 
7(B) = — [yp(R; — a)B}........... [45] 
This used in Equation [7] gives iel 


after choosing the coefficients of the solution of the comple- 
mentary equation equal to zero on the ground that any other 
choice allows an elliptical approximation to &(8), contrary to 
assumption. 

From this result it is seen that no freedom remains to 
choose £,;’; one simply must accept the value 


= 


implied by Equation [46]. 
In terms of the original physical variables, this spiral solu- 
tion is 


r(B) = e(Ri — 


with 6 closely equal to 


This corresponds to the spiral solution given as Petersen’s 
Equation [57], although it has been obtained by a somewhat 
different argument. 

Because this class of spiral orbits corresponds to such an ex- 
ceedingly small radial velocity, the transition from an elliptic 
to a spiral orbit must occur for a minescule value of eccentric- 
ity. Therefore, spiral orbits are of interest only if the initial 
eccentricity itself is very small or if the radial distance be- 
comes so small that the satellite is ready to fall precipitously 
toward the earth. Such small initial eccentricities are ex- 
tremely improbable, especially in view of the equivalent eccen- 
tricity which is induced by the oblateness of the earth (a topic 
treated by the author elsewhere). Therefore, one concludes 
that the class of spiral orbits is without practical importance. 


Effect of Aspherical Atmospheric Shell 


The analysis above has presumed a spherical earth and a 
spherical atmospheric shell, characterized by our representing 
altitude as simply h = r — a. Actually the earth departs 
significantly from a sphere, and one may expect the at- 
mospheric shell to be similarly oblate. Let us consider what 
effect it might have on the analysis to take account of this 
property. 

The oblateness of the earth really has two important effects 
on the air drag problem. Even neglecting the asphericity of 
the atmospheric shell, oblateness is important because it 
generally contributes to the eccentricity of the orbit. I have 
shown elsewhere that the effective geometric eccentricity of 
the orbit is compounded of the dynamical eccentricity used 
in the analysis above and terms proportional to the ellipticity of 
the earth. If this fact is ignored in establishing the satellite 
orbit, then the geometric eccentricity is nonzero even if one 
succeeds in achieving a dynamical eccentricity which just 
corresponds to a circular orbit. 

The second effect arises through the oblateness of the at- 
mospheric shell, and is the one considered below. To get a 
simple picture of the effect, assume that the density is con- 
stant over shells elie the same mene as the surface of 
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the earth. In this case, one should use in the density function 
an effective altitude given by 


h 
1+ e€cos8 + 2e* sin? 
where 
ay 
== (a*/6* — 1) gin? =... 
b = polar radius of the earth 
y = dihedral angle between the plane of the orbit and the 
earth’s equatorial plane, with 0 < y (y > 
corresponds to retrograde orbits) 
A = geocentric latitude of the satellite 


If Bp denotes the angle from the ascending node of the orbit to 
the perigee (at which 8 = 0), positive in the direction of 
satellite motion, then 


sin = sin y sin (@ + 
Therefore, the analog of Equation [23] for the case at hand 
is 


V/1+ 2e* 1+ €cosB 


O(e?, e€*, €*?) 


p(R — a) exp E cos B — sin? y sin?(8 + | + 
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Fig.3 Change of eccentricity and angular advance with altitude 
TUF 


= p(R (ae* ain? Y/2R exp E cos B + 
* 
sin? y cos 2(8 + a» [53] 
Also, to the same order of approximation 
G(R, €, B*) = [p(R — a)e~ y)/2R (1 — € cos B)] X 


* 
exp {nos B+ sin? cos 28p cos 2a x 


kae* ¢ 
[: OR sin? sin 28p sin 28 + | 


Here, the last square bracket is a convergent power scrics 
representation of 


* 
exp {ae sin? y sin 28p sin 2at 


In the Krylov and Bogolinboff averaging procedure, all of the 
terms involving e* and its powers in this square bracket will 
vanish because these terms are odd functions in 8, so they can 
be discarded at once and G can be represented simply as 


G = p(R — ae~ (2e* sin’ y)/2R (1 — € cos B) X 
kae* 
exp 4 xe cos B + OR sin? y cos 26p cos 2B ¢..... [55] 


The author is not aware of any result for this expression 
corresponding to Equation [29] which can be used to repre- 
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sent the averaged G-functions in terms of tabulated functions 
of xe and xe*. However, the fact that e* ~ 0.0033 sin? y 
implies that, especially for orbits which are not too close to 


meridian orbits, the expansion 


ef 00828 1+ €cos2B............... [56] 
57] 


is likely to be a tolerably good representation. Thus, the 
analogs of Equations [30, 31] are 


vR 
= —— - a)e ~ (2¢* sin? y)/2R x 
a 
= 


Which terms are dominant in Equations [58 and 59] de- 
pends on the relative magnitudes of e andé. Consider two 
special cases: € < Llandxe > €. 


whence 
R 
exp sin? y cos m aR ¢ 


Bean the result for the spherical earth and small eccentricity 
is just modified by the exponential function of R which can be 

- evaluated numerically without much difficulty. 

The magnitude of the exponential function is generally 
close to unity, but there is one case in which the change is 
significant. Consider Equation [61] for values of R extremely 
close to R;, so that 


QR; 4R2 


Suppose nowthat 


aa 
e* sin? y cos 


| 
| 
} 


as it may for lower altitudes (where p’/p is large) and orbits 
for which sin y and |cos 28p| are sufficiently close to unity. 
If the perigee is near the equator (8p near 0 or x), the co- 
efficient of (R; — R)/R is positive and the eccentricity decays 
as altitude decreases. (This is the case in which the perigee 
occurs in denser regions of the atmosphere, in as much as a 
fixed & corresponds to lower altitude near the equator.) On 
the other hand, if the perigee is above 45° N latitude or below 
45° S latitude, so that cos 28p < 0, the coefficient of (R; — 
R)/R; is negative and the eccentricity increases instead of de- 
creases. This is a new qualitative effect which is not obtained 
when the earth is considered spherical. 

If € < «xe, the dominant terms in Equation [58, 59] are the 


+ €[Zo(ke) — + + Is(ke)]}... . [59] 


2 
(1 + sint-y | X 


same as in the case treated earlier in which ellipticity is 
ignored. 

These considerations indicate that as e decreases, the effect 
of the oblateness of the atmospheric shell increases, but that it 
does not change the relative values of ¢ much even for small . 
However, in some cases an additional eccentricity is induced 
rather than the eccentricity monotonically decaying with 


altitude. 
Effect of Earth’s Rotation 


It has been assumed previously that the drag force acts 
toward the negative velocity vector of the satellite, con- 
sidered relative to inertial space. Actually, one should con- 
sider the velocity relative to the air mass producing the drag. 


_ If the atmosphere at satellite altitude partakes of the earth's 


rotation, another component of drag force exists normal 
to the plane of the orbit. 

Relative to a rotating atmosphere, the component of drag 
in the plane of the orbit can be treated exactly as in the second 
section, except that the force has a small oscillatory com- 
ponent as the satellite traverses its orbit. For y < 2/2, thie 
force is slightly less at the extreme excursions of latitude than 
near the equator, for the relative velocity component in the 
plane of the orbit is slightly smaller at these points. Con- 
versely, for y > 2/2, the drag force is slightly greater at thie 
extreme latitudes. 

A similar oscillatory force exists because of the ellipticity of 
the atmospheric shell, as discussed in the previous section. 
Over an oblate earth, the density at a constant radial dis- 
tance is less at high latitudes than near the equator. It 
follows that the effect in the orbital plane of the earth’s rota- 
tion may be analyzed in a manner very similar to the analysis 
of the earth’s oblateness effect. Moreover, for retrograde 
orbits these effects tend to compensate. The effect of a 
rotating atmosphere can be expected to be small, for the 
maximum difference between relative and absolute velocities 
is only about 5 per cent of the absolute velocity if the at- 
mospheric shell rotates as a rigid body with the angular speed 
of the earth. This almost certainly is not the case, and the 
atmosphere at altitudes of interest is likely to have a much 
lower velocity than the earth’s surface. 

The component of drag force normal to the plane of the 
orbit is another matter. It cannot be treated within the 
analytical framework of this paper: One must return to a 
complete set of equations of motion in which the possibility 
of motion out of the plane of the orbit isadmitted. Although 
the analysis of this “cross-drag’’ has not been attempted, it is 
expected that this effect will result in some rotation of the 
plane of the orbit and that it may have some influence on the 
orbital inclination. However, it is conjectured that any such 
effect is likely to be small enough to ignore in most applica- 
tions. 
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Investigation of The Flow 


Over Simple Two-Dimensional and Axial Symmetric Bodies 


Hypersonic Speeds 


on 


I. E. VAS,’ S. M. BOGDONOFF® and A. G. HAMMITT: 


fe 

4 study of the inviscid shock shape and pressure dis- 
tribution about thick two-dimensional and axial sym- 
metric bodies was conducted in the Princeton University 
Helium Hypersonic Tunnel. Details were obtained for 
flat and semicircular leading edges for Mach numbers be- 
tween 11 and 19. From the Schlieren pictures, the shock 
shape, detachment distance and sonic point on the shock 
were found. Similar pressure distributions were obtained 
for corresponding two-dimensional and axial symmetric 

bodies. 


Nomenclature 
a,b = defined by y/t = a(x/t)? 


maz 
= aM, 
= cylinder leading edge thickness, inches 
h = density ratio across a normal shock 
M = Mach number 
p = static pressure, psia 
R = radius of curvature 
Re = Reynolds number 
= surface distance measured from the center line stagna- 
tion point on the body 
= plate leading edge thickness, inches 
y = shock coordinates measured from the shock apex, inches 
) = ratio of specific heats 
5 = detachment distance, inches 
Ap =p-p 
M = angle between the tangent to the body surface and the 
free stream direction 
p = density 
Subscripts 
Q = stagnation conditions 
| = free stream conditions = 
b = body 
s = shock 
t = leading edge thickness (applied to Reynolds number) 
T = conditions behind a normal shock 
Introduction 


T HIGH Mach numbers, dissociation and ionization of 
the air may have a considerable effect upon the aerody- 


Received June 20, 1957. 

! The present study is part of a program of research on hyper- 
sonic flow under Contract No. AF33 (616)-2547 to Wright Air 
Development Center, U.S. Air Force. 

2 Research Associate. 

’ Professor; Head, Gas Dynamics Laboratory. 

* Senior Researe h Associ late. 


FEBRUARY 


Gas Dynamics Laboratory, Forrestal Research Center, Princeton University, Princeton, N. J. 


4 


namic characteristics of a body. As aerodynamic heating 
may cause serious problems at these speeds, the removal of 
heat near the front of the body must be considered. For this 
purpose, a blunt body would present certain aerodynamic ad- 
vantages. The study of the conditions on this body at simu- 
lated flight conditions at these Mach numbers would be ex- 
tremely difficult. However, the perfect gas aerodynamics of 
such bodies may be investigated as part of the over-all problem. 
Helium is the ideal gas to be used as its liquefaction tempera- 
ture is about 23 R at 0.0015 atm and its specific heat is ap- 
proximately constant over a wide range of pressures and tem- 
peratures. 

Newton has postulated that when discrete particles strike 
a surface the normal momentum component is eliminated 
and the tangential momentum component remains un- 
changed. The loss in the normal momentum causes an in- 
crease in the pressure. The value of the pressure coefficient 
predicted is 

p = 2 sin? 6, 


where 6, is the angle between the free stream flow direction and 
the body surface. This analysis holds only for y —~ 1 and 
M — o for which case the shock lies on the body surface. 
For y ~ 1 and M — o the pressure coefficient at the stagna- 
tion point, Cymax, isequalto2. For y = 5/3 and M = 15 the 
maximum pressure coefficient is 1.75. By accounting for the 
strong leading edge shock, one has 


a = sin? 6, 

Cpmax 
which can be used to predict the pressures on blunt bodies 
(1-4).® For blunt leading edges, the inviscid pressure on a 


surface parallel to the free stream has been predicted in (5-7) 
for two-dimensional and axial symmetric bodies. There seems 
to be some discrepancy between the values obtained from ex- 
periment and the predicted values. 

The shape of the shock, distance of the shock from the 
stagnation point and the sonic point on the shock have also 
received considerable attention (8-15). Reasonable agree- 
ment is obtained between the predicted values of the de- 
tachment distance between theory and experiment (16-19). 
Several experimental investigations have been made to deter- 
mine the pressure distribution on hemisphere cylinders at 
supersonic speeds and low hypersonic speeds (2, 3, 4, 20). 
experimental data are available of the pressures on the surface 
parallel to the free stream direction of a blunt flat plate or 
cylinder with the leading edge perpendicular to the free 
stream direction. There also seems to be a lack of such in- 
formation for the case of the flat plate with a hemicylindrical 
leading edge at hypersonic speeds, although some data (18) 
are presented for this type leading edge at a Mach number of 


. 


5 Numbers in parentheses indicate References at end of paper. 
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7 _ In some previous work performed at the Gas Dynamics 
Laboratory of Princeton University (21-23), the flow about 
flat plates was studied with the leading edge thickness vary- 
ing from 0.13 X 10~*in. to 44 X 10-3 in.; however, no details 
were sought near the leading edge. It is the purpose of the 
present paper to obtain the details near the leading edge of 
both two-dimensional and axial symmetric blunt bodies with 
square and rounded leading edges. The shape of the shock 
about the body, detachment distance, and sonic point are 
measured and compared to predicted values. The pressure 
distributions on these representative models are also pre- 
sented. 


Fig. 1 The pressure models 


a) Flat leading edge plate, = 13.7 


b) Hemicylindrical leading edge plate, 17 = 13.8 


(d) Hemisphere leading edge cylinder, = 12.7 i 
Fig. 2 Schlieren photographs of the bodies 


Experimental Facilities 


The experimental work was carried out in the Princeton 
University Helium Hypersonic Tunnel. The details of th: 
tunnel have been given previously (24), so only a résumé 
shall be presented here. An axially symmetric nozzle with an 
initial taper of 53 deg, joined to a 3}-in.-diam cylindrical sec- 
tion was used. The sides of the nozzle were milled parallel to 
each other and windows installed. Disturbances caused by 
the window-nozzle junction do not affect the flow in the test 
region of the tunnel, as the wave angle formed at these high 
Mach numbers is small. By changing the throat diameter. 
it was possible to change the Mach number. The Mach 
number range was 11 to 19 with stagnation pressures from 700) 
to 1500 psia. Because of the conical design, a Mach numbe: 
gradient existed along the center line of the nozzle. With th 
larger throat diameter, the Mach number range was 11 to 14 
with a gradient of 0.5 Mach number per inch; and with th: 
smaller throat diameter, the Mach number varied from 16 to 
19, with a gradient of 0.7 Mach number per inch. An attempt 
to correct for the pressure gradient in the tunnel was mad 
and is explained later. 


= Test Program 2 

A test program was set up to study in some detail the flow 
about thick two-dimensional and axial symmetric bodies with 
square and round leading edges. A single representativ: 
model was used for each case with pressure orificeslocated on 
the face and on the surface parallel to the free stream. The 
models are shown in Fig. 1. Several dummy models were em- 
ployed in the shock shape study for the determination of the 
detachment distance and sonic point on the shock. Schlieren 
photographs and static pressure distributions were obtained 
at Mach numbers between 11 and 19 at stagnation pressures 


between 700 and 1500 psia. ire oo 


Shock Shape Studies 


General 


Schlieren photographs were obtained for thick two-dimen- 
sional and axial symmetric bodies. For the flat plates, the 
leading edge thickness varied from 0.043 to 0.200 in. and for 
the axially symmetric bodies, from 0.063 to 0.5 in. A set of 
Schlieren photographs for the bodies is shown in Fig. 2. No 
boundary layer is noticed on either the two-dimensional or 
axial symmetric bodies. On thin flat plates (Re, < 4000) the 
edge of the boundary layer was well defined (21 and 22). 
From the photographs, the position of the shock with ref- 
erence to the body center line was measured on an optical 
comparator. The distance of the shock from the body on the 
body center line and the sonic point on the shock were also 
measured on the comparator. In the study of the shape of 
the shock about a flat plate with a square leading edge (21, 22) 
it was shown that the shape depended on both the geometric 
inviscid parameter x/t and on the viscous parameter (leading 
edge Reynolds number) Re, for small values of Re,, but only 
upon the geometric inviscid parameter for Re, > 4000. The 
leading edge Reynolds numbers for the present tests are about 
ten to fifty times larger than the above value. 

When the shock coordinates are nondimensionalized by the 
leading edge thickness and plotted on log-log paper, a straight 
line results. For the same body shape, the shock is independ- 
ent of the leading edge thickness and Mach number. As this 
is true, only a representative shock shape is shown in Fig. 3 for 
the four body configurations at a Mach number about 13. 
As these shocks form approximately straight lines on this 


plot, an equation of the form 
t 


The values of a and 6 remain ap- 


may be fitted to them. 
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Fig. 3 Shock shape of the four configurations i 


proximately constant for the entire range of Mach number and _ 
leading edge thickness. There is no effect of the leading edge 
Reynolds number on the shape of the shock. The shock shape 
obtained from the Schlieren photographs for the four configura- 
tions may be written as 


0.49 
= 2.0 () (for a flat plate with a square leading edge) 


0-50 
= 1.75 ; (for a flat plate with a round leading edge) 


( ) (for a cylinder with square leading edges) 


0.42 
1.4 
y \o-45 


(for a hemisphere cylinder) 


As y/t depended only upon 2/t for all values of Re,, the a 


viscous parameter did not enter into the relation, so one 
might assume that the phenomenon was primarily inviscid. 

The detachment distance 6 was measured from the Schlieren 
photographs for all bodies except the smaller ones. As this 
distance decreases, the inaccuracy of the measurements in- 
creases. In Fig. 4, the detachment distance nondimensional- 
ized by the leading edge thickness is shown as a function of 
the free stream Mach number. There is a certain amount of 
scatter in the points, but the average values are as follows for 
I2<M,< 19 


6/t = 0.48 (for the flat leading edge plate) 

6/t = 0.24 (for the round leading edge plate) __ 
5/d = 0.28 (for the flat leading edge cylinder) 
56/d = 0.08 (for the hemisphere cylinder) 


The sonic point on the shock was also obtained from the 
shock shape. In Fig. 5 the shock is sketched about the body 
for the four configurations. The sonic point and distance be- 


tween the shock and body are illustrated. 7 


Pressure Distribution Studies is 


Pressure distribution studies were made on a representative 
model for each of the configurations. The shock shapes found 
previously indicated that, for all the leading edges measured, 
the shape of the shock depended entirely upon the geometry. 
The pressures were obtained on one of the larger niodels used 
previously in the shock shape studies. Detail pressure dis- 
tributions were made with each of the configurations at Mach 
numbers between 11 and 14, and at a stagnation pressure of 
1000 psia. Pressures on the hemisphere cylinder were ob- 
tained to Mach 19 at a stagnation pressure of 1500 psia. 
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‘Fig. 4 (b) Detachment distance for the cylinder vs. 
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Due to the gradient in the tunnel, the static and total pres- 
sures are not constant along the tunnel center line. In an 
attempt to minimize the discrepancy that might be caused 
by this change in Mach number compared to a constant 
Mach number test, the static pressure p; and total pressure 
pr are assumed to be variables in the pressure coefficient ratio 
C,/Cpmax- This ratio is plotted against the distance parameter 
s/t where s is the surface distance measured from the midpoint 
of the body leading edge. The pressure coefficient, C,, is the 
difference between the measured pressure, p, and the static 
pressure p; (which would have existed at that point if no model 
were present), divided by the dynamic pressure; Cymax is the 
difference between the free stream total head pressure pr, and 
pi, divided by the dynamic pressure. For each configuration, 
two plots are made of ¢p/¢pmax, one near the leading edge where 
the pressure coefficient ratio is close to unity, and one on the 
afterbody, where the pressure coefficient approaches zero. 
On the latter plot, the approximate value of zero pressure in 
the tunnel is indicated. By this means, the details of the 
measured pressures near the leading edge and after the 
shoulder can be observed clearly. In Fig. 6 the pressure dis- 
tributions are shown for the four configurations at the 
various Mach numbers. In Fig. 6 (d) separation on the rear 
part of the plate at , = 13.9 causes the pressure coefficient to 
increase. 


Discussion 


In some previous work (21 and 22), interferograms were ob- 
tained of the flow about a plate with a flat leading edge varying 


Shock Shapes—Flat Plates 
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in thickness from 0.17 X 10-* to 59.5 X 107% in. From the 
interferograms of the thicker plates, a low pressure region was 
indicated immediately behind the corner From the Schlieren 
photograph of the two-dimensional bodies, Figs. 2(a), 2(b), 
it is possible to locate the shock clearly. The dark region at 
the corner of the plate is due to optical distortion through the 
high gradients that exist at that point. No details of the 
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ig. 5 (a) Shock about the flat leading edge plate 
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Fig. 5 (b) Shock about the hemicylindrical leading edge plate 
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(c) Shock about the flat leading edge cylinder 


Fig.5 (d) Shock about the hemisphere cylinder 
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boundary layer on the plate are visible because of the low 
density region that exists near the surface of the plate. 

Several analytical investigations have been made in an at- 
tempt to determine the shock shape, detachment distance 
and sonic point on the shock for a body with a round leading 
edge.® In the analysis of Hayes (10), the shock layer is as- 
sumed to be thin, so the radius of curvature of the shock is 
approximately that of the body. Between the shock and the 
body, incompressible, constant density flow is assumed. Ex- 
pressions have been obtained for the detachment distance for 
a plate with a hemicylindrical leading edge. Another solutio1 
has been obtained for the detachment distance (13) whicl 
is similar to that of Hayes. The shock is assumed to b 
parabolic and compressibility has been included in the calcu- 
lations. Chester has determined the flow characteristic: 
assuming a parabolic shock and calculated a body shape whic! 
is “nearly parabolic” (11). An exact solution is reached fo: 
7 — land M — o followed by a perturbation in (y — 1)/(* 
+ 1)and M-?, 

From the present series of tests of a plate with a flat leading 
edge, the shock can be approximated by a parabola of the 


form 
t 

For helium, (y = 5/3) at a Mach number of about 15, the angle 
that the shock makes with the free stream for sonic flow be- 
hind the shock is 63} deg. The sonic point on the shock is at 
y/t = 0.97 and lies above the body sonic point. The detach- 
ment distance is about 0.48. A sketch of body and shock show- 
ing the sonic point and detachment distance is shown in Fig. 
5 (a). 

A parabolic shock is also formed around a blunt two- 
dimensional body with a rounded leading edge The shock 


shape is of the for 
shape Is 0 orm 


(") = 3.06 (7) 
t 


The radius of curvature of the shock is 3.06 times that of the 
body or R,/t = 1.53 and the detachment distance 5/t = 0.24. 
The sonic point is at y/t = 0.66. 

The detachment distance as predicted in (10) for a hemi- 
cylindrical plate is 


kR 


Vi1 2k 


where & is the density ratio across the normal shock. In 
this reference, it is suggested that the radius of curvature of 
the shock R, be used instead of the radius of curvature of the 
body R. 


R,=R+6 


By this method, 6/t = 0.23. If, however, R, is obtained from 
the shock itself defined by the above equation, then the de- 
tachment distance is about 0.487. In (11) the radius of curva- 
ture of the shock on the axis of symmetry was calculated to be 
twice that of the body curvature. This value is less than the 
present one. The detachment distance predicted by this 
method is 0.21¢ which is reasonably close to the value of about 
0.24¢ obtained experimentally. From (11) the sonic point on 
the shock is at y/R, = 0.5 or y/t = 0.5 as the radius of curva- 
ture of the shock is twice that of the body. The sonic point 
on the shock is calculated to be outboard of the sonic point on 
the body. This value is again less than the experimental 
value, but has the same tendency. In Fig. 5 (b), the shock 
form, detachment distance and sonic point are noted. 


6 For blunt two-dimensional bodies a solution for the detach- 
ment distance and the y coordinate of the shock sonic point has 
been obtained (9). It is assumed that the shock is hyperbolic 
and asymptotic to the free stream Mach lines, and that the shock 
location is determined by continuity. 
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Pressure coefficient on the face of the flat leading edge 
plate 


Fig. 6 (a) 


Shock Shapes—Axial Symmetric 
The shock for the flat leading edge cylinder can be ap- 
proximated by the curve 


2.4 

(:) = 2.23 

) 

This equation holds for the measurements made to determine 
the shock shape for x/d > 0. If the radius of curvature were 
calculated using this equation, the value at the z-axis would 
be infinity. The sonic point is located at y/d = 0.67 and the 
detachment distance is at 6/d = 0.28. A sketch of the shock 
with the detachment distance and sonic point is shown in Fig. 
5 (ce). The sonic point on the shock is outboard of the body 
sonic point.? 

A considerable amount of experimental data is available 
on hemisphere cylinders or spheres at lower Mach numbers 
(3, 16 and 17). Reasonable correlations have been made with 
air (y = 1.4) and argon (y = 1.67) with some of the theories. 
For the present sets of tests, the shock may be defined as 


2.22 
0 


Again, the radius of curvature of the shock on the axis of 
symmetry is infinity. In an attempt to determine the radius 
of curvature, an equation of the form 


y \* y\° 
- = As 413 


was assumed and the coefficients evaluated. Using two terms 


y\? 
- = 0.709 { - 0.106 { 
(i) 


and with three terms 


x _fy\? _. fy\* 
— = 0.69 = 0.175 { - — 0.0 
q 397 @) + 0.175 (1) 0.074 (1) 


Omitting the higher order terms, this equation can be re- 7 


7 For a flat faced cylinder, Serbin (14) obtained a value for the 
detachment distance of 0.29d which is nearly equal to the ex- 
perimentally determined value. The detachment distance pre- 
dicted by Hayes is approximately twice the measured value. 
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Fig. 6 (b) Pressure coefficient on the surface of the flat leading 
edge plate parallel to the free stream direction 
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Fig. 6 (c) Pressure coefficient on the hemicylindrical section of 
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Fig. 6 (d) Pressure coefficient on the surface of the hemi- 
cylindrical leading edge plate parallel to the free stream direction 
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The radius of curvature of the shock on axis is about 1.4 times 
that of the body. 


NEWTONIAN 


The detachment distance was equal to 0.08d and the sonic 


point located at y/d = 0.37. It is noticed that the radius of 
curvature of the shock is approximately equal to the radius 
of curvature of the body plus the detachment distance. 


SONIC POINT 


Several of the theories give values for the detachment dis- 
tance with reasonable accuracy as are shown in Fig. 4. In 


(12), the detachment distance was calculated to be 0.067 


times the radius of curvature of the body and the radius of 
curvature of the shock about 1.25 times that of the body. 


It should be pointed out that many of the theories predict the 
value of the detachment distance for y ~ 1 andk—> 0. In 


helium the value of y is large compared to unity and k is large 
compared to zero. The results obtained from the theories 


‘an only give an indication as to the magnitude of this di-- 
tance. 


It has been suggested (9) that the characteristic features of 
the flow in the subsonic region behind the strong shock a1 


Fig. 6 (e) Pressure coefficient on the face of the flat leading dependent upon the choking phenomenon that occurs at the 


edge cylinder 


minimum area between the body sonic point and shock. A 
relationship has also been shown to exist between the sonic 
point on the shock and body sonic point, and that any changes 
in the body ahead of the body sonic point would have negligi- 


ble effect on the shock sonic point. For the two-dimensiona! 
bodies, the ratio of the y-coordinate of the shock sonic point 
to the y-coordinate of the body sonic point is about constant, 
but this is not true for the case of the axially symmetric 
bodies. The value of the detachment distance calculated 
by this method is about a third higher than the measured 
values. 


Pressure Distribution Study 


It has been shown experimentally that the pressure on a 
plate with a flat leading edge depends mainly on the leading 


q | | | 


edge Reynolds number Re, for thin leading edges or on the 
geometric parameter s/t for thick leading edges with the 


‘Fig. 6 (f) Pressure coefficient on the surface of the flat leading inviscid effects predominating for Re, > 4000 (21 to 23). De- 


edge cylinder parallel to the free stream direction 


tailed pressures were obtained for s/t > 5 for a plate with a 
flat leading edge. In the present tests, pressure distributions 
are presented for s/t < 10 for the inviscid model. 

On the face of the plate with the flat leading edge the pres- 
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sure has a maximum value in the center of the plate and de- 
creases as the corner is approached. The pressure at the sonic 
o\ dysin point is about one half the total pressure behind the norma! 
shock and is noted on the pressure curves. The pressure 
me bes measured at the nearest point behind the corner drops to 
. . I . 
a within 10 per cent of the static pressure in the tunnel and then 
ai ne increases to about Cp/Cpmax = 0.06, at s/t = 23. The pressure 
continuously decreases for s/t > 23. The pressure coefficient 
2 \ _sonic POINT 
Pmox | 
aan 
\ 
Fig. 6 (g) Pressure coefficient on the hemisphere section of the Fig. 6 (h) Pressure coefficient on the surface of the hemisphere 
cylinder leading edge cylinder parallel to the leading edge 
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increases with Mach number on the plate surface parallel to 


the free stream though there is no consistent variation of «aa 
It 


on the plate face which is normal to the flow direction. 


has been shown that for y = 5/3, a 90 deg expansion teens 5 


the Mach number to infinity with the pressure going to zero 
(22). The expansion waves around the corner strike the 
sonic line and return as compressions resulting in the pressure 
rise after the corner. The remaining part of the expansion fan 
hits the shock and returns as expansion waves causing the de- 
crease in the pressure on the plate. 

For the case of a plate with a hemicylindrical leading edge, 
the pressures on the cylindrical section are reasonably well 
predicted by the Newtonian approximation corrected for 
shock losses. However, near the junction of the cylinder sec- 
tion with the flat section, measured pressures are considerably 
higher than the predicted ones. For 6, = 0, (the plate 


parallel to the initial flow direction), the predicted value of — 


the pressure coefficient is zero. 

From the previous work it was shown that the viscous 
effects are unimportant for Re, > 4000. The pressures on the 
plate with the flat leading edge were obtained from two models 
with leading edge Reynolds numbers 37 X 10% (t = 0.0438 in.) 
and 105 X 10% (¢ = 0.124 in.). For the hemicylindrical] lead- 
ing edges, the leading edge Reynolds numbers were 37 X 10° 
(¢ = 0.0438 in.) and 170 X 10° (¢ = 0.20 in.). Good agree- 
ment is obtained between the same type model for the varia- 
tions of Re, These pressures, made at a Mach number of 
12.7, are presented as the incremental pressure ratio A p/p; 
(Fig. 7). On this plot, the pressure decays as the —0.5 power 
of s/t. The measured pressures for the two different leading 
edges are approximately equal for s/t > 3. 

Several experimental investigations have been made on 
axially symmetric bodies with round and square Jeading edges 
at various Mach numbers (2, 3, and 4). For the square lead- 
ing edges, Oliver has shown that at M = 5.8 in air, the pres- 
sure on the face is relatively constant and starts to drop off 
at about 0.35 s/d. No pressure measurements were made after 
the corner. The pressure after the shoulder on a flat leading 
edged cylinder behaves in a manner similar to the flat lead- 
ing edge plate. The expansion followed by compression waves 
returning from the sonic line and expansion waves from 
the shock wave cause the low pressure immediately after the 
corner, the increase in pressure and final pressure decay. 
The pressure coefficient has a maximum value of about 0.3 at 
s/d = 2. 

For the hemisphere cylinder, it has been shown that the 
pressure on the curved section is well predicted by the New- 
tonian pressures for Mach numbers between 2.0 and 6.8 in air. 
The pressures obtained in helium for Mach numbers between 
11 and 19 on the hemisphere section can also be predicted 
reasonably well by the Newtonian approximation, but the 
measured pressures diverge from the predicted ones near the 
junction of the hemisphere with the cylinder. The pressures 
on the cylindrical section approached the ambient pressure 
but were always greater than it. On the cylindrical section, the 
pressure coefficient increases with Mach number though there 
is no consistent variation of pressure coefficient with Mach 
number on the spherical portion of the body. The pressures 
on the cylindrical section of the hemisphere cylinder are 
lower than those on the flat leading edged cylinder for s/d > 1. 

On the surface parallel to the free stream, the pressures on 
the axial symmetric bodies are less than those on the two- 
dimensional bodies for both the flat and rounded leading 
edges. For the flat leading edge, the peak value of the pres- 
sure coefficient for the plate is twice that measured on the 
cylinder. For the round leading edges, the pressures on the 
plate after the junction were about two to three times those 
on the cylinder for the same vales of s/f. ~~ 6 


1 A parabolic shock is formed about plates with flat and 


Conclusions 
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Fig. 7 Incremental pressure ratio vs. dimensionless distance 
for the four configurations at 1f = 12.7 


hemicylindrical leading edges for z/t < 10. For cylindrical 
bodies, a “nearly” parabolic shock is formed for x/d < 10. 

2 The value of the detachment distance for the two- 
dimensional and axial symmetric bodies with a rounded lead- 
ing edge is within the order of magnitude of the predicted 
values. 

3 The sonic point on the shock depends both upon the 
body sonic point and body shape ahead of the sonic point. 

4 On the body surface parallel to the free stream direction, 
an increase in Mach number causes an increase in the pressure 
coefficient ratio for both two-dimensional and axial sym- 
metric bodies. 

5 For the same type of leading edge, the pressure co- 
efficient ratio on the surface parallel to the flow direction is 
higher for a two-dimensional body than for an axial sym- 
metric body. 

6 On the flat and hemicylindrical plates at the same free 
stream Mach number, the pressure coefficient ratios are equal 
for s/t > 23. 
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In order to obtain a realistic value of heat flux and tem- 
perature distribution within the walls of a rocket thrust 
_ chamber recourse must be taken to experimental and em- 
y pirical techniques. This paper discusses the experimen- 
- tally determined heat transfer data obtained on JP-4 and 

RFNA under conditions simulating those occurring during 
thrust chamber operation. An electrical analogy tech- 
‘nique is used to determine the temperature and heat flux 
distribution because of the irregular shape of the coolant 
passages. Using the experimental data for film coefficients 
and a two-dimensional electrical analog, the heat flux and 
temperature distribution within the walls of a tubular 
chamber cooled with JP-4 and within the walls of a drilled 
_ passage thrust chamber cooled with RFNA are obtained. 


q 


Nomenclature 


G = mass velocity, lb/in.?-sec 
hy = heat transfer coefficient 
ee = metal thermal conductivity 
Ly = length of paper representing h, ei 
Ly = length of contoured paper representing hy -— 
Li = length of paper representing gas film 
L's, = length of contoured paper representing liquid film 
Nu = Nusselt number 
p = perimeter of section for which the heat flux is 
calculated 
r = absolute pressure, psia 
ry = Prandtl number 
Q/A = heat flux, Btu/in.*-sec 
(Q/A)cony = predicted heat flux for convection at the same wall 
and bulk temperature that exists at (Q/A)ui 
(Q/A)u: = heat flux at the upper limit of nucleate boiling 
r = radius of curvature La 
Re = Reynolds number 
Ry = paper ratio = 
resistance/[_Jof paper representing h 

resistance/[_Jof paper representing metal 
S = scale factor size of model/size of prototype 
AT = temperature difference 


— in Thrust Chambers 


Aerojet-General Corp., Sacramento, Calif. 


_ characteristics under these conditions. 


Introduction 


HE problems of high coolant velocity, boiling, deposition 

of carbon and other impurities, high temperature and 
velocity of combustion gases, and odd geometry make a pre- 
cise analytical treatment of heat transfer in regeneratively 
cooled rocket engines very difficult or impossible. Therefore, 
in order to obtain realistic values of heat flux and temperature 
distribution within the walls of a rocket thrust chamber re- 
course must be taken to experimental and empirical tech- 
niques. 

Rocket thrust chambers designed for long duration firings 
must be cooled to protect the walls from the high tempera- 
ture, high velocity combustion gases. This is usually ac- 
complished by passing one of the propellants through coolant 
passages located within the walls of the thrust chamber. At 
the high heat fluxes and temperatures encountered in rocket 
thrust chambers the propellants may become corrosive, de- 
compose, or deposit impurities on the heating surface, thereby 
reducing their cooling effectiveness. It is impossible to pre- 
dict, with any degree of certainty, the propellant heat transfer 
To properly analyze 
the heat transfer conditions in this region of high heat flux 
and temperature, it is necessary to base the analysis upon ex- 
perimentally determined heat transfer data. 

The use of one-dimensional heat conduction equations will 
give erroneous results when applied to the irregular geometric 
configurations which are sometimes encountered in thrust 
chamber designs. The method proposed to solve the two- 
dimensional] heat transfer conduction problems involved (1)? 
is an electrical analogy technique. 

By making use of the experimentally determined liquid film 
heat transfer coefficients and the electrical analogy method 


Presented at the ARS Semi-Annual Meeting, San Francisco, 
Calif., June 10-13, 1957. 

1 Development Engineer, Liquid Rocket Plant. 

2 Numbers in parentheses callings References at end of paper. 
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Fig. 1 


for solving two-dimensional heat transfer problems, the heat 
flux and temperature distribution within the walls of the 
combustion chamber may be quite accurately determined. 
In this paper a heat transfer analysis is made to determine 
the heat flux and temperature distribution within the walls of 
tubular and drilled passage thrust chambers which are cooled 
with JP-4 (2) and RFNA (3), respectively. 


Solution of Two-Dimensional Heat Transfer 


Problems 


General Considerations 


The use of simplified conduction equations will give er- 
roneous results when applied to the solution of the heat trans- 
fer conditions which occur when irregular coolant passage 
geometries are involved. It has been shown in (1) that two- 
dimensional steady state heat transfer problems involving 
irregular geometries may be easily solved by the use of elec- 
trically conductive paper. The basic assumptions used in 
analyzing heat conduction problems by this method are: 

1 The similarities between electrical flow and heat flow 
are applicable to the problems under consideration. 

2 Two-dimensional heat flow predominates. 

3 The metal thermal conductivity is independent of 
temperature. 

4 The various film coefficients are correctly represented 
even though geometric similarity is not achieved. 


Electrical Analog 


The test installation used at Aerojet to solve these heat 
transfer problems is shown in Fig. 1. The base of the plotter 
is made from 3 in. Masonite; fastened to this base are two 
steel bars which serve as guides for the pantograph tracer. 
The pantograph tracer consists of an aluminum U beam; 
located inside the U beam is another slider which holds the 
voltage probe and the marker. The marker will contact 
the paper only when force is applied. A weight is fastened to 
the brass voltage probe so that a uniform pressure may be 
applied to the conductive paper surface at all times. 

A schematic of the power source and voltage measuring 
device is shown in Fig. 2. Power is supplied by a 6-volt lantern 
battery. Four 10-turn 200-ohm, 2-watt potentiometers are 
connected in parallel with the battery. By reading the dials 
on the potentiometers, the per cent of battery voltage may be 
determined. Leads from these potentiometers are connected 
to the model. The per cent of battery voltage is proportional 
to the temperatures imposed as boundary conditions on the 
model. A 1000-ohm 10-turn potentiometer is used for ac- 
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Electrical analogy test equipment 


curate measurement of the per cent of battery voltage. 
Zero deflections on the microammeter indicate when the value 
read on the dial of the 1000-ohm helipot is the per cent volt- 
age the probe is measuring. Since there is no current flow 
in the measuring circuit, the resistance of the lines will not 
affect the electrical measurement. Although not shown in 
the schematic, a potentiometer circuit has been installed 
in the measuring circuit which enables the resistance of the 
teledeltos paper to be determined. 

The electrically conductive paper was glued to a piece of 
} in. Masonite. This backing offers the necessary rigidity re- 
quired by the mode]. Rubber cement was used as the bonding 
material because it provided a good bond and did not react 
with the paper. It was clean to use and is easily removed from 
the Masonite after the test. 

Silver dots, made from electrically conductive silver paint, 
were used as the connecting medium at the boundary surfaces. 
A series of dots rather than bands were applied as the bound- 
ary connection in order to eliminate heat transfer along the 
boundary surface even though the exact size and spacing of 
these dots did not materially affect the isothermal lines. Dots 
of approximately 7g in. diam and } in. between the centers 
were used throughout these tests. 

To impress constant temperature (voltage) boundary con- 
ditions on the model, electrical terminals were provided by 
painting bands of silver paint along the boundaries. The gas 
temperature was impressed on a 3 in. silver band on the 
model’s gas film and the liquid bulk temperature was im- 
pressed on a silvered area located within the coolant passages. 

Voltage to the various terminals was adjusted by means of 
a variable potential source. The probe was placed in contact 
with each terminal in turn and the voltage adjusted by a 
variable resistor to a predetermined value. The form on which 
the plotting is to be done is fixed under pantograph tracer. 
The usual procedure in plotting was first to measure tempera- 
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Fig. 2 Electrical analogy voltage measuring circuit 7 
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ture (voltage) along the combustion wall by placing the probe 
on the silver connecting bands or dots. The potentiometer 
was set for the value desired and the probe was then used to 
bracket carefully that point on the model test section by ob- 
serving the point at which the microammeter registered zero 
current flow. The isothermal line was then plotted across 
the model before continuing with another isotherm. 

To obtain the required length of teledeltos paper which will 
represent the film heat transfer coefficients, the following 
formula was used 


I K,,8 
where 


resistance/(_]* of paper representing h; 


resistance/(_Jof paper representing metal 


In the case of applying a film to a curved surface, the L, 


calculated from Equation [1] should be applied to the surface 
as shown in Fig. 3 in order that one-dimensional heat transfer 
_ in the boundary layer may be simulated. 


It is possible to calculate a length of paper having the same 


“3 The term ohms/square is the resistance of a unit square of the 


_ electrically conductive paper which remains constant regardless 
of the size of the square. 
: silver strips located on opposite sides of the square. 


The resistance is measured between 
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resistance as the film but fitting the contour of the curved sur- 
face by the use of Equation [2a or 2b]. (Derivation of these 
equations is shown in the Appendix.) 


convex surface: 


convex surface: 
=Inr —In(r — L’).............. [2h] 


where L, is the length of paper calculated by Equation [1] 
and L’ is the length of contoured paper. 

In representing film coefficients by a curved piece of paper, 
slits must be cut in the paper to ensure that one-dimension:| 
flow occurs in the simulated film resistance. 

Once the isothermal lines are determined, the heat flow lines 
are then determined graphically. In the graphical analysis, 
the heat flow lines are drawn perpendicular to the isothermal 
lines. The space between two heat flow lines constitutes a 
heat flow tube. The amount of heat flowing in a heat flow 
tube can then be calculated and consequently the localize: 
value of heat flux can then be determined. 


~~ Liquid Film Heat Transfer Coefficients 


as General Discussion 


Rocket thrust chambers designed for long duration firings 
must be cool to protect the walls from the high temperature, 
high velocity combustion gases. This is usually accomplished 
by passing one of the propellants through passages located 
within the walls of the thrust chamber. At the high heat 
fluxes and temperatures encountered in thrust chamber opera- 
tions the propellants may become corrosive, decompose, or 
deposit impurities upon the heating surface thereby reducing 
the cooling effectiveness. It is impossible to calculate, with 
any degree of certainty, the propellant heat transfer charac- 
teristics under these conditions. Experimental test tech- 
niques have been developed which enable the heat transfer 
characteristics of various coolants to be evaluated at condi- 
tions simulating those encountered during thrust chamber 
operations. 

Shown in Fig. 4 is the heat transfer process which would be 
predicted for a pure nonreacting liquid like liquid nitrogen 
(4). 

In the nonboiling region A-B, the thermal resistance to heat 
transfer remains essentially constant. The heat transfer 


_ process will remain in the nonboiling region providing the 


wall surface temperature remains below its saturation or 
critical temperature. At increased rates of heat transfer, the 
mode of heat transfer will be influenced by the pressure im- 
posed upon the coolant. 

At subcritical pressures, nucleate boiling (B-C) will occur 


at the heated surface when the wall surface temperature 


slightly exceeds the saturation temperature of the coolant. 
The increased turbulence caused by the fast-moving bubbles 
greatly reduces the thermal resistance to heat transfer. In 
this region it is possible to increase the rate of heat transfer 
with only a slight increase in the heating surface temperature. 
As the rate of heat transfer is further increased (B) the bubbles 
start to coalesce and form a high resistance vapor film. If the 


_ thermal resistance offered by this vapor is high enough, the 


temperature of the wall may exceed its melting point, or the 
working pressure may exceed the stress limit of the wall at 
these high temperatures. 

At supercritical pressure, the nonboiling condition exists 
at heating surface temperatures below the critical (A-B). 
Beyond this point, a gradual transition to a stable film boiling 
boundary layer begins (B’E). A stable film boiling region 
exists from E to F. 
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The upper limit of heat flux occurs when the heating sur- 
face temperature reaches its melting point or the stress limit 
of the heating surface is exceeded at these elevated tempera- 
tures. 

The heat transfer conditions described are what might be 
expected when a pure, nonreactive and stable liquid is used 
as a coolant. Unfortunately, a rocket propellant may not 
always meet these requirements. At low rates of heat trans- 
fer, the heat transfer coefficients may be estimated quite 
adequately by either theoretical equations or empirical equa- 
tions based upon heat transfer measurements of other 
liquids. 

Since it is impossible to predict the heat transfer charac- 
teristics of rocket propellants in the nucleate and film boiling 
region, rocket designers are forced to base their heat trans- 
fer calculations upon experimentally determined heat transfer 
coefficients. 

Rocket Propellants 


Jet Fuel (JP-4) 

The heat transfer characteristics of jet fuels have been ob- 
tained on JP-3 and JP-4. If JP-4 was used as a regenerative 
coolant at a pressure below its critical pressure, the heat 
transfer process would have been in the nonboiling region pro- 
vided the wall temperature was below the saturation tem- 
perature of JP-4. The nonboiling heat transfer coefficient for 
JP-4 may be estimated by 


If the heat transfer rates were sufficiently high the heat 
would have been transferred in the nucleate boiling region. 
The heating surface wall temperature remained essentially 
constant in this region. The upper limit of nucleate boiling 
‘ould be estimated from 


(Q/A) ut 10° 
(Q/A )conv PG 

In regenerative cooling of rocket thrust chamber with JP-4 
the cooling is usually carried out under pressures which exceed 
the critical pressure of the coolant. In this case the heat 
transfer process would proceed as follows: A nonboiling con- 
dition exists up to a heat flux at which the wall in contact 
with the fuel reaches the critical temperature of the fuel. The 
film boiling data obtained on this hydrocarbon mixture may 
be compared with the film boiling data obtained for liquid 
nitrogen as follows: 

1 With liquid nitrogen, there was a gradual transition 
into film boiling, whereas in the case of JP-4, there was an 
abrupt transition. 

2 With liquid nitrogen there was a constant film cooling 
heat transfer coefficient which was less than the nonboil- 
ing heat transfer coefficient. With JP-4, the heat transfer 
coefficient did not remain constant in the film-boiling region. 

3 The upper limit of heat flux was obtained when the 
temperature of the heating tube reached a point where the 
metal was not able to withstand the external pressure of the 
system. Whether the film boiling heat transfer obtained on 
JP-4 would have followed the same trend as the Jiquid nitro- 
gen at increased pressures is not known. 

The trend of JP-4 heat transfer data at supercritical pres- 
sures is shown in Fig. 5. The solid line is the experimentally 
determined heat transfer data obtained from (2). The dashed 
line is the extrapolated heat transfer data used in the analysis 
of the heat transfer occurring in the tubular thrust chamber. 
It is felt that the errors introduced extrapolating the non- 
boiling heat transfer data in the range shown will be small. 
The errors introduced by extrapolating the 500-psia data to 
higher pressures are not known and further experimental work 
is warranted. 
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Fig. 6 Heat transfer characteristics of RFNA at 700 psia 


Nitric Acid 
Heat transfer tests conducted on nitric acid showed that 
the nonboiling heat transfer coefficient could be estimated by 
[3]. 
Nu = 0.027Re8Pr-4,. 


In the case of WFNA (6), nucleate boiling occurred at a 
surface temperature 50 to 100 F above the saturation tem- 
perature of the acid. As the heat flux increased, the tem- 
perature increased slightly. A light scale was formed on the 
heating surface. At the upper limit of nucleate boiling, burn- 
out occurred. In the case of red fuming nitric acid, containing 
only traces of impurities, and no HF inhibitor, a similar trend 
of data was observed. RFNA*‘ containing only traces of im- 
purities and 0.4 per cent HF deposited a large quantity of salt 
scale’ on the heating surface, thus decreasing its effective 
heat transfer coefficient.6 The trend of heat transfer data 


4RFNA HNO;—82%, NO.—14%, H,O—2.5%, HF—0.6%. 

5 The scale is composed of mainly iron, nickel and chromium 
salts. 

6 Refers to the calculated heat transfer coefficient based upon 
the temperature drop from the fluid bulk temperature to the metal 
surface temperature. 
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that can be expected for RFNA at 700 psia is shown in Fig. 6. 
: The heat transfer condition may be estimated in the non- 
S ing boiling region by Equation [5] up to a AT’ of 225 F. Scale 
_ then forms on the heated surface and the effective heat trans- 
fer coefficient decreases. With the absence of experimental 
_ heat transfer data the rocket engine designer would be apt to 
- conduct heat transfer calculation assuming nucleate boiling 
on a clean surface, which would lead to erroneous results. 
The effect of this scale in reducing the effective heat transfer 
coefficient at high rates of heat transfer can be seen in the 
- following data taken from (3). At a heat flux of 7 Btu/in.?- 
_ sec, the temperature drop across the liquid film was found to 
be only 239 F whereas the temperature drop across the scale 
was estimated to be 816 F. If the heat transfer coefficient for 
a clean surface had been 0.0292 Btu/in.*-sec, the presence of 
scale reduced the effective heat transfer coefficient 77 per cent 
to 0.0066 Btu/in.*-sec. 


Application to Thrust Chambers 

By making use of the experimentally determined liquid film 
heat transfer coefficients and the electrical analogy method 
for solving two-dimensional heat transfer problems, the heat 
and temperature distribution within the walls of the combus- 
tion chamber may be determined. 

The heat transfer condition existing in a 0.030-in. wall 
tubular aluminum thrust chamber being cooled with JP-4 
was investigated. The nonboiling liquid film heat transfer 
coefficient which was calculated according to Equation [3] 
was found to be 0.0082 Btu/in.*-sec-°F. Using the same 
weight flow of fuel, the nonboiling heat transfer coefficient 
existing in the finned tube was found to be 0.017 Btu/in.?-sec- 
°F. Details of the boundary conditions used in this investiga- 
tion are shown in Table 1 


Table 1 Boundary conditions for the tubular thrust 
chamber 
Unfinned Finned 
tube tube 
Gas film heat transfer coeffi- 
cient, Btu/in.?-sec-°F 0.00085 0.00085 
Liquid film heat transfer co- 
efficient, Btu/in.*-sec-°F 0.00823 0.0173 
Thermal conductivity of alu- 
minum (Btu/in.-sec-°F ) 0.00207 0.00207 
Coolant temperature, °F 150 150 
Combustion temperature, °F 5230 5230 
Scale factor 20/1 20/1 
Ry 1 1 
Ly,, inches 6.9 6.9 
Ly, inches 0.7 0.34 
L'y,inches 0.54 0.30 


FINNED ALUMINUM TUBE 


ALUMINUM TUBE 


Fig. 7 Temperature distribution within the walls of a tubular 
thrust chamber 


Shown in Fig. 7 are the isothermal profiles obtained from 
the electrical analogy test technique. These profiles represent 
a quarter section of a tube, which makes up the throat sec- 
tion of a thrust chamber. The hot combustion gases are as- 
sumed to have equal effectiveness over the total surface area 
exposed to the hot gases. 

In the case of the unfinned tube the gas side wall tempera- 
ture was around 700 F with an inside wall temperature of 650 
F. A maximum heat transfer rate of 3.7 Btu/in.?-sece was ob- 
tained on the liquid side. Because the isothermal lines were 
essentially straight and parallel, the heat flux of 3.7 Btu/in.?- 
sec was obtained on the gas side. In the case of the finned 
tube, where the JP-4 becomes a more effective coolant, the 
gas-side temperature is only 510 F and a liquid side wall 
temperature of 445 F was obtained. Accompanying this re- 
duced wall temperature was an increase in the rate of heat 
transfer. The maximum heat flux into the coolant passage 
was found to be 5.5 Btu/in.*-sec. As can be seen from Fig. 5 
this high rate of heat transfer can be easily absorbed by the 
coclant. 

It is noted that the nonboiling heat transfer condition was 
present in the coolant passage. This is indicated by the 
fact that the coolant side wall temperature was below the 
critical temperature of JP-4 and consequently below the heat 
flux where film boiling may be expected. 

The heat transfer conditions existing in a 0.165-in. drilled 
passage thrust chamber which is cooled with RFNA were in- 
vestigated. The results of this investigation are shown in 
Fig. 8. These sections represent a radial section cut from the 
thrust chamber. Hot gases are flowing past the bottom sur- 


ce eid. Table 2 Boundary conditions for the drilled passage thrust chamber 
Nonboiling Nonboiling Boiling aluminum 
stainless steel aluminum clean surface scaled surface 
Gas film heat transfer coefficient, Btu/in.?-sec-°F 0.0020 0.0020 0.0020 0.0020 
Liquid film heat transfer coefficient, Btu/in.*-sec-°F 0.011 0.011 = 0.011 and 0.011 and 
0.0060 
Thermal conductivity of metal, Btu/in.-sec-°F 0.00035 0.0025 0.0025 0.0025 
Coolant temperature, °F ee we 75 75 a 75 
Combustion temperature, °F al 4500 4500 4500 4500 
Scale factor 20/1 20/1 | 20/1 
33 jj|§§ 8.35 3.38 
Ls, inches ‘ a 0.529 0 506- 0.506- 1.545 to 
bine 
1.639 
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Fig. 8 Temperature distribution within the walls of a drilled passage thrust chamber ; 
face and RFNA is flowing through drilled passage. The Also shown in Fig. 8 are the heat transfer conditions en- 
boundary conditions used in this analysis are shown in Table 2. countered when heat is transferred in the nucleate boiling 
The nonboiling liquid film heat transfer coefficient was calcu- region. In the absence of experimental data, the most logical 
lated by Equation [5]. The heat transfer conditions in the means of estimating heat transfer conditions in the nucleate 
boiling region were estimated from the data shown in Fig. 6. boiling region is to assume that the coolant wall temperature 
A comparison between the heat transfer conditions existing remains at a constant temperature slightly in excess of the 
in an aluminum chamber and those existing in a stainless saturation temperature of the liquid, in this case, 380 F. This 
steel thrust chamber is shown in Fig. 8. The results of this was accomplished in this analysis by imposing a constant 
comparison are summarized in Table 3. temperature along the coolant where the heat flux is in excess 


of 3.3 Btu/in.?-sec. Based upon experimental data of (3), 
there is a considerable reduction in the heat transfer coefficient 
due to scale formation. As determined from Fig. 6, the heat 


Table 3 Heat transfer conditions in aluminum and transfer coefficient was equal to 0.0060 Btu/in.*-sec-°F at 
stainless steel thrust chambers heat fluxes above 3.3 Btu-in.*-sec. This condition was 
Stainless analyzed by applying coefficients on the liquid side varying 

Aluminum steel between 0.011 and 0.006 Btu/in.?-sec-°F. 


Gas side wall temperature, °F 850 1850 A comparison of the heat transfer conditions encountered 
Liquid side wall temperature, 650 650 with a scaled and unscaled surface is shown in Table 4. 


oF 
Q/A gas side, Btu/in.*-sec 
Q/A liquid side, Btu/in.*-see 


Table 4 Heat transfer conditions encountered with 


| 


a scaled and unscaled surface 

As would be expected, the use of a lower thermal conduc- Ls ite 3 es Clean Scaled 
tivity material, such as stainless steel, gives decreased rates of ‘| 25 se Bh! surface surface 
heat transfer across the combustion wall and increased wall Gas side wall temperature, °F 650 1050 

temperatures. The shape of the isothermal lines, curving Liquid side wall temperature, °F = 380 850 
Gas side heat flux, Btu/in.®-sec 7.35 6.95 


away from the coolant passage of the aluminum section, indi- Liquid side heat flux, Btu/in.? 

cates that the aluminum adjacent to the passage tends to aid ee eee eee 
in the cooling process, whereas with the stainless steel section 
the isothermal lines tend to curve around the coolant passage 
causing heat to flow into the coolant rather than around the 


passage. Nucleate boiling on a clean surface offered very little re- 

In actual practice, the heat transfer condition shown here is sistance to the flow of heat, resulting in extremely high heat 
not possible because as soon as the coolant side wall tempera- transfer rates. Considerable resistance to the flow of heat 
ture reaches the saturation temperature of the coolant, boiling was offered by the nucleate boiling ona scaled surface. The 
will occur. isothermal lines tended to bend away from the coolant pas- 
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Fig. 9 Heat flux distribution around the coolant passage 


sage, thereby resulting in low heat fluxes at the coolant 
passage and an extremely high combustion wall temperature. 

Shown in Fig. 9 is a plot of the heat flux distribution around 
the coolant passage for the various coolant characteristics 
considered. The extremely high heat fluxes occurring with 
boiling on a clean surface is above the maximum heat flux the 
coolant is capable of absorbing. This condition could cause 
failure during thrust chamber operation. 

The lowest heat flux obtained at the coolant passage oc- 
curred with boiling on a scaled surface. This heat transfer 
condition is also critical due to the high wall temperature, 
which exceeds 1000 F. The application of an aluminum oxide 
coating to the combustion wall would help alleviate these 
severe heat transfer conditions.. As can be seen from this 
analysis, an extremely high heat transfer coefficient may be 
detrimental to the cooling of a drilled passage thrust chamber. 
With irregular geometric configurations, a two-dimensional 
heat transfer analysis based upon experimentally deter- 
mined heat transfer coefficients is a necessity to insure a satis- 
factory solution. 


Appendix: Derivation of Equations [2a and 2b] 


b, the thickness of the conducting paper. Then if 


= resistance to electric flow _ 
resistivity of the paper 
length of paper calculated by Equation [1] 

area of paper to flow of current = 2mrr a/360 

= included angle of cylindrical segment of curved surface 


Consider a segment of a cylinder of radius r and height 7 = 


Then the mean area to flow of the contoured paper should be 
similar to the area of flow of the model calculated by Equation 
{1]. That is, the film obtained from Equation [1] has the 
same flow area as the model curved surface, but the contoured 
paper must be adjusted in length to L’ and area Am, to give 
the same flow resistance. 


flow area on curved surface: 7 
a 
As = 2ar7 —— 
360 


which is also the flow area of the film — —_ Equation 
(1). 
The mean flow area of the film from 


a 
2 Lr 
A, 
in —— 

If resistances are to be equal a 


where L is the length of film paper from Equation [1]. 


2Qr(r + 
360 
(r+ L’) 
L 


= In(r + L’) — Inr 


Having L from Equation [1] and radius of curved surface r, 
the length of contoured paper L’ giving the same resistance as 
L can be calculated. Then the single contoured piece can be 
cut out and placed on the model with more ease than the 
strips as shown in Fig. 3. It may be that L’ will be too large 
for practical application and strips will have to be used, but 
this method may be applied in many cases. 

If the contoured paper film is to be applied to the concave 
side cf the curved surface, by a similar proof it can be shown 
that 


L 
~ =inr—In(r —L’) 
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The characteristics of an inertial navigation system fora 
simplified guidance problem are described. Error equa- 
tions relating component performance to guidance accu- 
racy are derived. These equations exhibit the 84 min 


period oscillations which typify inertial systems. 


Introduction 


NERTIAL guidance is a relatively new method and means 
of navigation of very general application which is quite 
distinct in principle from older forms, such as celestial navi- 
gation, dead reckoning, piloting and radio systems. The 
significance of this new navigation technique is now generally 
recognized. 

Several recent papers (1, 2, 3, 4)? have discussed the general 
characteristics of inertial guidance. The present paper is an 
extension of such work to include performance analysis. 

A simplified guidance problem is considered, namely, navi- 
gation on a spherical nonrotating body. This has been done 
primarily so that significant features can be emphasized with- 
out the mathematical complexity associated with the more 
general problem. The general characteristics of this system, 
as well as the results of the analysis, are, however, only 
slightly modified when extended to a more realistic guidance 


Rectangular Coordinates on a Flat Earth 


As a first and greatly oversimplified example of inertial ,- 


guidance, consider the problem of locating a cart on the top of 
a level table such as that shown in Fig. 1. 

The cart is initially at rest at some known position. A 
force is applied starting it along some path on the table. Its 
position at any time can be determined by measuring and 
doubly integrating the applied acceleration. This is the basic 
concept behind an inertial guidance system. 

A mechanization of even this simple problem would bring 
out many of the significant features. First, and perhaps 
most important of these, is that acceleration is a vector quan- 
tity. Each accelerometer can sense only one component of 
the total vector acceleration. Therefore, each instrument 
must not only be accurate itself, but must be maintained in a 
precisely known orientation. A second fact is that no 
accelerometer can distinguish between acceleration and 
gravity. Thus if the output of one of the acclerometers is to 
give the x-position, then the accelerometer must be aligned 
along the z-axis at all times. Misorientation will cause it to 
sense the y-component of either acceleration or gravity. 

For this example the inertial system could be conveniently 
split into three subsystems: 

1 The accelerometers which sense x and y components of 
the acceleration of the cart. These instruments, in principle 
at least, can be identical to the many varieties of commercial 
accelerometers currently available. 
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Analysis of an Inertial Guidance System — 


Quantitative analysis 


shows, however, that the accuracy requirements imposed on 
these instruments require precision beyond that now generally 
available. An earlier paper indicates some cf the practical 
problems involved in acceleration measurement in inertial 
guidance (5). No attempt will be made in this paper to 
amplify that discussion. 

2 Astable platform which controls the orientation of the 
accelerometers so that each senses only the correct component 
of acceleration. The most straightforward way of mechaniz- 
ing such a stable platform is to use an assemblage of gyroscopes 
arranged so as to maintain each axis angularly fixed in inertial 
space (6). Again, these instruments are, or can be, identical 
in principle to standard gyroscopes, but must satisfy much 
more stringent accuracy requirements. 

3 A computer which doubly integrates the outputs of the 
two accelerometers to determine position. (In some cases, 
however, the accelerometers are designed to perform one or 
both steps of integration within the instrument.) The com- 
puter can be mechanized using either analog or digital tech- 
niques, or a combination thereof. 

The operation of an inertial autonavigator for more realistic 
guidance problems is remarkably similar to this example. The 
navigator can still be split into the same three subsystems, 
and the functions of each subsystem, conceptionally at least, 
remain almost identical. 


- Rectangular Coordinates in a Spherical Non- 


Rotating Body 


When the size of the “table top” becomes of the order of 
miles or hundreds of miles, the system outlined above would 
develop prohibitively large errors. 

As the distance from the center increases, the accelerometers 
which remain parallel to their initial orientation will sense a 
lateral component of gravity proportional to the sine of the 
angle traveled. The double integral of this input would pro- 
duce a large and increasing error. For reasonable navigation, 
some method of correcting for this lateral gravity component 
must be found. If the distance traveled from the center of the 
coordinate system is not too large, a simple method can be 


Fig. 1 Table configuration 
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Fig. 2. Rectangular coordinate system — 


used; the computer determines these corrections as functions 
of zand y. Referring to Fig. 2, where 


a = radius of the earth (assumed to be spherical) 

R = distance of the vehicle from the center of the earth 

h = altitude of the vehicle above the surface of the earth 

Z,y = rectangular coordinates on the “‘table’’ whose origin is 
at the point where the table is tangent to the earth 

go = gravitational acceleration at the surface of the earth 

9 = gravitational acceleration at the location of the ve- 
hicle 

Jz, 9, = x and y components of the gravitational acceleration 


Expressions for these lateral components of gravity are 


The accelerations sensed by the xz and y accelerometers be- 


eome in this case 


A simple case which illustrates the principles occurs when 
departures from the origin are sufficiently small to neglect the 
factor (h/a). Under these conditions the expressions for the 
accelerations become 


= 
9+ 


The gravity corrections are therefore simple linear functions 


of zand y. The system can be mechanized in essentially the 


same manner as that considered before if the value of the 
_ gravity correction is computed and fed back as a correction on 
the accelerometer. A transducer is required to convert the 
computed value to an equivalent force on the accelerometer. 
A block diagram for the system in this case is shown in Fig. 3. 

The acceleration of the platform with respect to inertial 
space in the z and y directions is sensed by accelerometers 


position as an output. The gravity correction term is com- 


puted as a function of the output and introduced as a feedback 


or bias on the accelerometer. 


The situation described above where motion is confined to 


7 die in the “table” or zy-plane is, of course, highly artificial. 
In the more general case, position would be defined by the 


_ three coordinates (x, y, z) and the corresponding expressions 
acceleration become 


STABLE 
| PLATFORM 
INERTIAL 
REACTION 
FORCE x DOUBLE |p 
x 
COMPUTED 
TRANSDUCER 
INERTIAL 
RE ACTION 
FORCE. Y DOUBLE 
Ay INTEGRATO y 
COMPUTED 
TRANSDUCER 


Fig. 3 System block diagram (rectangular coordinates) 


A, =#é+ Jz 
Ay + Gy 
A, 


For small values of z, y and z, linear approximations to the 
gravity correction terms can be made, and these expressions 
reduce to 


= | 


a 


A mechanization of these equations would be identical in form 
to that given in Fig. 3. 


Cylindrical Coordinates 


As the distance from center of the coordinate system shown 
in Fig. 2 becomes comparable with the radius of the earth, the 
size of the gravity correction becomes of the order of 1 g. 
Accurate mechanization of such a large acceleration correc- 
tion term places stringent requirements on the computer and 
the transducer used to convert the computed value to an 
equivalent acceleration. While these problems are tech- 
nologically solvable, it is generally more reasonable to rotate 
the platform so as to keep the accelerometers normal to the 
local gravity vector. For a general path this requires rotation 
about two axes. Often, however, the path of the vehicle may 
lie approximately in a straight line. In this case, only rota- 
tion about a single axis is required to keep the gravity correc- 
tion terms small. A coordinate system of this type is shown 
in Fig. 4. 

The XYZ axes are fixed in the body and hence in inertial 
space. The Y axis is chosen perpendicular to a reference great 
circle lying close to the desired path of the vehicle. The plane 
containing the XZ axes is referred to as the guidance plane. 
The orientation of the XZ axes in this plane can be arbitrarily 
chosen, but once defined remain fixed. The position of the 
vehicle is defined by coordinates corresponding to those of an 
ordinary cylindrical coordinate system: (a) a range angle, 0; 
(b) lateral distance from the guidance plane, y; and (c) dis- 
tance from the polar axis of the coordinate system, p. 

The z, y, 2 axes, referred to as true position coordinates, 
differ from the X YZ system by the rotation 6 about the Y-axis, 
and correspond to the desired platform orientation. As the 
vehicle progresses around the guidance plane, the platform is 
rotated through the angle @ so as to remain locally level. As 
the vehicle departs laterally from the guidance plane, the 
gravity correction is computed in exactly the same manner as 
that considered previously for the flat coordinate system. As 
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mentioned earlier, effects due to the rotation of the earth and 
ellipticity are being ignored. While these are important, they 
do not change the essential nature of the problem. 

The angular velocity of the coordinate system with respect 
to inertial space is given by 


[6] 


and the vector position of the vehicle is given by 


Where 1, is being used to define a unit vector along the y-axis. 
Using this expression, accelerations along the x, y, and z axes 


become 


The development can be carried on from this point with 
complete generality. Since, however, the purpose of this 
paper is to illustrate the principles involved, assume that the 
motion is restricted to the surface of the earth and that de- 
parture from the guidance plane is sufficiently small to ignore 
terms containing the factor (y/a)*._ Under these conditions 
the acceleration expressions become 


A; = 
72 
A, = go — ae? — ,-2| 
a ) 


Only the equations for range and lateral accelerations are of 
direct concern, since only two coordinates are required to de- 
fine position for the restricted motion being considered. 

Thus the output of the x or range accelerometer is to be 
doubly integrated and divided by the radius of the earth to give 
the range angle 6. The first integral, 6, can be used to give a 
signal for torquing the gyroscope controlling about the y axis 
to obtain the required angular rate with respect to inertial 
space about that axis. The output of the y or lateral acceler- 
ometer is doubly integrated to give the lateral position y. As 
before, however, a horizontal component of gravity must be 
computed and fed back to the y accelerometer as a correction 
term. A block diagram of the system is shown in Fig. 5. 

Other methods of achieving the required @ rotation are 
available. For example, the gyros could remain inertially fixed 
while one or both accelerometers could be mounted on a 
separate platform which is driven through the angle 6 by a 
mechanical drive. The logical extension of such a method 
leads to the “five-gimballed platforms” discussed by Wrigley 
(4). 


Error Analysis a 


‘Homogeneous Equations 


The guidance system illustrated in Fig. 5 can be used to 
determine position for the rather idealistic situation which is 
being considered. Any actual system, however, will not give 
exactly correct values of position. As a first step in deriving 
the errors in such a system, consider the equations which these 
errors must satisfy, ignoring for the moment the sources of 
these errors. 

Five basic error parameters are needed which can be defined 
in the following way: 

1 Errors in position, A@ and Ay. The outputs of the 
range and lateral channels should indicate the true position 
of the vehicle, (6, y). Defining A@ and Ay as the differences in 


position of the vehicle as indicated by the guidance system 


and true position, these outputs are actually 0 — A@and y — 


Ay. Do Assuming for the moment that the components are per- 
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Fig. 4 Cylindrical coordinate system 
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Fig. 5 System block diagram (cylindrical coordinates) 


2 Platform error angles, $,, ¢,, ¢:. The stable platform 
should be oriented along the z, y, z axes as shown in Fig. 4. 
Again there will actually be errors. Define the z,, y,, 2, axes 
corresponding to actual platform orientation. The trans- 
formation between x, y, z and x,, y,, 2, coordinates can be de- 
fined by means of three Eulerian angles. Since these two co- 
ordinate systems differ from each other by small angles, the 
resulting transformation has the following characteristics: 

(a) The transformation can be defined by giving three angles 
(@z, $,, $2) corresponding to the rotations positive in the 
usual sense about the xyz axes required to go from the xyz 
coordinate system to the z,, y,, 2, coordinates. 

(b) The transformation is independent of the order of these 
rotations. 

(c) The angles ¢;, ¢,, @. may be regarded as components of a 
vector ¢. 

The equation described by propagation of errors in the 
autonavigator can now be derived. The inertial reaction force 
in the x, direction which will be sensed by the range ac- 
celerometer is given by 


Az, = Ar + — = 6 + GA, — [10] 
and the acceleration in the y, direction by 


Ayp = Ay — + = — + [11] 
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Fig. 6 System block diagram including error parameters 
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Fig.7 Accelerometer 


fectly mechanized, a block diagram of the system including the 
error parameters is given in Fig. 6. 

From the block diagram, two of the error equations can be 
obtained immediately 


= f f [a6 + — $,A,.]dr'dr = 
a 0 0 


t T 
0 0 a a 


(y + Ay)]dr'dr = y + Ay 


6+ 


The other three equations are obtained by equating two 
expressions for the angular velocities of the platform. One 
expression for these angular velocities, which follows directly 
from the definition of the ¢ angles, is 


_ where @, is the angular velocity of the platform with respect 
to inertial space and @, is the angular velocity of the xyz co- 
_ ordinate system with respect to inertial space. The second 
- expression for angular velocity of the platform is obtained 
_ from the system block diagram which shows that the only 
_ angular rate applied to the platform is about the y,-axis and 
equal in magnitude to the integral of the output of the range 
- mechanization which may be used to produce this rotation. 
The second equation for the angular velocity of the platform 
becomes 


@, = 1,0 + Ad)................. [14] 
Expressing 1,, in the zyz coordinate system, this equation be- 
comes 

(—dilz + 1, + $:1,)(6 + A6) 


The resulting error equations can be simplified by introduc- 
ing W angles defined as 


‘These correspond to platform misalignment from the position 
6 — AQ@, as defined by the computer rather than true posi- 
tion 6. Since platform torques are based upon computer 
position rather than true position, the y angles more properly 
_ measure performance of the platform components. For more 

complicated coordinate systems a similar transformation can 
7 be carried out and is essential to obtaining manageable 


equations. Substituting these relationships, the error equa- 
tions become 


aAé + (aA@) = — WyA,.......... [17] 


bg + Ay = + [18] 

[19] 

¥.+ by. =0...... 21] 


Note that the y equations can be solved independently, which 
immediately implies that the behavior of the platform as de- 
fined by W angles is independent of any position errors. 


Introduction of Mechanization Errors 


In the previous section, five differential equations for the 
basic error parameters were derived. Since perfect compo- 
nents were assumed, there are no driving functions. The solu- 
tion will therefore give the propagation of initial errors in an 
otherwise perfect guidance system. It is now necessary to re- 
place the idealized components with reasonable approximations 
to the actual instruments. Again, since the purpose of the re- 
port is to describe the principles of inertial navigation, only a 
limited number of the important component errors will be 
considered. For the analysis of any actual system, the list 
would be much more extensive. 

1 Acceleration errors. Until this point it has been assumed 
that the accelerometer gave a perfectly accurate measure of 
the inertial reaction force acting along the sensitive axis of the 
instrument. For any actual instrument, this could not be 
the case. A reasonable approximation to many instruments is 
indicated in Fig. 7. There are two errors in the output: a bias 
error 6A independent of the imposed acceleration and a scale 
factor error ka which depends linearly on acceleration. 

2 Stabilization errors. In deriving the error equation, the 
angular velocity of the platform with respect to inertial space 
was assumed to be exactly determined by the control torques 
applied to the gyroscopes controlling about the xyz axes. 
For the particular case considered, these rates were 0, 6 + A6, 
and 0. Actually there will be uncertainty torques acting on 
the gyroscopes which will produce platform drift. Defining 
the drift rate of a particular gyroscope as being e, the drift 
rate of the platform about zyz axes becomes e,, 6 + A6é + 
€,, and 

3 Initial errors. In addition to component errors, there 
will be errors in initial position A@, Ayo; errors in initial 
velocity 46, Ayo; and errors in initial platform orientation 
and 

When mechanization errors of the accelerometer are included 
the output of the range accelerometer instead of being given 
by Equation [10] becomes equal to 


(1 + kz)Az, + 6Az = (1 + kz) (08 
+ gAy — byAs) + [22] 
and the output of the lateral accelerometer by 
a + y — + 6.4.) 


It is straightforward to repeat the process which led to the 
position errcr equations which now become 


(1 + ky)Ay, + 6A, 
+ 6A,.... [23] 


a 
Ay = + Weds + kyAy + Sy... [25] 
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In deriving the platform error equations, it is only necessary 
to replace Equation [14] by the expression 


= Iz,€: + 1,,(6 + Ab + €,) + 


This can again be equated to Equation [17] and the platform 
error equations become 


These are the complete error equations; their solution with 
the proper initial conditions defines the behavior of an inertial 
guidance system for the restricted case considered here. 


Solution of Error Equations 


The most significant characteristic of errors in inertial 
navigation can be seen from Equations [24 and 25]. These 
show that the position errors satisfy the equation of an un- 
damped harmonic oscillator of period T = 2x(a/g)!/2 = 
S4.4 min, assuming values of a and g at the earth’s equator. 
This is the so-called Schuler period, named after its discoverer 
who explained its significance for acceleration-sensitive 
inertial systems operated on the earth (8). All of the physical 
intuition which has been developed from considerations of 
simple pendulums immediately applies. For example, an im- 
pulse driving function would give rise to undamped sinusoidal 
oscillations. 

The only difficulty in obtaining a complete solution to 
these equations rests in assumptions on the path of the 
vehicle and the form of the driving functions. The following 
set of assumptions form a reasonable approximation to many 
cases of interest and lead to simple exact solutions. 

1 All mechanization errors are constant. 

2 Motion is in a straight line along the guidance plane. 
Velocity is constant except for a rapid change at ¢t = 0, and is 
small in comparison to (ago)'/ or satellite velocity. This as- 


sumption is equivalent to rT * 
Az = v0(t) 
A, = 0 


v/a< (go/a)'/2 


Defining w as (go/a)', Ax = aAé, Az = aA, the solution to 
the error equation is 


Ax . 
Ar = Az cos wt + — sin wt + ae,\ ¢ — — sinwt} + 
w w 


6A, 
w? 


aYyo(1 — cos wt) + kz = sin wt + (1 — cos wt)... . [29] 


Ayo . a . ut 
Ay = Ayo cos wt + — sin wt + aeéz ( - sin — — — sinwt 
v a w 
a? ut 
1 — cos — ~, (1 — cos wt) + 


vl . ul 
00s — cos wt) + sin — 4+ 
a a 


6A 
(1 — cos wt)... . [30} 
w 
t 
a 2 
v a 
[32] 


vt 
a 


aé, . wu 


vt 
Yo cos — + sin — sin — + 
a v a 


Ar Az € 
dy = — cos wt + =) sin wt + cos wt + 
a aw w 


6A, 
kz — sin wt + — (1 — cosw)..... [34] 
aw g 


It is straightforward to obtain the error from any particular 
source by substitution in these equations. The general 
characteristics of some of these errors are probably of more in- 
terest than the exact form. 

1 Position errors. In general the position errors are 
sinusoidal in form with oscillations having the characteristic 
84 min period. Acceleration impulse errors lead to oscilla- 
tions about zero; acceleration bias errors give rise to oscilla- 
tions about an offset position error equal to the product of the 
radius of the earth and the ratio of the acceleration error to the 
acceleration of gravity; gyro drift produces an error increasing 
with time which oscillates about a value equal to the product 
of the radius of the earth and the angle through which the 
gyro has drifted. 

The extreme precision required of both accelerometers and 
gyroscopes is made apparent by the presence of the large 
quantity a, the radius of the earth, as a factor. It should be 
noted, however, that the error is always less, and for large 
times very much less, than an error of the form 1/2 (6A)t? 
which might have been expected. 

2 Platform alignment errors. For some applications, plat- 
form error angles are of considerable interest in themselves. 
It should be noted that gyro drift produces oscillatory plat- 
form errors rather than ones which increase continuously with 
time as might be expected. 

The extension of this analysis to include other important 
sources of error can be made using these same techniques. In 
some cases the results are quite obvious. For example, a bias 
error 6C in the first integrator produces oscillatory errors of 
exactly the same form as the acceleration bias error 6A. It 
has been stated that a system involving the double integration 
of acceleration, such as that considered here, is subject to 
position errors increasing parabolically with time from compu- 
ter bias errors and that the mechanical complexity of a “five- 
gimbal” system is required if high accuracy over long time is 
to be achieved. This conclusion was apparently reached based 
upon a system in which the integrators for each channel were 
repeated, once to determine the feedback terms and again for 
indicating position. The above analysis and confirming system 
results show that it is not true for the case considered here. 


Conclusion 


— might be expected, the sources of error in inertial 


guidance are seen to be exclusively of dynamical or of trigono- 
metrical nature, as distinguished from such sources of error as 
refraction of radiant energy, distortion of magnetic fields, etc. 
This fact leads to the supposition—borne out so far by ex- 
periment—that improvements in inertial instruments, in 
accuracy of initial alignments, etc., will pay off directly in 
enhanced guidance accuracy, without any foreseeable limit. In 
general, the gyroscope presents the more formidable problem, 
since defects thereof can result in cumulative position error; 
defects associated with accelerometers and integrators can 
usually be caused to result merely in bounded errors, due to 
the oscillatory nature of the error functions. In general, also, 
the nature of errors, in a properly designed inertial guidance 
system using the same basic components, is largely inde- 
pendent of the particular mode of mechanization or system 
of coordinates used. A similar analysis to that given here can 
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be carried out for any practical coordinate system including nautical Engineering Review, vol. 15, no. 1, Jan. 1956, p. 344. 
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portant. In some cases RMS time averages of the trigono- 290 
metric functions are often used. More exact solutions of the 5 Slater, J. M. “Measurement and Integration of Accelera- 
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Technical Notes. 


Bore-Surface Temperature Variation 
During Rapid Firing of a 
4.0-mm Gun 


Introduction 


URING rapid fire of an automatic gun, heating of the 
barrel by hot propellant gases is a very severe problem. 


ant It is well known that an extended burst will cause the bore- 

- <I : W. H. GIEDT: AND D. L. RALL? surface temperature to reach the mann point of the gun 
; - metal, and severe erosion will occur. In addition to this 
Detroit Controls, Redwood City, Calif. effect, distortion of the barrel, structural weakening of the gun 


components and the danger of auto-ignition of ammunition 
Experimental bore-surface and inwall temperature (‘‘cook-off”’) may result. Although consideration is given to 
measurements at five locations along a 40-mm gun barrel these effects in gun design, incomplete knowledge of heat 
during a burst firing of eight rounds are presented. Values transfer phenomena makes it almost impossible to proceed 
for heat transfer from the propellant gases to the barrel with complete confidence. 
were calculated as described in a previous paper (1).’ Re- The recent development of surface thermocouples (1-2), 
sults indicate that heat flow into the barrel is essentially has provided a very useful tool for studying this problem. 


: complete about 120 millisec after ignition, and heat flow 
4 


out before the firing of the subsequent round is negligible. 

The total heat transfer per round decreases with each suc- 

Controls Div. of American Standard. Associate Professor of 

of a simple power function. Use of this power function to 2 Research Engineer, Research Department, Detroit Controls 

estimate the minimum surface temperature between Div. of American Standard. 

rounds is also demonstrated. 3’ Numbers in parentheses indicate References at end of paper. 


Received Aug. 5, 1957. 
1Senior Research Engineer, Research Department, Detroit 


BORE - SURFACE THERMOCOUPLES a = 
i 8900" i - * 
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OSCILLOSCOPE 4 FIRING SWITCH 


: Fig. 1 Instrumentation of 40-mm automatic gun 


Epitor’s Note: This section of Jer Proputsion is open to short manuscripts describing new developments or offering comments on 
papers previously published. Such manuscripts are published without editorial review, usually within two months of the date of receipt. 
Requirements as to style are the same as for regular contributions (see masthead page of this issue). 
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Such devices were utilized in the work described herein‘ to 
measure the bore-surface temperature variation during burst 
firings in a 40-mm rifled-barrel gun. The nature of the heat 
transfer variation was then investigated, using these records 
in conjunction with simultaneous inwall (barrel) temperature 
records. 


rw 
Experimental Equipment and Procedure 


40-mm Automatic Gun Instrumentation a 

\ standard antiaircraft gun was instrumented, as shown in 
Fig. 1, with five bore-surface thermocouples along the barrel. 
These units indicate the temperature variation at the interface 
between a thin plated layer of nickel and the barrel steel at a 
distance of approximately 0.0002 in. below the bore surface. 
To supplement bore-surface temperature measurements, an 
inwall thermocouple was installed a depth of 0.038 in. from 
the bore surface and displaced 90 deg circumferentially 
from the No. 3 surface thermocouple. The junction was 
formed by percussion-welding an insulated 0.020-in.-diam con- 
stantan wire to the bottom of a small hole drilled radially into 
the barrel. A thin layer of insulation on the wire was pro- 
vided by oxidizing it in a low temperature flame. 

The outputs of the five surface thermocouples were indi- 
cated by a Maryland Electronics Corporation six-channel 
oscillograph and recorded by a 35-mm strip-film drum camera. 
The output from the inwall thermocouple was registered on a 
Sanborn recorder. Velocities were determined by using mag- 
netized projectiles together with two appropriately spaced 
magnetic pick-up coils. 


Experimental Results 


Bore-surface temperature records were obtained during an 
8-round burst (890-gm projectiles with a 293-gm charge of 
SPDN 7184 propellant, fired at the rate of 115 rounds per 
min). Inwall temperatures were also recorded 0.038 in. from 
the bore surface at the No. 3 thermocouple location. 

An approximate transcription of these bore-surface and in- 
wall temperature variations at the No. 3 location, showing the 
maximum and minimum values, is presented in Fig. 2. Al- 
though a general increasing trend is apparent, a rather large 
variation in maximum temperature occurred, a fact that can 
be attributed to differences in the burning of individual 
rounds. The minimum between-round bore-surface and in- 
wall temperatures, on the other hand, rise consistently, since 
they depend on the total heat transferred to the barrel with 
each round, which is not so markedly influenced by variations 
from round to round as is the peak temperature. 

Fig. 3 shows the distribution of the maximum and mini- 
mum bore-surface temperatures along the barrel measured 
during the firing. Attention is called to the results for Round 
No. 1 to emphasize the extent of burning variation that can 
occur. It appears that burning was retarded initially, result- 
ing in a maximum bore-surface temperature of only 430 C 
ut the No. 1 thermocouple location (as compared to around 
900 C for successive rounds), and higher temperatures at the 
No. 2 and 8 locations. It is interesting that in spite of the 
apparent combustion variation, the muzzle velocity was only 
slightly lower than for the other rounds. Fig. 2 also emphasizes 
the significantly higher surface temperatures reached in the 
region near the origin of rifling. This is the primary reason 
for the high erosion rates observed in this area. 


Thermal Properties Based on Comparative Study of 


Bore-Surface and Inwall Temperature Measurements 

It is generally permissible, in heat transfer calculations, to 
consider a gun barrel as a semi-infinite solid and to neglect 
heat conduction parallel to the bore surface. 


* Performed for the Bureau of Ordnance of the U. 8. Depart- 
ment of the Navy under Task II of Contract NOrd 12102. Their 
— in and approval of its presentation are gratefully acknowl- 

ged. 
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Temperature 


variation with time at planes of the order of several hundredths 
of an inch from the surface can be satisfactorily approximated 
by assuming that the heat transfer from the propellant gas 
occurs instantaneously as a pulse of heat, Q (3). With average 
values for the thermal constants assumed, the temperature 
variation is given for this case by 
Q 
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Fig. 2 Bore-surface and inwall temperature variations during 

8-round burst in a 40-mm gun (maximum and minimum tem- 


peratures indicated) 
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Fig. 3 Maximum and minimum bore-surface temperature 
distribution along a rifled 40-mm barrel during an 8-round burst 


INWALL THERMOCOUPLE LOCATION: 
0.038 INCH FROM BORE SURFACE wu 
27.5 INCHES FORWARD FROM ORIGIN OF RIFLING 
X TOTAL HEAT TRANFERRED- EQUATION | 
re 
| | 
t + t t 
°. MEASURED 
O__ ____ PREDICTED - EQUATION 2 
SING GORE- SURFACE TEMPERATURES 
HEAT INPUT- EQUATION 5 
—+e 
} 
Ls 
| 
re 
a: os os os oe or 
Time - Seconds 


Fig. 4 Inwall temperature variation vs. time during first round 


of an 8-round burst in a 40-mm gun 
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Fig. 5 Heat transfer to a 40-mm gun barrel 1'/; in. from origin 
of rifling during an 8-round burst 


Note that if the temperature 7 at any time ¢ and distance 
from the bore surface z is measured, this equation can be 
used to calculate the total heat transfer Q (per unit area) 
assuming average values for the density (p), specific heat 
(c), thermal conductivity (4), and thermal diffusivity (a). 

A more precise expression for the temperature variation 
with time and distance is (3) 


tp(A)e7 7/42 dr 
(t — 


where the function $(A) represents the temperature variation 
with time at the bore surface, \ being the time variable of 
integration. This equation is applicable for all values of z. 
By differentiation of Equation [2] with respect to xz and inte- 
gration with respect to time at z = 0, the heat transfer up to 
any time t is found to be (1) 


De (t — 
In using any of the above equations, it is necessary to know 
the appropriate values of the thermal conductivity (4) and 
diffusivity (a). Fortunately, an effective value of a can be de- 
termined from a measured inwall temperature record and 
Equation [1]. Differentiating this equation with respect to 
time and equating the result to zero gives the time at which 
the temperature at any distance z reaches a maximum as 


By using this expression, a value for a = 0.0352 cm2/sec was 
determined. 

A check on the validity of this result was possible by using 
it in Equation [2] along with the measured surface-temperature 
history to predict the temperature variation at the location of 
the inwall thermocouple. This result is compared in Fig. 4 
with the experimental curve. The agreement was regarded 
as reasonably satisfactory, being within about 10 per cent at 
the maximum temperature point. Slight adjustments in the 
value of a did not yield improved results, suggesting that 
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variations with temperature or possibly other factors may be 
involved. 

For calculation of the total heat transfer to the barrel, an 
average value of k (0.0500 cal/sec-em-C) was determined 
from the above value of a, using values of p = 7.85 gm/cm? 
and c = 0.181 cal/gm-C (4) . Substituting these values and 
the experimental data at the time of maximum inwall tem- 
perature in Equation [1] gave 


Q = 6.8 cal/em? (at No. 3 thermocouple location) 


Analternative calculation using Equation [3], and substituting 
the experimental surface temperature variation, yielded a 
value of 


Q = 6.9 cal/em? 


It is significant to note that the effective value of therm: 
conductivity determined and nominally confirmed by the 
above calculations is only of the order of 60 per cent of that 
indicated for the barrel metal in the literature (4). This 
difference may be due to material variations but is more 
likely attributable to the altered metal layer which is eventu- 
ally formed at the bore surface by the combined chemical an‘| 
thermal effects of the hot propellant gases. . on 

Heat Transfer During Rapid Fire -_ 

To investigate the nature of the surface heat transfer varia- 
tion with time during a burst firing, measured temperatures at 
the No. 1 thermocouple location (14 in. from the origin of rif!- 
ing) weresubstituted in Equation [3]and the integration carrie: 
out numerically. Average values of k and a (0.050 cal/sec- 
em-C and 0.0352 cm?2/sec, respectively), determined as de- 
scribed above, were used. The results have been plotted in 
Fig. 5, where it is shown that heat from the propellant gas 
flows to the barrel for a period of from 0.1 to 0.2 sec after « 
round has been fired. The heat rate decreases very rapidly, 
however, so that most of the transfer occurs soon after the 
firing. Thus, the curve has the appearance of a series of steps. 
The amount of heat transfer per round can be clearly observed 
to decrease with each successive round in the burst, and afte: 
the second round some reverse flow actually occurs. 


Prediction of Minimum Surface Temperature 


Since it may be of interest to know what the bore- and com- 
bustion-chamber surface temperatures will be near the be- 
ginning of any round in a burst (to estimate danger of cook- 
off), consideration was given to using the information in Fig. 
5 to predict values for these temperatures. It was first 
noted that if the bore-surface temperature variation ¢(é) is 
unknown and Q(t) is known, expression [3] may be regarded 
as an Abel integral equation. Its solution (5 )is 


sin vr d Q(A) 


dt (O<v<1).... [5] 


g(t) = 


where vy is the exponent of (¢ — A) in the denominator of the 
integrand of Equation [3]. With interest focused on the sur- 
face temperature just before the firing of each round, the possi- 
bility of approximating Q(t) by a simple power function of ¢ 
was investigated. Because of the decreasing heat transfer per 
round in a burst, an expression of the form 


[6] 


with n < 1, was found to give a curve passing through the 
points on the actual curve corresponding to the firing of each 
round (the dashed curve in Fig. 5). The surface temperature 
variation resulting from a heat input in accordance with 
Equation [6] might therefore be expected to coincide with 
the actual surface temperature just prior to the firing of 
each round. 

If the expression for Q given by [6] is substituted in [5] and 
the integrand expanded in a series, integration yields 
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For values of n < 1, this series converges rapidly. Use of only 
the four terms indicated was estimated to give an error of less 
than 3 per cent. 

The applicable values of C and n were determined from a 
log-log plot of heat transferred per round vs. round number. 
The equation of the curve through the experimental points 
for the No. 1 thermocouple location was found to be 


where K is the heat transferred during the first round (23 
cal/em?), NV the round number, and the exponent n = 0.837. 
If R denotes the rate of fire, VN = Rt, and Equation [8] be- 
comes 


Q(t) = KR" = [9] 


The rise in bore-surface temperature is then : 


/ 
T, = o(t) = 24.4 Rosse [10] 

The bore-surface temperature at the No. 1 thermocouple 
location at the end of the 8-round burst, as calculated from 
Equation [10] was approximately 15 per cent higher than the 
experimentally observed value. In view of the assumptions 
involved, this sort of agreement may be regarded as satis- 
factory. 

In applying Equation [10] to other locations along the bar- 
rel, the appropriate value for the heat input during the first 
round (K) must of course be used. Exponent n (Equation [8]) 
would, however, be expected to depend primarily on the 
thermal properties of the barrel and should be reasonably 
constant for a given gun. A limited check on this hypothesis 
was obtained by substituting the value of K = 8.3 cal/cm? 
determined for the No. 3 thermocouple location in Equation 
[10], using the value of n found for the No. 1 location. This 
vielded a value of 140 C at the end of the 8-round burst, which 
compares very favorably with the measured value of 136 C. 
An additional check was made using the bore-surface and in- 
wall temperature data given in (2) for a 5-round burst in a 
0.50 caliber machine gun. Agreement between predicted and 
observed results was also within 15 per cent. 
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The effect of extent of completion of a radical recombina- 
tion reaction on the specific impulse of a simple system 
has been considered. The kinetic impulse between the two 
conventional limiting conditions of frozen and equilibrium 
flow has been obtained as a function of recombination rate 
constant, chamber pressure and motor diameter. The 
most recent data (obtained from shock tube measure- 
ments up to 2800 K) on recombination rate constants are 
used. These values are smaller by a factor of 10~-10* than 
those previously considered in this connection. For the 
best value of the rate constant, 20-80 per cent of the dif- 
ference between frozen and equilibrium flow is realized, 
depending on the motor size and chamber pressure. 


N CONVENTIONAL propellant thermochemical calcula- 
tions, specific impulses are cited for both frozen and equi- 
librium flow. These are the two limiting conditions, cor- 
responding to zero and infinite reaction rates, respectively, in 
the nozzle expansion. Chemical rate constants are, of course, 
never either zero or infinite and the true kinetic impulse must 
lie between the two limiting conditions. Penner and Altman 
have given this problem extensive consideration and have 
derived approximate criteria for the realization of near-frozen 
and near-equilibrium flow (1).2- It appears desirable to study 
the actual kinetic impulse as a function of rocket and kinetic 
parameters for a simple model. Krieger (2) has carried out an 
exact numerical solution of hydrogen gas expanding in a 
nozzle, starting at 3500 K (about 12 per cent dissociation), 
for a particular set of conditions. He found that about 85 per 
cent of the difference between frozen and equilibrium flow 
was realized for a motor diameter of 10 cm, a recombination 
rate constant of 10liter?/mole? sec, and a chamber pressure 
of 300 psi. 

Another use for this study of nozzle kinetics is the observa- 
tion of the effect of the new data on recombination rate con- 
stants (3-5). Penner and Krieger used recombination rate 
constants in the range 10'°-10"'liter?/mole? sec, which were 
derived from rather old room-temperature data. The later 
values (for recombination of iodine and bromine) were ob- 
tained in a shock tube up to 2800 K, and were not only smaller 
than previous values, but also showed a decrease with tem- 
perature. 

The present precision of the data is only sufficient to indi- 
cate that the rate constant is inversely proportional to the 
first or perhaps three-halves power of the temperature. A 
number of rudimentary theoretical treatments of the recom- 
bination problem have been made (6-8), generally indicating 
an inverse one-half power dependence of rate constant with 
temperature. Recombination rates are larger for larger third 
bodies, such as benzene, but only small molecules and atoms 
are available at flame temperatures. All things considered, 
the best order of magnitude of atomic recombination rate 
constants under nozzle conditions (2000-3000 K) now appears 
to be 10%liter?/mole? sec, a factor of 10?-10* smaller than pre- 
viously thought. 

Even if we restrict the problem to a single recombination 
reaction and make other noncritical simplifications, such as 
constant heat capacities and zero chamber velocity, a numeri- 
cal solution is still necessary. An approximate closed analyt- 
ical solution can be obtained if, following Penner, we use a con- 
stant average value of the cooling rate in the nozzle, —dT /dt. 
The result for the true kinetic impulse should still reflect the 
variation with the three parameters under study: Rate con- 
stant, chamber pressure and motor diameter. In order to find 
the true impulse we must find the gas composition at 7’, the 
exhaust temperature; that is, we must assess the extent to 
which reaction [1] proceeds 


2 Numbers in parentheses indicate References at end of paper. 
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_ where R is a radical or atom, X is a third body, and k, and k, 
are the forward and back rate constants, respectively. The 
___ parameters which are fixed at typical values in this study are 


3500° K 
P,/P, = 1/68 
chamber velocity = zero nay 


Do = dissociation energy of R. = 100 kcal 
Ye = mole fraction radicals at 7’, = 0.10 
AS —23 cal/mole® K 


C, = 8cal/mole® K 


The general plan of approach is to find the exhaust tem- 
perature and composition from the constant entropy condi- 
tion and the relation provided by the kinetics of reaction [1]. 
The specific impulses for frozen, equilibrium and kinetic flow 
are then calculated by standard formulas. The dependent 
variable will be a reduced impulse J, defined gnete 


nek, 


1, = Tirue Ttrozen 
Tequil = Ttrosen 
which expresses the fractional difference realized between 
frozen and equilibrium flow. Thus, when y, = y, (frozen flow), 
I, = 0; and when y, > y, (equilibrium flow), J, = 1. A plot 
of I, vs. y, — y, is nearly linear. 
Returning to the kinetic problem, the time rate of change of 
concentration of R is given by (1) 


rex | | 3] 
da 


where R, R:, and X are the concentrations of the respective 
species in mole/cc, K, = k,/k, = R,TK, is a concentration 
equilibrium constant and y is a suitable mean heat capacity 
ratio. The first term on the right expresses the change in R 
due to expansion. It disappears when Equation [3] is trans- 
formed to mole fractions 

E p 


Equation [4] may be solved analytically by using certain 
approximations which are at least as good as the assumption 
of a constant average cooling rate. The details will not be 
given here, but the result is an equation relating exhaust tem- 
perature and composition. With the aid of Equation [1] the 
results in Figs. 1, 2 and 3 were derived. 

An inspection of the figures shows that with the new order 
of magnitude of the recombination rate constant k, (10%) the 
recombination is not complete, but that 20-80 per cent of the 
difference between frozen and equilibrium flow is realized, de- 
pending on the parameters. For example, with a motor 
diameter of 60 cm and a recombination rate constant of 10°, 
67 per cent of the difference in impulse between equilib- 
rium and frozen flow will be realized at 300-psi pressure and 
86 per cent will be realized at 1000-psi pressure. However, 
if the recombination rate constant is 10’, only 20 per cent 
of the difference in impulse will be realized at 300-psi and 65 
per cent at 1000-psi pressure. 
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Transfer Between Vehicles in Circular | 


Orbits 


BERNARD H. PAIEWONSKY! 


Wright Air Development Center, Ohio 
Nomenclature 
A = initial vehicle 
B = terminal vehicle 
y = gravitational constant 
M, = mass of the earth 
Rx, = radius of circular orbit of vehicle A 
Ry = radius of circular orbit of vehicle B 
ws = angular velocity of A 
wp = angular velocity of B . 
ry = time of flight in Hohmann orbit 
= period of Hohmann orbit 
@ = central angle between vehicles A and B 
Angles ABCDE: 
Angle D: angle between true vertical and horizon as seen by 
vehicle A 

Angle E: angle between horizon and terminal vehicle as seen by 


initial vehicle A 
ugle Ho: angle E at time of launch 


Introduction 


— purpose of this investigation is to examine the orbital 
transfer of a commuter vehicle between two co-planar 
satellites. Suppose we wish to go from one satellite vehicle to 
another by means of a commuter rocket using chemical pro- 
pellants. In order to conserve fuel, an orbit requiring the 
minimum expenditure of impulse should be sought. The op- 
timum paths for transfer between circular orbits have been 
found by Hohmann (1)? to consist of the family of cotangential 
ellipses. In practice, it will usually be necessary to go from 
one vehicle to another, rather than merely changing orbits. 
In order to do this without applying corrective thrust in 
flight. the launching of the commuter rocket from the initial 
orbit must be timed so that the commuter arrives at the 
terminal orbit coincident with the second orbital vehicle. D.F. 
Lawden (2, 3) has derived in general terms the equations for 
orbital transfer between two planets. This note presents a 
simple approach to the restricted case of co-planar satellites 
in concentric circular orbits. 


Transfer Between Satellite Vehicles - 

The analysis presented here considers only Hohmann trans- 

fer orbits requiring minimum expenditure of fuel. In order to 

use a Hohmann orbit (a cotangential ellipse) there is a definite 
angular relation that must exist between the initial and 


terminal vehicles. This relationship is based on the require- _ 


ment that the orbital vehicle and the commuting vehicle 
meet. A lead-collision calculation is required. 

Let vehicle A be the initial vehicle and let vehicle B be the 
terminal vehicle. For convenience, assume a geocentric polar 
coordinate system (r, 6) with 6 = 0 taken through the initial 
vehicle at time t = 0. Let Ra, Rx be radii of orbits of vehicles 
A and B respectively. The case of Ra = Ry will not be con- 
sidered. The angular velocities w, and wx of the satellites in 
these orbits are given by 


VyMe 


=VyM. 


Ty in the Hohmann orbit, from apogee 


WA 


WB 


The time of flight, 
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Hohmann orbit showing positions of initial and terminal 
vehicles 


Fig. 1 


to perigee is given by half the period, 7, of the Hohmann 


orbit. 


Ty 
H 2 


where (Ra + Rx)/2 is the semimajor axis of the orbit of 
transfer (4). 

The commuter vehicle will reach the terminal orbit at 6 = 
am. The problem now is to find a position of vehicle B at the 
time of launch such that the commuter vehicle from A will 
meet it at @ = 7. In the time ¢, vehicle B has gone an angular 
distance A@ given by 


=wpt 


where Ax is taken congruent modulo 27.3 
From Equation [2] we know that 


2 VyMe 


Thus 


= 


2 


The angular location (@3),.0 of vehicle B at t = 0 can be 
found by subtracting from 


3/3 
(On) = 7 — = (1 ) [4] 


~The geometry of the maneuver is shown in Fig. 1. 

If we wish to go from an outer orbit to an inner orbit and if 

Ry is very much greater than Ra, we must take A@ congruent 

modulo 27 and count the number of revolutions of the vehicle 

in the inner orbit. That is, if the time of transfer exceeds the 

period of the inner or terminal vehicle, integral multiples of 

the period must be subtracted from the transfer time. The 

remainder ¢ is then used to compute A@ and to determine the 


3 This means that integral multiples of 27 are subtracted from 
wpl. The final remainder, Aég, is thus < 27. 
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proper initial angular relation between the two satellites at the 
time of launch from Equation [4]. 


In-Flight Determination of Angular Relationships 
Between Satellites 


It will be necessary for the commuter rocket to know the 
angle between the terminal vehicle and the initial vehicle 
el om to launching. This information may be obtained by 

communication with an outside observer or by direct meas- 
urements. By means of measurements of the angle between 
the horizon and the terminal vehicle, and the known quanti- 
ties, Ra, Rg, and (@),-0, the time to launch may be determined. 

Tt may happen that the required angular relationship be- 
tween the vehicles is such that visual observations cannot be 
ok = This method works only when there is an unob- 
___ structed line of sight between the observers. Fig. 2 illustrates 
measurements required. 

a _ Formulas from elementary trigonometry provide the solu- 
tion. +The lettered angles are shown in Fig. 2. iene 


SnB SnC Sn(D+EZ) 


Rz sin 63 


. D E = 
sin (D + E) + Rp? — 2RaRz cos 


E = (sin! K) — D 


The value of angle Z at the time of launch can be calculated 
and is called Ep. 

When the angle between the horizon and the terminal 
vehicle is determined to be equal to Ey the separation angle 
between the vehicles is the required 6,3 and the commuter 
rocket should be launched.‘ 

The effects of errors in measurements and navigation are not 
considered in this paper although they certainly will be of im- 
portance in the practice of space travel. 


Results 


The relationship between 6, and Rs/R, is illustrated in 
Fig. 3. As Rs/R, gets very large it is seen that 63/7 ap- 
proaches a maximum value of 0.646. Although the case 
Rx = Rs was not considered in the derivation, the 63 curve 
goes through zero when this occurs. This means that if the 
satellites have the same initial radii then they must coincide, 
otherwise one cannot travel between them using only one 
Hohmann orbit. 

These results will now be applied to some problems. For 
a trip between two artificial earth satellites with radii Ra = 
R, + 500 km and Rs = R, + 1000 km, R, = 6370 km. Fig. 3 
shows 63 = 9.08 deg. 

The conditions necessary at the start of a trip between 
space station staging satellite and the moon are also obtained 
from Fig. 3. If the space station R, is 7000 km from the 
earth’s center and the moon is assumed to be in a circular orbit 
with radius Rg = 3.88 X 10°km, then 63 is found to be 115 
deg. The time of flight for this trip is 


Ty = ———— = 122 hours 
/yM | 2 


or a little over five days. 

This method can also be used for rough calculations for in- 
terplanetary travel. Suppose we are planning a trip between 
the earth and mars, using chemical propellants. The elliptical 
planetary orbits are approximated by circular orbits about the 


sun. 


The initial angle between the planets, 6, is seen to be 
nearly 45 deg. The time of flight for this journey is 264 days. 
The time of flight for this trip using 1-g acceleration has been 
computed by Jones (5) to be 2.08 days. The most economical 
flight path is not the shortest in time of travel. 

A return trip from the Moon to the space station circling 
the earth may be examined using Equation [4]. Using the 
radii given in the previous example we can obtain 6,’ and 

Let the staging vehicle be at a radius Rs = 7000km. The 
moon is assumed to be at radius Rg = 3.88 X 105 km. 


Ra + Rs |”? 


d Aés = (mod 27) 
(84) = 7 — = —155.9° ~ —156° 


Ra (earth) = 1.41 X 108 km 
Rs (mars) = 2.28 X 108 km 


Ads’ = 0.93(2r) = 335.9° 


‘It has been pointed out by a reviewer that the angle Ey may 


correspond to two distinct values of 6g in the case of transfer from 
a larger to a smaller orbit. A cure for this ambiguity is to measure 
the distance between vehicles A and B. The correct range has 
already been computed and is equal to a, with angle A = 6Breq. 
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Another question may be raised at this point. Suppose 
vehicles A and B are initially in arbitrary positions with 
separation angle e«. How long must wait be the 
found to be 


tdelay 


The maximum time delay occurs if (Oreq — €) 


2r 
tmax delay = | 
Wa — 
Summary and Conclusions 


A simple method is developed for calculating the angular 
relationship required between vehicles desiring to use Hoh- 
mann orbits for orbital transfer. This angular relationship is 
computed and is shown to depend only on the radius ratio of 
the orbits. A method is also presented for determining the 
relative angular positions of the satellites from in-flight 

isurements. 
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Some Observations of Flame 
Stabilization in Sudden 
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Introduction 


REMIXED flames stabilized by sudden expansions ex- 

hibit oscillatory burning within the region of stabiliza- 
large. An experimental study of such a configuration was 
tion on the limits of stabilization. 
herein for consideration. 

Combustion instability as applied to afterburners has been 

studied previously by Rogers and Marble (1),? and Noreen 
and Kaskan (2). Although the proposed mechanisms are 
controversial, there seems to be general agreement that the 
oscillations are defined acoustically by the chamber geometry. 
The configuration under investigation exhibited an organ pipe 
type of oscillatory combustion in which the observed fre- 
quency corresponded to the longitudinal resonant frequency 
of the chamber. 
The oscillatory combustion had relatively no effect on the 
lean limits, but had a large influence on the rich limits of 
stabilization. 
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= 27. Thus 


“Die Erreichbarkeit der Himmelskorper,” 


2 Lawden, D. F., Journal of the British Interplanetary Society, — 


tion when the length-to-diameter ratio of the chamber is — 


conducted to determine the effects of the oscillatory combus-_ 
The results are presented — 


A schematic diagram of the variable length sudden expan- 
sion combustion chamber used in this investigation is shown 
in Fig. 1. From this figure it is apparent that the desired 
variation in length is facilitated by the installation of a metallic 
plunger which can be positioned within the Vycor tube com- 
busion chamber. The plunger is mounted by means of an O- 
ring and asbestos seal ring to prevent leakage and heat trans- 
fer within the annular cavity between the concentric cylindri- 
cal sections. Not shown in this figure are two 200-mesh 
screens placed upstream in the feed line to assure adequate 
mixing of the propane with the air. It is noted that the 
expansion ratio can be represented by the quotient A?/B?. 

The frequency spectrum of the sound emitted from the 
exhaust is measured with a General Radio sound level meter 
and analyzer. The corresponding sound wave is determined 
from oscilloscopic recordings, the oscilloscope being used in 
the place of the analyzer in the primary measuring system. 
The wave shape is either saw-tooth sinusoidal or flat, depend- 
ing on whether the burning is oscillatory or smooth. The fre- 
quencies obtained from the wave diagrams are in excellent 
agreement with those measured from the analyzer. 


Qualitative Observations 


7 Fig. 2 is a sketch of the various zones observed with a stable 
flame. The sharp increase in axial pressure near the sudden 
expansion (i.e., note flaring of the flame), along with the 
velocity gradient present at the jet interface, create a recircu- 
lation in the region around the mixing zone near the face of 
the plunger. As in the case of bluff body stabilization this re- 
circulation zone provides the necessary ignition energy to the 
premixed fuel-air entering the mixing zone. The length of 
the flame decreases slightly with increasing Reynolds number, 
but is relatively insensitive to the length of the chamber. How- 
ever, at any given approach velocity, the mixture strength 
greatly affects the length of the flame. 

As the lean limit of combustion is approached, both the 
length and width of the flame decreases until, immediately 
prior to blowout, the flame appears as a faint brush around the 
jet. With rich blowout the length of the flame again de- 
creases, but its width increases, filling the tube. Striations 
appear in the middle of the flame, extending upstream into the 
mixture jet; and the volume of the recirculation zone is 
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Fig. 1 Variable length sudden expansion chambers 
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noticeably reduced. In this condition the flame blows off 
immediately with the incidence of any flame pulsations. 
The flame photographs shown in Fig. 3 illustrate the change 
in the flame configuration which occurs with oscillatory com- 
bustion. The recirculation zone is no longer definable; the 
jet, previously sharply outlined, is now obscured; and the 
whole flame is pushed up against the face of the piston. Ac- 
companying this change is a three- to fourfold increase in 
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static pressure, indicating that flame blockage is significantly 
increased with the onset of oscillatory burning. oy 

The stability limits for the 16-, 18- and 20-in.-length cham- 
bers with expansion ratios of 3.92:1 are presented in Fig. 4. 
Also shown in the figure are the regions of oscillatory burning 
encountered in the 18- and 20-in. chambers; combustion in the 
16-in. chamber was smooth throughout its region of stabiliza- 
tion. The upper portion of the blowout curves could not be 
definitely established because of the rough burning en- 
countered at the high velocities. Rough burning is distinctly 
different from oscillatory burning, the former occurring in- 
termittently at a random frequency and the latter exhibiting 
sustained operation at constant frequency. 

The lean limits of oscillatory burning were obtained by de- 
creasing the fuel-air mixture ratio of an oscillatory system 
until smooth burning is achieved. The data obtained in this 
manner are remarkably consistent and reproducible. The 
rich limit of oscillatory combustion is not as sharply defined, 
however. As the mixture fuel-air ratio is increased, the os- 
cillatory burning gradually changes to either smooth or 
rough burning, depending on the magnitude of the approach 
velocity. The region of oscillatory burning is indicated by 
the cross-hatched area in Fig. 4. 

Further inspection of Fig. 4 indicates that the rich blowoff 
limits for the 18- and 20-in. chambers are coincident through 
the range of velocities in which oscillatory burning occurred, 
and no definite rich limits could be determined at velocities 
greater than 180 fps. In addition, it is noted that the region 
of oscillatory combustion extends up to the stability limits of 
rich burning for the 18- and 20-in. combustors. From these 
observations it is evident that the existence of oscillatory com- 
bustion restricts the range of stable burning and imposes an 
upper bound on the velocities at which rich blowoff can be de- 
termined. It is noted that Polanyi and Markstein (4) re- 
ported similar narrowing of the stability region by the im- 
position of high frequency (3-12.5 ke) sound on a premixed 
flame. 

High speed motion pictures of the oscillatory phenomenon 
in a 20-in. chamber with an expansion ratio of 3.92:1 are 
shown in Fig. 5. These pictures were obtained with a 16-mm 
Fastex camera operating at 3000 frames per sec with the 
necessary luminosity provided by the introduction of mag- 
nesium oxide into the flame. The frame sequence proceeds 
from top to bottom and from left to right. The region photo- 
graphed is that area near the edge of the expansion. These 
pictures clearly show the longitudinal character of the pulsa- 
tions, approximately three cycles appearing in the sequence. 
Nicholson and Field (3) obtained similar pulsation pictures of 
a flame stabilized on a V-gutter which indicate that the charac- 
ter of the pulsations, once they are established, is independent 
of the mode of stabilization. 

Further study revealed that the frequency of oscillatory 
burning is insensitive to expansion ratio while varying in- 
versely with the length of the chamber. The measured fre- 
quencies for the 18-, 20- and 22-in. chambers with expansion 
ratios of 2.85:1 and 3.92:1 were correlated by the equation 


1 


with a maximum deviation of 5 per cent. In this equation, 
which was obtained from a simplified one-dimensional per- 
turbation analysis of the chamber flow process, a and M are 
the speed of sound and Mach number of the hot gases in the 
chamber, respectively; and L is the length of the chamber. A 
similar analysis was performed previously by Woolstan and 
Runyon, the details of which can be found in (3). The cor- 
relation presented by Equation [1] further establishes the 
longitudinal character of the oscillations which were ob- 
served. 


Results 
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Fig. 5 High-speed motion pictures of oscillatory burning; 3000 
140 fps, ¢ = 


frames per sec, , = 18in., V = 1.35 


Conclusions 


1 Oscillatory burning restricts the range of stable burning 
for rich flames and imposes an upper bound on the velocities at 
which rich blowoff can be determined, but it had no effect on 
the lean limits of stabilization. 

2 The lean limit of oscillatory burning is sharply defined 
whereas the rich limit is vague. 

3 The flame blockage is significantly greater with use 
tory combustion. 

4 The oscillations are longitudinal and correspond to 7 
resonant frequency of the chamber. After oscillatory com- 
bustion has been established, the character of the — 
is independent of the mode of stabilization. 
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On Radiation From Combustion Gas _ 


SYOGO MATSUNAGA! 


‘High-Frequency Os- 
ASME, vol. 


University of Osaka Prefecture, Sakai, Osaka, Japan _ 


HEN the temperature of combustion gas in the furnace 

is held constant, combustion gas burning perfectly radi- 
ates weaker than combustion gas burning imperfectly by the 
effect of unburned carbon particles illuminating in the com- 
bustion gas. 

If we burn gas imperfectly for the purpose of obtaining 
more radiation from combustion gas by the effect of illuminat- 
ing carbon particles, then the gas temperature falls by the 
effect of imperfect combustion. 

In this paper, the author presents the method of burning 
gas for obtaining more radiation than the usual combustion 
method, by which we can burn gas perfectly for obtaining 
high temperature combustion gas under the condition of com- 
bustion gas having metallic ion or other ions in gas, and these 
ions give the radiation effect of the illuminating carbon 
particles to perfect combustion gas 

The experimental method can be seen in Fig. 1; the ex- 
periments were performed by adjusting the air flow rates into 
the Bunsen burner for obtaining the gas flame at the condi- 
tion of perfect combustion. For the two cases, i.e., the case of 
free combustion of gas and the case of combustion of gas 
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Radiation from flame 

having no wire im it. © 
Fig. 2. Experimental result 


- colored by the small steel wire held horizontally in the gas 
flame, the author read the temperature scale of the radiation 
meter. 
The experimental results can be seen in Fig. 2, by which we 
can understand the following facts: 
1 The temperature drop caused by the latent heat of 
_ vaporization of metal can be neglected. 
2 The radiation effect of the perfectly burned gas will be 
’ made stronger by adding the ions used for coloring the 
flame. 


~ On the Thermal Resistance of the Water 
Droplet on the Metallic Surface 


SYOGO MATSUNAGA! 


University of Osaka Prefecture, Sakai, Osaka, _— 


= ge THIS note, the author presents the new electric tank 
method for study of heat conduction in solids. The solid 
mixture of paraffin wax and carbon black powder can be used 
— the electric tank filled with city water. 
__ Using this method, the author presents an example. The 
_ thermal resistance of the water droplet on the metallic surface 
ean be obtained by the following process: 
The experimental method is shown in Fig. 1, where 


 ~-R4 = electric resistance (i.e., thermal resistance for the metal 

sheet and water droplet) 

_ q Rp = electric resistance (i.e., thermal resistance of metal sheet) 
- \, = thermal conductivity of metal sheet (given by city water 

in the electric tank) 

Xe = thermal conductivity of water droplet (given by the solid 
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aa mixture of paraffin wax and carbon black powder in the 
electric tank) 

LaLg = distance between two electric poles adjusted its 

position in the water tank. 


Fig. 1 Experimental method ~ 
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Fig. 2 Experimental results obtained by the electric tank method 
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Fig. 3 Thermal resistance of the water droplet on the metallic 
surface 


When we adjust L, and Lg to the balanced state of the 
Wheaston Bridge, the relation of 


Ra La 


can be obtained. 

The experimental results can be seen in Fig. 2 

When fA has the value of nearly zero, the ratio of R;/Ra 

must be nearly equal to the inverse ratio of the contact 7 ace 
i.€. 

Rs 

so the curves in Fig. 2 must have the 

tt h/D = 0. 


area, 
_ — 


= 0.89 
1d? 


= 0.89 
By this fact, the curves of the thermal re- 
sistance can be written as shown in Fig.3. 


value of Rp/R4 


Calculated Viscosity of a Solid 
Propellant Rocket Exhaust Gas Mixture 


W. GIN! 
Hughes Aircraft Company, Culver City, Calif. 


The absolute viscosity of the exhaust gas mixture from 
the rocket combustion of a composite-type solid propel- 
lant has been calculated according to rigorous kinetic 
theory (Reference 2). Seven species including trace con- 
centrations were considered, and the temperature of the 
mixture was taken to be 2160 R. The computed value of 
the viscosity is 9.76 (10-7) Ib-sec/ft?, which is shown to be 
approximated well by a value obtained from less rigorous 
theory. 


Discussion 


HE purpose of this note is to report the result of a calcula- 

tion of the absolute viscosity associated with the laminar 

shearing stress of a composite-type, solid propellant exhaust 
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gas mixture. This mixture is assumed to be under the thermo- 
dynamic conditions prevailing near the exit plane of a super- 
sonic rocket nozzle. There is considerable interest in the heat 
transfer characteristics near the exit of the nozzle because of 
the possible location in that section of various devices like jet 
deflectors or parasitically ignited tracking flares. Although 
the flow conditions may be turbulent rather than laminar in 
the nozzle exit region, the value of viscosity due to molecular 
momentum exchange is nevertheless of practical application 
(1).? 

For this computation, the temperature of the gaseous mix- 
ture is 2160 R, the assumed local static temperature of the 
free stream at the nozzle exit. A total of seven molecular 
species in the concentrations shown in Table 1 was con- 


sidered. Chemical equilibrium is assumed. 

Table 1 Exhaust gas composition of composite- type 
== solid propellant 
C nt Mole fraction 
H.0 (0.2 

HCl (0.23000 
co 
CO, 
Ne 
so, 
Se 
Note: There remains less than 1 per cent of other gases. 


At first, the rigorous kinetic theory of Hirshfelder, Curtiss, 
Bird and Spotz (2) was applied to this calculation. The pro- 
cedure entailed the determination of the force constants for 
each species in like as well as unlike molecular collisions and 
the formulation of a somewhat unwieldy eighth-order deter- 
minant. The evaluation of the determinant was performed 
by the staff with the IBM 650 digital computer. The result- 
ing value for the absolute viscosity of the mixture is 9.76 
Ib-sec/ft?. 

Reference (2) indicates a correction to be applied to the 
viscosity calculation in the event that a strongly polar 
molecule like water occurs as one of the components of the 
mixture. Table 1 shows that H,O turns up as one of the 
major components. The correction is to be made in the 
evaluation of the interaction force constant of water vapor 
with the other constituent gases. However, when the cor- 
rection was made, it was found that the changes on the force 
constants were comparatively insignificant, so that H.O in 
this calculation may be treated as a nonpolar molecule. 

Although there is no published experimental information 
available to the writer for verification of this result, a check 
of some data (3) concerning the viscosity of some of the pure 
component gases indicates that the calculated result for the 
mixture is very credible. As the method of Hirshfelder and 
his associates is based on the Lennard-Jones model of the 
molecule, there may be some question as to the conformance 
of the components of a high temperature gaseous mixture to 
that model. Whalley and Schneider (4) have pointed out that 
the viscosity is not greatly affected by departure from the 
Lennard-Jones model and that extrapolation of the viscosity 
vs. temperature may be accomplished up to 2700 R with 5 per 
cent or better accuracy. It would seem that for the tempera- 
ture range generally encountered in the exit section of rocket 
nozzles, viscosity calculations by this method should lead to 
satisfactory results. 

However, it is realized that the comparatively weak de- 
pendence of heat transfer coefficients on the value of viscosity 
and other usually gross assumptions made in heat transfer 
analyses may not justify the careful and time-consuming 

calculation of viscosity by the rigorous method (assuming that 


2 Numbers i in parentheses indicate References at end of paper. 
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Dynamic Pressure Measurement 
(Continued from page 85) 


able for high temperature operation, but 
are widely used as blast gages.6 They 
were difficult to calibrate statically be- 
cause of the lack of extremely high input 
impedance amplifiers until the recent 
introduction of electrometer tube cir- 
cuits and dielectrics which make possi- 
ble the retention of charges for static 
calibration. 

The research work done by Brosa, 
Nier and Mauer in Germany has led to 
the development at Freiburg of several] 
successful models of water-cooled and 
air-cooled quartz pickups (16). The 
over-all designs appear to be similar to 
that of the Li-Liu pickup and are be- 
lieved capable of operation at elevated 
temperatures. The Swiss Locomotive 
and Machine Works (in Switzerland) 
piezo-electric (SLM) pickup has a num- 
ber of novel design features. A sealed 
steel tube, not unlike the strain tube, 
tightly compressing two quartz pieces of 
semicircular cross section constitutes the 
sensing element. Axial loading of the 
cylinder results in both axial and trans- 
verse stresses in the quartz to produce 
the piezo-electric effect. Although they 
are highly promising transducers, the 
SLM pickups nevertheless suffer from 
appreciable hysteresis (about 2 per cent 
maximum) and substantial temperature 
drift. At very high neutron flux a 
considerable amount of error is in- 
troduced. 

In using the high output impedance 
type of sensing element, like the piezo- 
electric and capacitance types, the con- 
necting cable is an important factor in 
the dependability of the measurement. 
Perls has pointed out that large errors 
can be caused by the migration of elec- 
tric charges when a coaxial cable is sub- 
jected to vibration (17). His research 
led to the development of the low noise 
coaxial cable which, when used in a short 
length adjacent to the pickup, can elim- 


5 Newer ceramic-type piezo-electric ma- 
terials such as Piezite II, lead metanio- 
bate, and lead zirconate-lead titanate 
show promise as sensing elements for lower 
temperature operation. 


inate these vibration generated noises. 

The magnetic type of sensing ele- 
ments, such as the Ultradyne variable 
reluctance pickup, has good zero sta- 
bility and high signal-to-noise ratio. 
Many magnetic pickups, however, are 
sensitive to temperature effects and 
their responses are limited by the rela- 
tively low carrier frequency—usually 
less than 20 ke—and by the core loss. 
The inductive type of gage, which is less 
susceptible to temperature effects and 
has a relatively large output, can be 
made linear over a given range but still 
suffers from core losses at high carrier 
frequencies, and the relatively large 
moving mass prevents its use for meas- 
urements above 1000 cps. A few spe- 
cially built types, among which is one 
German pickup, are known to be able to 
operate up to 3 ke. 

Generally speaking, pickups which use 
changes in inductance and capacitance 
do not have wide linear range. By re- 
stricting their use to a small part of the 
general curve or through use of external 
linearizing devices, closer linear opera- 
tion can be obtained. 

Pressure gages using the optical inter- 
ference principle originally reported by 
Buck and Barkas in 1948 (18), have 
gained increasing acceptance in recent 
years although they are bulky and can- 
not be used at elevated temperatures. 
The data presentation form also is quite 
awkward to use. 


Diaphragmless Devices 


The possibility of sensing pressure 
without resorting to diaphragm-type 
transducers has been investigated re- 
cently by aeronautical researchers. 
Some of their investigations in the low 
pressure measurement field have in- 
cluded the well-known Alphatron, x-ray 
and electron beam methods (19). 
These are primarily density measuring 
devices, but for simple gas systems at 
constant temperature their outputs can 
be calibrated in terms of pressure. The 
Alphatron uses a radioactive source 
which emits alpha particles at a speed 
0.05 that of light. When the alpha 
particles collide with gas molecules, the 
latter are ionized by the loss of an elec- 
tron. Under standard atmospheric con- 


ditions, a single particle will ionize 10‘ 
air molecules per cm. Thus, an elec- 
tric field imposed between the electrodes 
generates a current proportional to the 
number of ionized particles produced per 
unit time. Over wide ranges of pres- 
sure, however, the output of the Alpha- 
tron is by no means linear. 

The electron beam method depends 
partly upon the absorption and partly 
upon the scattering of an electron beam 
by the molecules in gases (20), and the 
x-ray technique makes use of the absorp- 
tion of soft rays, i.e., x-rays with wave 
length longer than 2 A to measure the 
density (21). These two methods are 
based upon the familiar absorption law 


I = 


where yu is the absorption coefficient, and 
is a function of the density and the 
atomic number of the absorbing gas; | 
is the thickness of the layer of matter 
to be measured. Arnold, Dimeff and 
Kistiakewsky were among the workers 
who have applied the x-ray techniques 
to measurement of transients (22). 
Except in special applications, however, 
none of these systems has, at the present 
stage, attsined the degree of simplicity 
and versatility of the diaphragm-type 
sensing element systems. A knowledge 
of the status of development in dia- 
phragm design is, therefore, essential to 
dynamic pressure measurement engi- 
neering. A discussion of this topic is out- 
side the realm of the present paper. The 
classic writings in the field are well 
known and have been reviewed in (23). 
Analytical information on practical de- 
sign aspects of various pressure receiv- 
ing elements is presented in (24). 


Very High Speed Transient 
Pressure Measurement 


Electromechanical pickups cannot, in 
general, by themselves fulfill the re- 
quirements of ultra high speed measure- 
ment because: 

1 The ringing of pickups at their 
natural frequencies tends to obscure the 
waveform of the input pressure phe- 
nomena. 

2 They do not have the requisite 
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Fig. 3 (a) Response of transducer to 
shock wave using GN. (sweep speed: 
10 » sec/cm) 


Fig. 3 (b) Response of a Rutishauser 

transducer to a shock tube employing 

pneumatic feedback (sweep speed: 200u 
sec/cm) 


Fig. 3 (c) Response of Fig. 3 (a) equal- 
ized by DADEE (sweep speed: 10 u 


Fig. 3 (d) Response of Fig. 3 (b) ‘equal- 
ized by DADEE (sweep speed: 200. 
sec/cm) 


Fig. 4 Response to a step pressure using 
DADEE (sweep speed: 2 » sec/cm) 


Fig. 5 Driving function, simulated trans- 
ducer output, and equalized DADEE out- 
put (sweep speed: 10 u sec/cm) 


transient response and, therefore, can- 
not indicate the true pressure magnitude 
without overshoot or attenuation and 
phase shift. 

These inadequacies are immediately 
apparent when one examines Figs. 3(a) 
and 3(b) which show the response of 
two pickups respectively to a step func- 
tion and to a complex pressure function. 
The true waveforms of the same driving 
functions are shown in Fig. 3(c) and 
3(d). The driving functions shown in 
Fig. 3(c) and 3(d) were measured and 
recorded instantaneously by the DAD- 
KE system (Dynamic Analog Differen- 
tial Equation Equalizer) developed by 
the authors at Rocketdyne (25). This 
system and the Dyna-electronic Tran- 
sientgraph are two new methods by which 
fast pressure phenomena can be meas- 
ured without any postmeasurement 
correction for dynamic errors except an 
extremely small amount of time delay. 
With the present model of the DADEE, 
pressure transients with rise times as 
short as 0.5 microsec have been meas- 
ured. New models promise even faster 
response. Fig. 4 shows a 1-microsec 
pressure front taken with a Rutishauser 
pickup whose natural frequency is 80 
ke. 

Briefly, the DADEE system automat- 
ically compensates for the transducer 
dynamic characteristics during any 
transient or steady-state measurement. 

An exhaustive series of tests has been 
carried out at Rocketdyne in order to 
prove the validity of this system. To 
illustrate the results, Fig. 5 shows a driv- 
ing function in the form of an electronic 
pulse whose rise time is estimated to be 
of the order of millimicroseconds. The 
ordinary response of a simulated pickup 
is shown next to its driving function. 
Note that this response function has 
very little resemblance to the original 
function. On the contrary, the output 
of the DADEE system as shown in the 
bottom of the picture is a close repro- 
duction of the original function. 

With DADEE designed on the basis 
of a linear differential equation, the 
question naturally arises as to the ap- 
plicability of this technique to nonlinear 
pickups. Nonlinear pickups are known 
by such response characteristics as jump 
phenomena and ultra- and subharmonic 
resonances. In some well designed 
pickups these nonlinear characteristics 
are minimized through design skill but, 
when this is insufficient, we must resort 
to the technical solution of the nonlinear 
problem. An effective technical tool in 
approaching the nonlinear pickup 
problem is the Transientgraph system 
which was developed by Liu and is based 
upon the phase-plane method (Fig. 6). 


6 A detailed description of the DADEE 
system will be published separately. The 
systems have been submitted for patent 
application by Rocketdyne under the 
authors’ names. 


A x (normalized time derivation of 2) 


Fig. 6 Dyna-electronic Transientgraph 
(nat. frequency = 60,000 cps) 


As a tool for solving nonlinear 
systems, this approach has been studied 
by Andronow and Chaikin (26) and 
others. When modified to apply the 
Ritz-Klotter method, it is possible to 
extend it to the solution of steady-state 
nonlinear vibration problems (27). This 
presentation provides, in a single pic- 
ture, all the essential information in the 
time and space domain from which 
natural frequency, damping ratio, etc., 
can be determined. With this infor- 
mation available, one can use simple 
graphical methods to reconstruct ac- 
curately the true waveform from the 
phase-plane diagram’ The present 
models of the Transientgraph are prac- 
tical systems which can be used to deal 
with transducers having nonlinear char- 
acteristics and even those with more 
than one degree of freedom. 
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Increased activity in the 
design and production of sys- 
tem electronics has created 

openings for engineers in the 
following areas: 


ELECTRONIC AND AIR DATA 
SYSTEMS Required are men of proj- 
ect engineering capabilities. Also 
required are development and 
design engineers with specialized 
experience in servo-mechanisms, 
circuit and analog computer design 
utilizing vacuum tubes, transistors, 
and magnetic amplifiers. 


SERVO-MECHANISMS 

AND ELECTRO-MAGNETICS Complete 
working knowledge of electro-magnet- 
ic theory and familiarity with materi- 
als and methods employed in the design 
of magnetic amplifiers is required. 


FLIGHT INSTRUMENTS AND 
TRANSDUCER DEVELOPMENT 
Requires engineers capable of analyz- 
ing performance during preliminary 
design and able to prepare proposals 
and reports. 


FLIGHT INSTRUMENTS 
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low friction, freedom from vibration 
effects and compensation of thermo 
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Electrical input 10 milli-amperes 
Operating pressure__200 to 4000 psi 
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New AiResearch Time Dwell 
Servo Valve overcomes performance and reliability 


limitations of present hydraulic control systems 


Acceleration switching is a new contro] technique ¢ Less than 5% null shift between —20°F 
which provides for positive control of spool velocity, and +700°F 7 
allowing the servo system engineer to achieve resolu- _—* ‘Infinite resolution 7 
tion, reliability and response previously unattainable. — * High spool control forces ( 100 pounds) — 
Significant advantages include: at all signal levels . 
© Unimpaired operation in contaminated oil ¢ High pressure gain (10 times linear valve) 
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Astrophysical Quantities, by C. W. Allen, 


John de Graeff, Inc., New York, 1955, 
263 pp. 


Reviewed by Donatp H. Loucuripce 


Reviewed by 8S. F. Sincer 
University of Maryland 


This is probably the one indispensable 
handbook for an astrophysicist or anyone 
interested in the earth, solar system and 
stars. The unusual feature about the 
book is that it has been prepared by a 
single individual who is himself an out- 
standing astrophysicist and this leads to 
an important element of coherence which 
is so often lacking in handbooks. The 
volume is unusually complete in giving 
data connected with spectroscopy and 
atomic physics. Such matters as ionization 
potentials, electron affinities and atomic 
cross sections for collisions are well dis- 
cussed, as are spectroscopic quantities, 
such as oscillator strengths, line broadening 
and strength of multiplets. One whole 
chapter deals with radiation phenomena, 
others with the earth, the sun, planets and 
satellites, interplanetary matter, normal 
stars, stars with special characteristics, star 
populations, clusters and galaxies, and in- 
terstellar space. There are of course very 
useful tables of general constants and 
units, and special tables which simplify 
astronomical calculations, e.g., tables for 
conversion from equatorial to galactic 
coordinates. An especially valuable fea- 
ture of the book are the many good refer- 
ences. Every piece of information given 
is backed up by a reference to an original 
paper. 


Nuclear Chemical Engineering, by Man- 

son Benedict and Thomas H. Pigford, 
McGraw-Hill, New York, 1957, vii + 
594 pp. $9.50. 


General Motors Corp. 


q The need for a comprehensive text on 
nuclear chemical processing has been 
realized by those working in nuclear en- 
gineering for some time. The present 
volume most adequately fulfills such a re- 
quirement. Both authors have had ex- 
tensive experience in the fields covered by 
the present work, and in particular the 
senior author was one of the leaders during 
the war in the development of the gaseous 
diffusion method of uranium isotope sep- 
aration. 

The text is divided into twelve chapters, 
the first two of which provide a brief sur- 
vey of nuclear properties and those im- 
portant chemical elements entering into 
nuclear power applications. A clear pres- 
entation of the theory of radioactive de- 
cay chains and the neutron reactions en- 
tering into the fission process is presented. 
Next, the various fuel cycles appearing in 
thermal nuclear reactors and their quanti- 
tative considerations are fully developed. 
Illustrative problems are provided at the 
end of each chapter, thus greatly increas- 
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Ali Bulent Cambel, Northwestern University, Associate Editor 


ing the usefulness of the book as an ad- 
vanced textbook for senior or graduate 
use. 

Production processes utilized in pre- 
paring uranium feed material and the prop- 
erties of those important elements, zir- 
conium, thorium and beryllium, are pre- 
sented with the necessary thermodynamic 
data to make such sections useful as refer- 
ence material. 

Theory and practice of solvent extrac- 
tion of metals is here presented in the ful- 
lest detail (outside of AEC reports), and 
the properties of radiated fuel naturally 
lead the chemical engineer studying the 
text into the basic problem of the separa- 
tion of reactor products. As might be 
expected in a text from such authorities, 
the subject of isotope separations of both 
light and heavy elements, the various 
theories underlying such processes, and 
the practical results obtained by the 
various methods to date, along with the 
economics of the various processes, oc- 
cupy about the last third of the book. 

The Appendix contains a section on 
fundamental physical constants, tables of 
the separation potential, tables of mass 
and energy equivalents, and a long table 
giving the known properties of the nu- 
clides. 

It is to be hoped that soon each chemi- 
cal engineering department in leading 
universities will be adequately staffed to 
present the material covered by this book, 
and without doubt instructors in the field 
will feel a great debt of gratitude to the 
authors for their obviously tremendous 
task of preparing a text and reference 
book of such high standard. The print- 
ing is large and clear, the drawings ade- 
quate in size, and the mathematical typog- 
raphy well arranged. 


New Horizons in Astronomy (vol. 1, no. 
1, Smithsonian Contributions to Astro- 
physics), edited by Fred L. Whipple, 
Smithsonian Institution, Washington, 
D. C., 1956, 181 pp. $1.50. 

Reviewed by S. F. SiIncer 
University of Maryland 


This is a first issue of Smithsonian 
Contributions to Astrophysics, which is to 
consist of short monographs on astronom- 
ical topics. This first issue contains a 
large number of short contributions which 
give an impression of where astronomy 
stands in various fields. The main head- 
ings are: Techniques and Instrumenta- 
tion, which includes. optical and radio 
methods for studying astronomical ob- 
jects. The next heading. is Related 
Sciences, which includes a short discus- 
sion of rocketry and fields of physics 
which are of importan¢e to astrophysics. 
The third main section is the Solar Sys- 
tem and includes the earth and smaller 
planets and meteors. The final section 
deals with Stars and the Galaxy. As ex- 


pressed by the editor, Fred L. Whipple, 


who is the Director of the Smithsonian 
Astrophysical Laboratory, one of the main 
purposes of the volume is to encourage es 
tablished astronomers to pause for 
moment and reflect on the most advan 
tageous planning in their own field of re 
search, and supply them with a broade: 
insight into the activities in other fields 
of astronomy, possibly related to thei 
own. The first volume is perhaps less 
useful to a non-astronomer, since in som 
cases the discussion is very cursory and 
without references, in other cases too 
specialized. It is, however, a book which 
is easily read and should therefore be 
studied by anyone who has an interest 
in astrophysics. 

Subsequent issues are being written by 
individual authors who will cover many 
important astrophysical topics of current 
interest. 


Experimental Design: Theory and Ap- 
plication by Walter Federer, Mac- 
millan, New York, 1955, xix + 544 
pp. with 47 planographed problems 
bound with the book. $11. 

Reviewed by Frep M. WESTFIELD 
Northwestern University 


This is a text in statistics which will be 
used as a reference book in many labora- 
tories. It deals with schemes for plan- 
ning experimental investigations so as to 
permit efficient application of the tools of 
statistical inference. 

While the subject matter of this book 
has become especially important for those 
many areas of modern research where 
costs of observations are high and the 
extent of laboratory control decidedly 
limited, it is to be regretted that most of 
the specific examples and problems are 
taken from biological and agricultural ex- 
perimentation. But for this the author 
can hardly be blamed. Historically, this 
is largely where the techniques of experi- 
mental design were studied and developed, 
and the title page tells us that the author 
is professor of biological statistics in the 
Department of Plant Breeding at Cornell 
University. 

The book is readily accessible to the 
reader familiar with basic statistical tech- 
niques as taught in a good first-course in 
statistics at the college level. Knowledge 
of the calculus or matrix algebra is not 
required. The bibliography is excellent. 


Light Scattering by Small Particles, 
edited by H. C. van de Hulst, Wiley, 
New York, 1957, 470 pp. $12. 

Reviewed by S. F. SincER 
University of Maryland 


This book has already become a classic 
in its field and should be in the library of 
anyone who deals with scattering prob- 
lems. It is particularly concerned with 
the scattering of electromagnetic waves, 
light and radiowaves, and is by far the most 
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PRECISION ACTION, SMALL SIZE, LOW COST are features of 
Bristol’s new pressure switches for aircraft electrical circuits. 
Miniature size shown here. _ 


New! Miniature pressure 
switch for aircraft use 


Accurate, reliable, repeatable performance in any position 
and under MIL-spec environmental requirements is the design 
aim for Bristol’s ® new line of pressure switches. 

Designed specifically to meet aircraft requirements, the 
switches are precision devices for switching electrical circuits 
in response to pressure changes in gases and liquids. They’ re - 
available i in both regular and miniature sizes. 

Specially-designed capsular elements—Al stainless steel or Ni-Span 
C, welded construction. Exclusive design assures maximum re- 
sistance to vibration. 

Absolute, gage, and differential pressure models available. 
Hermetic sealing affords maximum 
protection. All exterior joints are 
metal-to-metal or metal-to-glass. 
Outstanding over-range protection — 
Built to take high over-pressures 
without damage or loss of calibra- 
tion. 

Variety of mounting arrangements 
possible— Clamps, studs, or pres- 
sure fittings can be utilized to 
mount switch in any appropriate 
location. BRISTOL'S REGULAR-SIZE pres- 

Write for bulletin AV2004. switch. Both 
Bristol Company, 175Bristol Road, *" 


made in absolute, gage, and 
Waterbury 20, Conn. 6.44 differential pressure models. 


TYPICAL SPECIFICATIONS 
Miniature Regular 
Size Size 
Pressure Setting: Any pressure be- Any pressure be- 
tween 2 psi and tween 5 psi and 
100 psi, absolute, 150 psi, absolute, 
gage, ordifferen- gage, or differen- 
tial, as specified tial, as specified 
Contact Ratings: 5S amp and 10 amp 5 amp and 10 amp 
resistive at resistive at 
30v d-c or 115v a-c d-c or 115v a-c 


Weight (ounces): 1% 
Size (inches): 1% long 3% long 
15/16 diam. 2% diam 
7 Ambient Range: —65°F to 250°F —65°F to 250°F 
ag Vibration & Shock: Designed to meet Designed to meet 
Spec. MILE-5272A Spec. MIL E-5272A 


FINE PRECISION INSTRUMENTS 
FOR OVER 67 YEARS a 


comprehensive treatise on this subject. 
Fortunately it has been written by a 
single individual who is an expert in all of 
the applications as well as the theory, and 
who has made many important contribu- 
tions himself, particularly in astronomy. 
Van de Hulst first covers basic scattering 
theory, applicable to any problem; most 
of the book, however, is devoted to scat- 
tering by special types of particles. 
There are simplifications which can be 
introduced, either when the particles are 
small compared to the wave length or 
when the particles are large compared to 
the wave length. The rigorous theory 
can be applied for spheres of arbitrary 
size and has been developed by Mie. 
There one distinguishes between non- 
absorbing spheres which introduce an 
important simplification, and spheres 
with refractive index near 1 which can be 
handled in a special way. In the case of 
very large spheres, other approaches can 
be applied. The case of absorbing spheres 
falls again in a different category. Van 
de Hulst discusses in detail the optics of 
the rain drop, a classical problem. Then 
there is also a discussion of scattering by 
circular cylinders and particles of other 
forms. The applications of scattering 
and extinction experiments are many. 
They are applied to chemistry to deter- 
mine the size of small particles and macro- 
molecules. Scattering experiments are of 
importance in aerosol physics. In meteor- 
ology scattering is of importance in con- 
nection with atmospheric haze and dust. 
Scattering phenomena occur in clouds and 
fog, and the subject is of great importance 
for radar meteorology. 

It is in astronomy, however, where 
scattering phenomena give us the greatest 
information about inaccessible objects. 
The author takes up applications to 
planetary atmospheres and to his particu- 
lar specialty, the determination of the 
size of interplanetary dust particles from 
the zodiacal light, which is really scat- 
tered sunlight. 

For the practical user the book is es- 
pecially valuable because it tells him just 
where and how to apply scattering for his 
measurements; the book summarizes 
numerical results and graphical data 
which are spread throughout the literature 
and are now made accessible for the first 
time in usable form. 7 


Book Notices 


Engineering Enrollment in the United 
States, edited by Norman N. Barish, New 
York University Press, 1957. This vol- 
ume presents statistics of enrollment in 
engineering education. It discusses and 
interprets the data. Future requirements 
for engineering education are suggested. 

Energy, by Sir Oliver Lodge, John F. 
Rider Publisher, Inc., 1957. This is a 
new version of Sir Oliver Lodge’s publica- 
tion first written in 1929. 

Principles of Engineering Heat Trans- 
fer, by Warren H. Giedt, D. Van Nostrand 
Company, Inc., 1957. This volume is ex- 
cellently suited as a textbook for the first 
course in engineering heat transfer. It 
covers conduction, convection and radi- 
ation. 
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AMF has missile experience you can use. 


© Auxiliary Power Supplies 


FEBRUARY 1958 


Looking for solutions? 


~e@ Building a reliable missile system introduces problems at every step... the | 
; kind AMF solves daily. From drawing board to target, AMF is constantly — 


bridging the gap between missile concept and performance. * This wealth — 


of experience, gained as a leading contributor to numerous major missile 


programs, can produce the solution to your particular problem. ¢ For a full — 
~ description of AMF engineering and production facilities in the missile field, 
as well as a review of their performance, contact the AMF Defense Products 


Manager in any of the cities listed below. 


@ Armament 
© Ballistics 
@ Radar Antennas 
Guided Missile 
Support Equipment 


# 
. 
Defense Products Group 


AMERICAN MACHINE & FOUNDRY COMPANY 


1101 North Royal Street, Alexandria, Va. 


DEFENSE 


PRODUCTS 
© Control Systems 


Asbury Park @ Atlanta @ Boston @ Brooklyn @ Dallas @ Dayton @ Los Angeles @ Seattie ® Tucson @ Washington, D.C. 
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ENGINEERS 


Aerodynamics & Propulsion 
Information 
manual about 
APL and its 


programs 


now available 


The Applied Physics Labora- 
tory (APL) of The Johns 
Hopkins University is unique 
in that we are neither an in- 
dustrial nor an academic 
organization, but rather a 
composite, having drawn 
freely from the methodolo- 
gies of each. 

For thirteen years APL has 
pioneered in guided missiles. 
Today we are engaged in a 
broad program of R & D for 
the Navy; in addition, we are 
responsible for technical di- 
rection of industrial and 
academic contractors in de- 
veloping the Terrier, Talos 
and other major weapons and 
weapons systems. Our staff 
members enjoy not only the 
stimulus of association with 
their immediate colleagues at 
APL, but also with those in 
other organizations of con- 
siderable stature. 


NEW 30-PAGE PUBLICATION 


A few positions for senior engi- 
neers and scientists are now 
open. Information on our ac- 
complishments and goals is 
available in a new 30-page pub- 
lication, just off the press. 

In it staff leaders representing 
each of the various disciplines 
and fields outline the nature of 
their programs. Information on 
our new laboratory in Howard 
County, Md. (equidistant be- 
tween Baltimore and Washing- 
ton) is also included, together 
with facts on the outstanding 
communities in which our staff 
members live. 


Quantity 1s somewbat limited. 
May we suggest you send now to: 
Professional Staff Appointments, 


The Johns Hopkins University 
Applied Physics Laboratory 


8617 Georgia Avenue, Silver Spring, Md. 


New Patents 


Jet action rotors for use in a jet propulsion 
system (2,812,898). Ernest H. Buell, 
Albany, N. Y. 

Turbo compressor system having a 

first hollow outer shaft and containing 
an inner shaft. An outer compressor 
rotor is formed in two sections, one in- 
tegral with the hollow shaft. A second 
hollow outer shaft is integral with the 
other of the outer compressor rotor sec- 
tions. This shaft receives an end of the 
inner shaft, forming a combined unit for 
rotation in the direction opposite to the 
outer shafts and rotors. 
Intake sprinkler for gas turbine engines 
(2,812,899). R. G. Meschino, 
Canada, assignor to A. V. Roe Canada, 
Ltd. 

A rotatable spinner fairing into the 
contour of the wall of the central fairing 
adjacent the intake. A passage for liquid 
extends through the spinner shaft and 


head opening into the intake through a 
hole in the surface of the head. 


Inclined engine installation for jet air- 
craft (2,812,912). Wm. P. Stevens and 
Arthur Nitikman, Inglewood, Calif., as- 
signors to Northrop Aircraft, Inc. 

Engine assembly extending forwardly 
and downwardly with the longitudinal 
axis at an acute angle to the fuselage axis. 
The opening of the air intake duct is 
adjacent and below the lower surface of 
the fuselage. 

Fuel injection system for ramjet aircraft 
(2,812,978). Louis 8. Lancaster, 
Calif., assignor to the U.S. Navy. 

Flush skin fastening of ‘fuel injection 

nozzles. A cylindrical sleeve is receivable 
into the inner body so that its retaining 
collar is even in surface with the diffuser 
inner body surface. 
Jet deffecting device (2,812,980). M. 
Kadosch, F. G. Paris, J. Bertin and R. H. 
Marchal, Paris, France, assignors to 
SNECMA Co. 

Laterally extending guide vanes posi- 
tioned downstream of the nozzle outside 
the normal flow path of the jet. Vanes 
define successive passages having inlet 
ends spaced along the flow path, and out- 
let ends opening into the atmosphere at 
a substantial angle with the nozzle axis. 
Sight mount for a rocket launcher (2,814,- 
118). P. I. Evans and H. Harvey, 
Redondo Beach, Calif., assignors to the 
U.S. Army. 

Range plate fixed on a bracket coaxially 
of the elevation axis. The plate is op- 
erable to hold a sight mount body in a 
selected elevation adjustment about the 


George F. McLaughlin, Contributor 


axis of a post journaled in a normally 
horizontal pivot bearing. 

Return burning motor (2,814,179). L. 13. 
Edelman and C. F. Miller, Sacramento, 
Calif., assignors to the U. 8. Navy. 

Rocket motor comprising a front plate 

from which a first cylindrical wall ex- 
tends forwardly from the rear of the 
rocket, terminating to the rear of the front 
plate, forming an annular chamber and a 
cylindrical chamber inwardly. A singie 
propellant charge fills the chambers, 
burning from the front of the annular 
chamber rearwardly, and from the rear of 
the cylindrical chamber forwardly. 
Flow responsive dump valve for gas tur- 
bine fuel systems (2,814,180). J. 1’. 
Hession, Indianapolis, Ind., assignor to 
General Motors Corp. 

Bypass having an inlet and an outlet 
with a fluid conduit on opposite sides of « 
pressure drop creating means. A con- 
trol connected to a servo actuator is pro- 
vided for initiating and discontinuing 
fluid flow. 

Aircraft propulsion apparatus (2,814,349). 
D. W. Berry, Kansas City, Kans., as- 
signor to Westinghouse Electric Corp. 

Jet propelled helicopter rotor blades 
with ducts extending from the hub to ex- 
haust nozzles at blade tips. A gas tur- 
bine is disposed in each duct. A driving 
means, interposed between the compressor 
and the turbine power take-off shaft, de- 
livers compressed air to the turbines. 
Aircraft duct screen arrangement (2,814,- 
454). F. W. Atkins and H. W. Elkin, 
Pasadena, Calif., assignors to North 
American Aviation. 

Foldable screen covering the power- 

plant air intake duct and carried by the 
retractable landing gear. The screen 
moves from a stowed position within a 
well, to an operative position in which it 
screens air entering the intake duct. 
Flap type variable area thrust nozzle 
(2,813,395). Robert E. Meyer, Glaston- 
bury, Conn., assignor to United Aircraft 
Corp. 

Jet engine thrust nozzle surrounded by 
an elongated ring with flaps pivoted at 
the discharge end. Openings in the for- 
ward end of the ring effect an injector 
action at the flaps to draw cooling air 
through the nozzle. 

Tail cone and resilient mounting therefor 
(2,813,396). Ralph Kress, La Mesa, 
Calif., assignor to Solar Aircraft Co. 

Packing of braided metallic gaskets to 

dampen vibrations tending to pass from 
the inner tail cone to the outer shell dur- 
ing operation of a jet engine. 
Thermal expansion means for combustion 
chambers (2,813,397). Earl R. Fisher, 
Charles D. Robin and Richard C. Hurd, 
Manchester, Conn., assignors to United 
Aircraft Corp. 

Fuel manifold with an air swirler unit 
comprising an inner sleeve telescoped in 
sliding fit into the fuel nozzle exterior. 
An outer sleeve is spaced externally from 
the inner sleeve, with vanes projecting 
between. Outer and inner cups at the 
forward end of the combustion chamber 
are free to expand without stress to the 
parts. 


Epitors Norte: 


Patent Office. Printed copies of Patents m 


Patents listed above were selected from the Official Gazette of the U.S. 
may be obtoined from the Commissioner of 


Patents, Washington 25, D. C., at a cost of 25 cents each; design patents, 10 cents. 
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push-button 
warfare 


The RCA Talos Defense System is the first completely 
automatic land-based system for launching and guiding 
missiles, and utilizes the Talos Missile developed by 
the Applied Physics Laboratory and produced by 
Bendix. The Defense Unit receives target signals from 
remote outposts, analyzes them with regard to number 
of attackers, location, course and speed. Next, com- 
puters determine the logical points of interception, 
order the missiles loaded on launchers, guide them at 
supersonic speei to the vicinity of the target, after 


RADIO CORPORATION of AMERICA 
DEFENSE ELECTRONIC PRODUCTS 
CAMDEN, N. J. 


This Talos launching from the RCA Defense System on Dec. 13, 1957, resulted in direct hit on distant drone plane 


which the missiles ‘‘lock”’ on the target and close in 
for its destruction. All without even the touch of a 
button! The RCA Talos Defense System, with its 
electronic equipment and guidance systems, was 
designed, developed and built by RCA as prime 
contractor, aided by many subcontractors. It was 
turned over to the U. S. Army on October 15, 1957, 
and is a missile milestone, exemplifying the continuing 
determination of American enterprise to secure peace 
with honor and justice. 
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Wind tunnel pressure pick-ups (2,813,- 
421). Gordon H. Cheney, John H. Rus- 
sell and Phillip C. W hitener, Seattle, 
Wash., assignors to Boeing Airplane Co. 
Wall with apertures aligned at inter- 
vals lengthwise of the airflow. A carriage 
guided along the line of apertures has 2 
part which registers with successive 
eo apertures during its advancement. 
A stylus, electrically operated from the 
carriage, records pressure as the part 
comes in registry with each aperture. 


Shaped propellant charges for solid fuel 
rocket type motors (2,816,418). A. C. 
Loedding, Princeton, N. J., assignor to 
Unexcelled Chemical Corp. 

An annularly disposed charge having an 
axial opening of hill and valley shape for 
regulating the burning characteristics of 
the charge radially outward. The charge 
consists of four quadrants with triangular 
stringers between. 

Hydronamic bullet catcher (2,813,422). 
Karl E. Schuessler, Philadelphia, Pa., 
assignor to the U.S. Navy. 

Trap for high speed projectiles consist- 
ing of a tube with an inlet for projectiles 
and water. A solid conical surface of 


water adjacent to the inlet faces and im- 
pedes projectiles entering the tube. 
Bali and socket joint with internal bearing 
means (2,813,731). W. D. Tracy, A. L. 
Mager and C. A. Dasney, Arcadia, Calif., 
assignors to Southwest Products Co. 

Duct assembly with a gasket ring re- 
cessed in an outer tubular member having 
an annular concave bearing surface. The 

gasket abuts a convex surface on an inner 
‘linine member. Members are held to- 
gether at a hub supported on spiders, per- 
mitting angular movement of the inner 
ring. 
Leak recovering device for pulse jet units 
(2,814,930). H. L. P. Meulien, M. J. 
Baraut and J. Le Foll, Pre Saint-Gervais, 
France, assignors to SNECMA Co. 

Jet propulsion plant in which combus- 
tion gases flow both rearwardly through 
the exhaust and forwardly through the in- 
take portion, to issue from both ends of 
the pulsejet units. Tubular members 
extending forward are spaced to allow 
free passage of ambient air into the intake, 
Gas-turbine power plant installations 
with means for disposal of drainage fuel 
(2,814,931). C. L. Johnson, Allestree, 
England, assignor to Rolls-Royce, Ltd. 


Fuel drain arrangement including a 
venturi passage having an inlet con- 
nected to a source of air pressure. Fuel 


in a collector tank is entrained in the air 
flowing through the venturi passage and 
delivered into the exhaust gas stream, 
carrying the fuel clear of the installation. 
Fuel system for gas-turbine engine reheat 
combustion stage (2,814,928). D. O. 
Davies and A. Jubb, Bradford, England, 
assignors to Rolls-Royce, Ltd. 

Device located between the compressor 
delivery and the turbine inlet, and be- 


HIGH PRECISION 


MACHINE PARTS 
SINCE 1908 


The equipment you are designing or building 
may well include a number of small high pre- 
cision machine parts, and you will be anxious to 
see that they are produced by a firm in which 
you can have complete confidence. A firm com- 


\ chure will be mailed promptly upon request. 


petently staffed, well equipped, with years of 
experience, and a reputation for high quality pro- 
duction and complete dependability. LaVezzi, 
with a history of many years in this specialized 
field is the answer to this need, and earnestly 
solicits your quote requests. An illustrated bro- 
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tween the turbine and the propelling noz- 
zie. The device is sensitive to a_pre- 
selected value of the ratios of pressures at 
the various points, and operates the 
throttle to maintain a preselected value of 
the ratio. 


exit door 
Pacifie Pali- 


actuated missile 
R. T. Johnson, 
Calif., assignor to Northrop Air- 
craft, Ine. 

Conduit in an airplane wingtip pod 


Exhaust 
(2,816,483). 


sades, 


rearwardly of a rocket launching tube 
providing a passageway for exhaust gases. 
An actuating flap in the conduit is moved 
from closed to open position upon im- 
pingement of exhaust gas on the flap. 

Fuel supply control for rocket-type jet 
propulsion units (2,814,929). A.W. Mor- 
ley and A. R. Mortimer, Ickenham, 
England, assignors to D. Napier & Son, 
Ltd. 

Bullet-like member extending into and 
axially movable within the nozzle passage 
of a reaction chamber to vary its cross- 
sectional area. As liquid propellant con- 
trol valves are opened, the member is 
moved progressively to increase the are: 
of the passage, and the supply working 
fluid is increased correspondingly. 

Ejection seat thruster (2,815,008). (. 
E. Hirt, New Haven, Conn., assignor to 
Taleo Engineering Co., Inc. 

Firing of a cartridge within a housing 
actuates a piston to move hydraulic fluid 
in a cylinder. An orifice plate restricts 
the flow of fluid from the cylinder to an 
adjacent buffer casing. 

Fuel (2,815,270). E. A. Lawrence and 
F. L. von Brecht, Monrovia, Calif., as- 
signors to Aerojet-General Corp. 

Fuel comprising 20 to 40 per cent (by 
volume) of nitropropane dissolved in 
nitromethane, and containing chromium 
acetony] acetonate. 


Fuel (2,815,271). F. Zwicky, F. J 
Ewing, J. M. Carter and A. J. Stosick, 
Pasadena, Calif., assignors to Aerojet- 


General Corp. 

Self-combustible, shock-proof stable fuel 

comprising a mixture of nitromethane 
with 25 per cent to 70 per cent (by volume) 
of the mixture of nitroethane and chrome 
acetony] acetonate. 
Vertically airborne aircarft design (181,- 
599). N. EK. Nelson and J. B. Reichert, 
Los Angeles, Calif., assignors to Doak 
Aircraft Co. 

Ducted fan type of aircraft supported 
in vertical position on the ground by 
shock-absorbing casters at the tips of di- 
rectional fins. 

Remote control fuze removing device 


(2,815,565). E. M. Johnson and R. 
eg Vallejo, Calif., assignors to the 
U.S. Navy. 


A vise which may be raised and lowered 
mounted beneath a drive head and fuse 
puller for holding ammunition to be de- 
fuzed. Clutch fingers engage when the 
drive head is rotated, unscrewing the fuze 
from the projectile. 
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Variable area nozzle (2,815,643). F. L. | 
Geary and W. G. Taylor Jr., Glaston- | 
bury, Conn., assignors to United Aircraft 
Corp. 

Flaps pivotally mounted to vary the 
area of a duct outlet. Pivots are located 
inwardly from their ends. Sealing means 
extend around the duct at the forward 
ends of flaps, forming a chamber in com- 
munication with the space between the 
duct and flaps. 

Diffuser (2,815,770). A. G. Thorp II, 
Wallingford, Pa., assignor to Westing- 
house Electric Corp. 

Preformed section for use between the 
compressor and combustion sections of an 
axial flow gas turbine. It comprises a 
pair of converging and a pair of diverging 
concentric tubular members spaced to 
form inner and outer annular fluid conduit 
sections. 

Control system for dual rotor pump drives 
(2,816,417). D. J. Bloomberg, Newton, 
Mass., assignor to General Electric Co. 

Gas turbine arrangement of separate 
means for supplying fuel to part of the are 
of one tier of blading on each rotor. 
Pump generated pressures govern opera- 
tion of the fuel supply means. 

Propellant displacement gas generators 
(2,816,419). H. E. Mueller, Kenmore, 
N.Y., assignor to Bell Aircraft Corp. 

Container of ammonium carbonate 

compound adapted to progressively de- 
compose when heated, evolving large 
quantities of gas to the oxidant and fuel 
tanks. 
High speed motor vehicle (2,816,615). 
L. Kraus, Stuttgart-Rotenberg, Germany, 
assignor to Daimler-Benz Aktiengesell- 
schaft. 

Jet producing device including a rotary 
Llower with its axis at an angle to the ver- 
tical central plane of the vehicle. When 
driving through a curve, a jet of air is dis- 
charged to counteract the effect of centrif- 
ugal force. 


2,816,720 


Incremental remote radio control sys- 
tem (2,816,720). W. G. Hodson, W. E. 
Peterson, H. H. Kenkel, J. J. Hagopian 
and S. R. Hennies, Englewood, Calif., as- 
signors to Northrop Aircraft, Inc. 

System in which the control signals are 
in increments, or pulses, rather than be- 
ing continuously proportional. Zero-drift 
is eliminated, rendering the attitude of 
the controlled vehicle exactly in corre- 
spondence with the indications of the con- 
troller. 

Rocket powered aerial vehicle (2,816,721). 
R. J. Taylor, Buffalo, N.Y., assignor to the 
U.S. Navy. 

Spaced forward and rearward combustion 
chambers located symmetrical with the 
vehicle. Propellant burns rearwardly in 
the forward chamber and forwardly in the 
rearward chamber. Thus, as the propel- 
lant in both chambers is consumed at an 
equal rate, the center of gravity of the 
rocket remains constant. 
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Performance Data 
Recorded on-Board 


| 


Economically, Faithfully 
with 
Model 409D 
Recording Oscillograph 


Accurate, reliable on-board collection of 
missile performance and environmental 
data is accomplished with the compact, 
lightweight Century 409D. Near 100° 
record recovery from hundreds of flights 
plus low per-channel cost gives a high 
data/dollar ratio. 


@ 12 trace recording on 50’ or 100’ record 
@ 5” x 5” x 11%” size; 12 lbs. (50’ magazine) 
@ Wide range of recording speeds 


@ Hi-contrast, permanent record 


Century Electronics & Instruments, Inc. 


Ask for Bulletin CGC 312 


1333 No. Ut T klahoma 
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Technical Literature Digest___ 


H. Smith, Associate Editor, and M. H. Fisher, Contributor | 
The James Forrestal Research Center, Princeton University | 


Combustion, Fuels and 
Propellants 


Hydrogen Peroxide and the Peroxide 
Compounds, by M. E. Pozina, Moscow, 
1951, 475 pp. (Mass. Inst. Tech., Dept. 

_ Chem. Engng., Rep. 52, April 1957). 

Phase Relationships in the Three Com- 
ponent Liquid System CO.-H,O-C,;H; OH 
at High Pressures, by Louis C. W. Baker 
and Thomas F. Anderson, J. Amer. Chem. 
Soc., vol. 79, May 5, 1957, pp. 2071-2074. 

Electronic Spectra of Cyclic Aromatic 


Hydrocarbon Radicals and Ions by H. C. 


Longuet-Higgins and K. Lenore McEwen, 
J. Chem. Physics, vol. 26, April 1957, pp. 

719-723. 

Gaseous Detonations in Mixtures of 
Cyanogen and Oxygen, by H. Milton Peek 
and R. G. Thrap, J Chem. Physics, vol. 
26, April 1957, pp. 704-745. 

Infrared Spectrum of Solid Nitrous 
Oxide, by David A. Dows, J. Chem. 
Physics, vol. 26, April 1957, pp. 745-747. 

Statistical Mechanics of Liquid-Vapor 
Equilibrium, by B. Widom, J. Chem. 
Physics, vol. 26, April 1957, pp. 887-893. 

Thermodynamic Formulas for Two- 
Phase Systems, by James C. M. Li, J. 
Chem. Physics, vol. 26, April 1957, pp. 
909-913 

Research on the Mechanism of Oxyda- 
tion of Single Iron Crystals at High Tem- 
perature and Under Low Oxygen Pressure, 
by Jean Bardolle, France, Ministere de 
V Air, Pub. Sci. Tech. 327, 1957, 62 pp. 
(in French). 

Titanium Cut J54 Weight by 16 Per 
Cent, by D. C. Goldberg, Aviation Age, 
vol. 27, June 1957, pp. 112-115. 

Researches on Corrosion and Inhibition. 
Adsorption, Inhibition and the Langmuir 
Equation, by George 8S. Gardner, J. Frank- 
lin Inst., vol. 263, June 1957, pp. 523-535. 

On the Activation Energy of High Tem- 
perature Creep in Metals, by Paul Felt- 
ham, Philos, Mag., 8th Series, vol. 2, May 
1957, pp. 584-588. 

Rocket Propellants for Underwater 
Missiles, by Alfred J. Zaehringer, Missiles 
and Rockets, vol. 2, May 1957, pp. 91-93. 

Reaction Bomb for Use at Elevated 
Temperatures and Pressures, by Edwin 
M. Larsen and James J. Leddy, Rev. Sci. 
Instrum., vol. 28, April 1957, pp. 292-293 

A Theory of Inflammability Limits and 
Flame-Quenching, by D. B. Spalding, 
Proc. Royal Soc., vol. A240, April 24, 
1957, pp. 83-89. 


Instrumentation, Teleme- 
_ tering, Data Recording 
; Data Reduction Center Clears Way for 
Missile Research, Jndustrial Labs., vol. 8, 
April 1957, pp. 48-49. 

Accuracy of the Cutler-Hammer Re- 
cording Gas Calorimeter when Used with 


Epitor’s Nore: Contributions from E. R. 
_G. Eckert, J. P. Hartnett, T. F. Irvine Jr. and 
P. J. Schneider of the Heat Transfer Labora- 
tory, University of Minnesota, are gratefully 
acknowledged. 
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Gases of High Heating Value, by John H. 
EKiseman and Elwin A. Potter, J. Res. 
Nat. Bur. Standards, vol. 58, April 1957, 
pp. 213-226. 


Instrumentation Used in Measurement 
of the Three-Dimensional Flow in an 
Axial Flow Compressor, by J. Horlock, 
Gt. Brit., Aeron. Res. Council, Curr. Pap. 
321 (ARC TR 17518), 1957, 12 pp. 

NAMTC Instrumentation Status Report 
as of 1 March 1950 with Addendum to 1 
March 1951, Naval Air Missile Test 
Center, Tech. Mem. Rep. 48, May 1951, 
106 pp. (Declassified from Confidential 
by authority of ASTIA TAB 123, p. 40, 
4/15/57.) 

Preset Gating Unit for Aeroballistic 
Testing, by Samuel E. Dorsey, Electronics, 
vol. 30, May 1957, pp. 164-169. 

Missile Telemeter Uses Transistor 
Amplifier, by John H. Porter, Electronics, 
vol. 30, May 1957, pp. 170-175. 

Ultra-High-Speed Flash Chemicro- 
graph, by J.S. Courtney-Oratt and Charles 
M. Huggins, Rev. Sci. Instrum., vol. 28, 
April 1957, pp. 254-255. 

Construction and Performance of a 
High-Speed Cinemicrograph, by J. 5S. 
Courtney-Pratt and Charles M. Higgins, 
Rev. Sci. Instrum., vol. 28, April 1957, pp. 
256-262. 

Metal System for Chemical Reactions 
and for Studying Properties of Gases and 
Liquids, by H. V. Neher and Satya Prak- 
ash, Rev. Sci. Instrum., vol. 28, April 1957, 
pp. 267-270. 

Device for Measuring Friction, by 
James L. Lauer and Patrick J. Friel, Rev. 
Sci. [nstrum., vol. 28, April 1957, pp. 294- 
295. 

Mercury in Quartz Thermometers for 
Very High Accuracy, by H. Moreau, J. A. 
Hall and Vera M. Leaver, J. Sci. Instru- 
ments, vol. 34, April 1957, pp. 147-151. 

Telemetering Ground Equipment for 
the Corporal E, by N. R. Silk, Calif. Inst. 
Tech., Jet Propulsion Lab., Progress Rep. 
4-32, July 1947, 31 pp. (Declassified 
from Confidential by authority of 
Calif. Inst. Tech., Jet Propulsion Lab., Pub. 
38, p. 7, 1/3/55.) 

Telemetering Nose Equipment for the 
WAC Corporal, by C. I. Cummings, 
Calif. Inst. Tech., Jet Propulsion Lab., 
Progress Rep. 4-33, July 1947, 42. pp. 
(Declassified from Confidential by author- 
ity of Calif. Inst. Tech., Jet Propulsion 
Lab., Pub. 38, p. 7, 1/3/55.) 

Recent Modifications of Corporal E 
Telemetering Systems, by F. W. Lehan, 
Calif. Inst. Tech., Jet Propulsion Lab., 
Progress Rep. 4-60, April 1948, 28 pp. 
(Declassified from Confidential by author- 
ity of Calif. Inst. Tech., Jet Propulsion 
Lab., Pub. 38, p. 8, 1/3/55.) 

Corporal E Telemetering Antenna 
Pattern Based on Model Tests at 400 
Megacycles, by Henry Curtis, Calif. Inst. 
Tech., Jet Propulsion Lab., Mem. 4-30, 
June 1948, 8 pp. (Declassified from 
Confidential by authority of Calif. Inst. 
Tech., Jet Propulsion Lab., Pub. 38, p. 10, 
1/3/55.) 

Air Borne Telemetering Equipment for 
the Corporal, by Alvin W. Newberry, 
Calif. Inst. Tech., Jet Propulsion Lab., 
Rep. 20-63, Jan. 1952, 38 pp. (Declassi- 


est Engineering 
at Marquardt | 


Roy E. Marquardt 
President 


Unique among air-breathing engines, 
the ramjet cannot run independently 
on the ground. Because ramjets depend 
on their velocity through surrounding 
air masses for compression of combus- 
tion air, intricate test facilities must be 
provided. 

Marquardt Jet Laboratory simulates 
these high-speed, high-altitude flight 
conditions during ground testing. A 
USAF facility, combined with special 
USN test facilities, MJL is one of the 
most extensive in existence today. 
Valued at $18 million and occupying 8 
acres of land, this dynamically func- 
tional engineering tool minimizes the 
number of costly trial and error flights 
of new engines. 

Complex test cells, with associated 
instrumentation and computers, per- 
mit sea-level and altitude testing of 
full size and test-scale models. In our 
two sea-level test cells, high pressure 
air is directed over the test engine 
through subsonic, sonic, and supersonic 
nozzles, and exhausted to the atmos- 
phere. Turbojets are ingeniously uti- 
lized in two full scale engine altitude 
test chambers to simulate heated and 
rarefied air conditions encountered by 
high-flying supersonic ramjets. 

Special configurations in both sea- 
level and altitude cells permit angular 
motion of the free jet nozzles. This 
allows evaluation of the effects of vary- 
ing angles of attack likely to be en- 
countered by operational engines. 

A new production acceptance test 
facility now being built in conjunction 
with our Ogden production plant will 
greatly supplement the up-to-the- 
minute equipment furnished the Mar- 
quardt Test Team. 

Within this Division, engineering 

openings exist for: 
MECHANICAL ENGINEERS AEROTHERMO ENGINEERS 
ELECTRICAL ENGINEERS CONSTRUCTION ENGINEERS 
For information about these positions 
and the professional engineering en- 
vironment at Marquardt, we invite you 
to write Jim Dale, Professional Per- 
sonnel, today. 


Marquardt 
; VAN NUYS, CALIFORNIA OGOEN, UTAH 


| FIRST IN RAMJETS 
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Marquardt Means Opportunity — Test 
engineers no longer need feel grounded by the lack of 
stimulating projects. At Marquardt Aircraft, the com- 
pany where an ENGINEER/BARRIER* has never existed, 
your work will span a broad range of supersonic 
propulsion problems. Look to your future by looking to 
Marquardt, today. Address your inquiries to Jim Dale, 


Professional Personnel, 
16556 Saticoy Street, marquardt 


Van Nuys, California. 
VAN NUYS, CALIFORNIA ~*~ OGDEN, UTAH 


Pictured above: Leigh Dunn, Chief Engineer, Test Sub-Division 
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from Confidential by authority of 
Gay Inst. Tech., Jet Propulsion Lab., 
2 Pub. 32, p. 15, 1/3/55.) 


s,s Theory of Free Molecule, Orifice Type, 
Pressure Probes in Isentropic Flows, 
Th G. N. Patterson, Toronto University, 
Tnot Aerophysics, Rep. 41, Nov. 1956, 
pp. 

Data Mechanization Speeds Ramjet 
Engine Testing, /ndustrial Labs., vol. 9, 
1957, p. 95. 

A New Adventure into Micro-Minia- 
turization. Instruments and Automation, 
vol, 30, April 1957, pp. 667-669. 


Programming a Computer, by Martin 
, L. Klein, Frank K. Williams and Harry C. 
Morgan, Instruments and Automation, 
vol. 30, April 1957, pp. 671-676. 

Recording 10 TC/SEC, by E. N. Kauf- 
man, Instruments and A utomation, vol. 30, 
April 1957, p. 703. 

Rapid Measurement of Thermal Con- 
ductivity by Transient Heating of a Fine 
Thermo-Junction, by R. A. W. Hill, 
Proc. Royal Soc., vol. A239, April 9, 
1957, pp. 476-486. 

High Speed Systems of Wind Tunnel 
Data Handling, by J. Lukasiewicz, J. A. 


Van der Bliek, and J. G. Scott, North 
Atlantic Treaty Organization, Advisory 
Group Aeron. Res. and Dev., Rep. 17, 


Feb. 1956, 24 pp. 

Absolute Intensities and Line Width 
Measurements, by D. Weber, Calif. Inst. 
Tech., Daniel and Florence Guggenheim 
Jet Propulsion Center, Tech. Rep. 23, 
March 1957, 6 pp. 

Cytac: Air Force’s New Position- 
Fixing System, by Henry P. Steier, 
American Aviation, vol. 20, April 8, 1957, 
p. 52. 

The Application of Rocket Sled Tech- 
niques to Flutter Testing, by W. R. Laid- 
law and V. L. Beals Jr., Inst. Aeron. Sci., 
Prepr. 666, Jan. 1957, 38 pp. 

Technique for Flutter Tests Using 
Ground Launched Rockets, with Results 
for Unswept Wings, by W. G. Molyneux, 
F. Ruddlesden and P. J. Cutt, Gt. Brit. 
Aeron. Res. Council, R&M 2944 (Formerly 
ARC TR 13562 and 14606, Roy. Aircr. 
Estab., Rep. Structures 72 and 118), 1956, 
19 pp. 


Materials of Construction 


Properties and Applications of Ferrites 
in the U. H. F. Faraday Effect, by M. 
Vassilev, Fusées et Recherche Aéronautique, 
vol. 1, no. 2, Oct. 1956, pp. 143-149 (in 
French). 


Aerodynamics of Jet 
Propelled Vehicles 


Analysis of Supersonic Wind Tunnel 
Tests on a 7g Scale Model of the Shore to 
Ship Guided Missile, by J. B. Kendrick, 
Calif. Inst. Tech., Jet Propulsion Lab., 
Mem. 4-15, Aug. 1947, 31 pp. (Dee lassi- 
fied from Confidential by authority of 
Calif. Inst. Tech., Jet Propulsion Lab., 
Pub. 38, p. 10, 1/3/55.) 

Characteristic Methods for Quasilinear 
Hyperbolic Differential Equations and an 
Application to Axially Symmetric Super- 
sonic Flow Past an Ogive, by H. K. Fors- 
ter, Calif. Inst. Tech., Jet Propulsion Lab., 
Rep. 4-34, Aug. 1947, 31 pp. (Declassi- 
fied from Confidential by authority of 
Calif. Inst. Tech., Jet Propulsion Lab., 
Pub. 38, p. 6, 1/3/55 .) 

Some Stability Considerations for a 
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Rocket Launched at Great Altitudes, by 
H. K. Forster, Calif. Inst. Tech., Jet Pro- 
pulsion Lab., Rep. 4-32, March 1947, 28 
pp. (Declassified from Confidential by 
authority of Calif. Inst. Tech., Jet Propul- 
sion Lab., Pub. 38, p. 6, 1/3/55.) 

Experimental Investigations with Jet 
Control Vanes, by W. B. Powell, Calif. 
Inst. Tech., Jet Propulsion Lab., Rep. 4-30, 
July 1948, 20 pp. (Declassified from 
Confidential by authority of Calif. Inst. 
Tech., Jet Propulsion Lab., Pub. 38, p. 6, 
1/3/55.) 

A Comparison of Long Range Surface to 
Surface Rocket and Ramjet Missile, by 
R. W. Krueger and J. E. Lipp, Rand Corp., 
R-174, May 1950, 37 pp. (Declassified 
from Secret by authority of AST/A TAB 
U-119, p. 46, 2/18/57.) 

The Erosion of Meteors and High-Speed 
Vehicles in the Upper Atmosphere, by C. 
Frederick Hansen, NACA TN _ 3962, 
March 1957, 38 pp. 


Heat Transfer and 
Fluid Flow 


Study of Chemical Relaxation Behind 
Shock Waves, by M. Steinberg and T. F. 
Lyon, Wright Air Dev. Center., Tech. Rep. 
56-394 (ASTIA AD 118050), Jan. 1957, 
18 pp. 

The Design of Supersonic Nozzles for 
Rockets, by R. P. Fraser and P. N. Rowe, 
Gt. Brit., Directorate Guided Weapons 
Res. and Dev., D.G.G.W. Rep. /EMR/ 
55/3 (Imperial College Sci. Tech., Rep. 
JRL 28), Oct. 1954, 22 pp., 33 figs. 

On Jet Separation in Supersonic Rocket 
Nozzles, I. The Characteristics of Flow, 
by P. Eisenklam and D. Wilkie, Gt. Brit., 
Directorate Guided Weapons Res. and Dev. 
D.G.G.W. Rep./EMR/55/4 (Imperial 
College Sci. Tech., Rep. JRL 29), May 
1955, 55 pp., 29 figs. 

Matrix Analysis of Heat Transfer Prob- 
lems, by Louis A. Pipes, J. Franklin 
Inst., vol. 263, March 1957, pp. 195-206. 

Experimental Study of Stagnation Tem- 
perature in Free Molecular Flow, by Mar- 
cel Devienne, Fusées et Recherche Aéro- 
nautique, vol. 1, no. 3, Dec. 1956, pp. 261- 
267 (in Frene h). 

The Numerical Integration of the Lam- 
inar Compressibility Boundary Layer 
Equations with Special Reference to the 
Position of Separation When the Wall Is 
Cooled, by G. E. Gadd, Gt. Brit., Aeron. 
Res. Council, Curr. Pap. 312, (formerly 
ARC TR 15101), 1957, 17 pp. 

Water Cooling Combustion Gases for 
Cubicle 3 Exhaust Chamber, by S. E. 
Matthews and C. R. Walton, Johns Hop- 
kins Univ., Appl. Phys Lab., CF 2619, 
Feb. 1957, 6 pp. 

Analysis of Heat-Driven Oscillations of 
Gas Flows, I. General Considerations, 
by H. J. Merk, Appl. Sci. Res., vol. 6, 
Sect. A, no. 4, 1957, pp. 317-336. 

Stability of Laminar Flow in Curved 
Channels, by Chia-Shun Yih and W. M. 
Sangster, Phil. Mag., vol. 2, 8th Ser. 
March 1957, pp. 305-310. 

Fluctuations and Irreversible Thermo- 
dynamics, by Laszlo Tisza and Irwin Man- 
ning, Phys. Rev., vol. 105, March 15, 1957, 
pp. 1695-1705. 

Spectrum of Turbulent Fluctuations 
Produced by Convective Mixing of Gra- 
dients, by Albert D. Wheelon, Phys. Rev., 
vol. 105, March 15, 1957, pp. 1706-1716. 


Intensity, Scale and Spectra of Turbu- 
lence in Mixing Region of Free Subsonic 


Jet, by James C. Laurence, NACA Rep. 
1292, 1956, 27 pp. (supersedes NACA 7'N 
3561 and 3576). 

Survey of the Acoustic Near Field of 
Three Nozzles at a Pressure Ratio of 30, 
by Harold R. Mull and John C. Erickson 
Jr., NACA TN 3978, Apr. 1957, 32 pp. 

The Cooling of Walls by Parietal In- 
jection, by Jean Berger (Aerothermal Se- 
minar of the Faculté des Sciences de Paris, 
Année 1954-1955), France, Ministére de 
VAir, Pub. Sci. Tech., Notes Tech. 63, 
1957, 110 pp., pp. 1-11 (in French). 

Heat Transfer from Liquid Metals 
Flowing in Tubes of Circular Cross Sec- 
tion, by Marcel Robin (Aerothermal Se- 
minar of the Faculté des Sciences de Paris 
Année 1954-1955), France, Ministére ce 
V Air, Pub. Sci. Tech., Notes Tech. 63, 1957, 
110 pp., pp. 13-27 (in French). 

Skin Heating During Re-entry of Satel- 
lite Vehicles to the Atmosphere, by T. h. 
F. Nonweiler, J. Brit. Interplan. Soc., vol. 
16, Jan.-March 1957, pp. 10-21. 

The Surface Temperature of a Missile 
in High Speed Flight, by P. L. Chambre, 
Calif. Inst. Tech., Jet Propulsion Lab., 
Progress Rep. 4-13, May 1945, 17 pp., 5 
figs. (Declassified from Secret by author- 
ity of Calif. Inst. Tech., Jet Propulsion 
Lab., Pub. 38, p. 7, 1/3/55.) 

The Skin Temperature of the Bumper 
WAC in High Speed Flight, by J. Lorrell, 
Calif. Inst. Tech., Jet Propulsion Lab., Rep. 
4-37, July 1947,22 pp. (Declassified from 
Confidential by authority of Calif. Insi. 
Tech., Jet Propulsion Lab., Pub. 38. p. 6, 
1/3/56.) 

Heat-Transfer Problems of Liquid- 
Cooled Gas-Turbine Blades, by Henry 
Cohen and F. J. Bayley, Proc., Inst. Mech 
Engrs., vol. 169, no. 53, 1955, pp. 1063 
1080. 

Oblique Flow Losses in Screen Matrix 
Heat Exchangers, by W. H. Comtois, 
Stanford Univ., Dept. Mech. Engng., Con- 
tract Nonr 225 (23), Tech Rep. 29, June 
1956, 72 pp. 

Research on Application = Cooling to 
Gas Turbines, by J. B. Esgar, J.N.B. 
Livengood, and R. O. Hickel, Trans., 
ASME, vol. 79, April 1957, pp. 645- 
652. 

Generalized Optimal Heat-exchanger 
Design, by D. H. Fix and R. R. Mills Jr., 
Trans. ASME, vol. 79, April 1957, pp. 
653-662. 

Review of Industrial Applications of 
Heat Transfer to Electronics, by J. Kaye, 
Proc. Inst. Radio Engrs., vol. 44, Aug. 
1956, pp. 977-991. 

Criteria for Validity of Lumpted-Param- 
eter Representation of Ducting Airflow 
Characteristics, by T. R. Stalzer and G. J. 
Giedler, Trans. ASME, vol. 79, pp. 833- 
840. 

Heat Transfer in the Critical Region. 
by R. P. Bringer and J. M. Smith, J. 
Amer. Inst. Chem. Engrs., vol. 3, March 
1957, pp. 49-55. 

Theoretical Investigation of the In- 
fluence of Cooling Wall Height on the 
Heat Transfer of Condensing Droplets, by 
R. Gregorig, Chemie-Ing.-Technik, vol. 28, 
Aug.-Sept. 1956, pp. 551-553 (in German). 

Heat Transfer to Boiling rd Liquid 
Mixtures, by W. R. van Ww ijk, A. S. Vos, 
S.J.D. van Stralen, Chem. Engng. Sci., vol. 
5, no. 2, April 1956, pp. 68-80. 


Heat Transfer in Pipe Flow at High 
Speeds, by J. E. Bialokoz and O. A. Saun- 
ders Combustion and Boilerhouse Engng., 
vol. 10, no. 11, Nov. 1956, pp. 412-420. 

Heat Transfer from Turbulent Flowing 
Gas in a Tube Inlet, by G. Grass, Allgem. 


JET PROPULSION 


~ A 
cy! 
A 
2 / 
} 
+ 
> 
= 


"A Rep. 
CA TN 


field of 
0 of 30, 
rickson 
pp. 

tal In- 
mal Se- 
e Paris, 
stére de 
ch. 63, 


Metals 
ss Sec- 
nal Se- 
e Paris 
tére cle 
1957, 


Satel- 
vol. 


Missile 
umbre, 

Lab., 
pp., 5 
uthor- 
ulsion 


Warmetechnik, vol. 7, no. 2, 1956, pp. 58- 
64 (in German). 

Pressure Loss in Pipes of Cross-Section 
at High Velocities, by A. Naumann 
Allgem. Warmetechnik, vol. 7, no. 2, 1956, 
pp. 32-41 (in German). 

Effect of Axial Fluid Conduction on Heat 
Transfer in the Entrance Regions of Paral- 
lel Plates and Tubes, by P. J. Schneider 
Trans. ASME, vol. 79, May 1957, pp. 765- 


i74. 


Jet and Rocket - 


Propulsion Engines 


The Problem of Goddard and Optimum 
Thrust Programming, by Frank D. Faulk- 
ner, Naval Postgrad. School. Res. Pap. 11 
March 1957, 10 pp. 

Underwater Propulsion, by Jack W 
Hoyt, G. G. Gould, J. F. Brady, S. Wolf 
and R. M. Dunlap, Missiles and Rockets 
vol. 2, May 1957, pp. 71-72, 75-77. 

Lycoming T53 Gas Turbine, by Anselm 
Franz, SAE J., vol. 65, May 1957, p. 25. 

Turboprops, by Brooke E. Allen, SAE J. 
vol. 65, May 1957, p. 31. 

Orion Gas-Generator Turbocompound 
Engine Provided Unusual Development 
Problems, by Ralph J. Hooker, SAE J., 
vol. 65, May 1957, pp. 51-56. 

Boeing 502-10c Gas Turbine, by Wal- 
lace E. Skidmore, SAF J., vol. 65, May 
1957, pp. 57-59. 

Processes in Pulse Jet Engines, by 
Fritz Staab, Z. Flugwissenschaften, vol. 5, 
March 1957, pp. 77-87, 22 ref. (in Ger- 
man). 

Investigation of Acceleration Character- 
istics of a Single-Spool Turbojet Engine, 
by Frank L. Oppenheimer and George J. 
Pack, NACA RM E53H26, Oct. 1953, 29 
pp. (Declassified from Confidential by 
authority of NACA Res. Abstracts 115, p. 
12, 5/22/57.) 

Altitude Performance of Compressor, 
Combustor, and Turbine Components of 
XT38-A-2 Turboprop Engine, by Frederick 
W. Schulze and William R. Prince, NACA 
RM E54E04, Aug. 1954, 37 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abstracts 115, p. 12, 5/22/ 
57.) 

Temperatures in a J47-25 Turbojet- 
Engine Combustor and Turbine Sections 
during Steady-State and Transient Opera- 
tion in a Sea-Level Test Stand, by C. R. 
Morse and J. R. Johnston, NACA-RM 
E54K30a, June 1955, 45 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abstracts 115, p. 12, 5/22/57). 

Analysis of a Form of Peak Holding 
Control, by G. J. Delio, NACA RM 
E56B10, March 1956, 57 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abstracts 115, p. 12, 5/22/57.) 

Removal of Secondary-Flow Accumula- 
tions in a Two-Dimensional Turbine 
Nozzle Passage by Boundary-Layer Bleed, 
by Robert Y. Wong, NACA RM E55E11, 

June 1955, 22 pp. (Declassified from 
Confidential by authority of NACA Res. 
Abstracts 115, p. 12, 5/22/57.) 

An Experimental Investigation of the 
Effect of High-Pressure Tail-Pipe Length 
on the Performance of Solid-Propellant 
Motors for Rocket-Powered Aircraft, by 
Charles J. Rodriguez, NACA RM L52E- 
12a, Aug. 1952, 37 pp. (Declassified from 
Confidential by authority of NACA Res. 
Abstracts 115, p. 16, 5/22/57.) 

Seal Leakage in the Rotary Regenera- 
tor and Its Effect on Rotary-Regenerator 
Design for Gas Turbines, by D. B. Har- 
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*The only valve that improves with use! 


VALCOR’s patented ‘‘Floating Seal’’ solenoid valve actually 
improves with use, due to its unique self-lapping action. 

The ‘‘heart’’ of this valve is a precise, optically flat, carbon- | 
graphite disc which floats in the plunger. A slight pressure, 
from either direction, moves the disc against 

an equally flat, stainless steel seat, 

sealing perfectly. 


The imaginative engineering which created the ‘‘Floating Seal’ has been 
responsible for the use of over a half million VALCOR solenoid valves in 
Military Aircraft and Guided Missiles. Hundreds of thousands more are in 
constant use throughout industry. Tell us your valve problems... . let our 
engineering staff offer its recommendations. Write today for Catalog 105 


and other valve information. 1a 


VALCOR ENGINEERING CORP. 
SOLENOID VALVES 


5371 CARNEGIE AVENUE, KENILWORTH, NEW JERSEY 
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TURBINES, PUMPS, COMPRESSORS 


Men who know 


oe Get in on the development of the most a 
advanced high-speed rotary equipment ever built 


Here’s your opportunity to step up to a new, higher level of 
turbo-machinery technology—the Large Rocket Engine. 

Whatever type of rotating machinery you know best, your 
experience could be extremely valuable in the important develop- 
ments now going on at Rocketdyne. The seasoned and ambitious 
man who has cut his teeth on jet engines, steam or gas turbines, 
or other elements of rotating machinery, is urgently needed to 
apply mechanical principles to meet the increasing demands of 
power plant performance. 

The combination of high speed, light weight, heavy loadings 
and exceptional pressures required in rocket engine work is lead- 
ing to an entirely new breed of high-performance rotating 
machinery...and a new breed of engineer. You can be one of 
this advance guard of the turbo-machinery field—if you have the 
desire to build your professional status by accepting new 
challenges. 

Rocketdyne is building high-thrust rocket engines for the 
nation’s major missiles. You'll work with the leading producer in 
the nation’s fastest growing industry. You and your work will 
be recognized as a vital part of the overall achievements. Testing 
facilities are among the world’s finest. The power produced is 
beyond anything ever before thought possible. If you would like 
to tackle new assignments working alongside some of the finest 
minds in turbo-engineering, write and tell us about your back- 
ground: A. W. Jamieson, Rocketdyne Engineering Personnel 
Dept. R-2, 6633 Canoga Avenue, Canoga Park, California. 


ROCKETDYNE 


A DIVISION OF NORTH AMERICAN AVIATION 


BUILDERS OF POWER FOR OUTER SPACE 


per, Trans. ASME, vol. 79, Feb. 1957, 
pp. 233-245. 

Studies for New Hot Air Engine—3, by 
H. A. Havemann, N. N. Narayan Rao, P. 
Jayachandra, and G. C. Garg, J. Indian 
Inst. Sci., Sec. B, vol. 38, July 1956, pp. 
172-202. 

Experimental Investigation of Turbine 
Stator -Blade- Outlet Boundary - Layer 
Characteristics and a Comparison with 
Theoretical Results, by Warren J. Whit- 
ney, Warner L. Stewart and James W. 
Miser, NACA RM _ E55K24, March 
1956, 24 pp. (Declassified from Confiden- 
tial by authority of NACA Res. Abs. 316, 
p. 29, 7/5/57.) 

Distribution of Losses Behind a Com- 
pressor Rotor as Measured by a Rotating 
Rake, by William R. Godwin, NACA 
RM L55F29, Jan. 1956, 66 pp. (Declassi- 
fied from Confidential by authority of 
NACA Res. Abs. 116, p. 23, 7/5/57.) 

Method for Estimating Combustion 
Efficiency at Altitude Flight Conditions 
for Combustor Tests at Low Pressures, 
by Walter T. Olson, J. Howard Childs 
and Wilfred E. Scull, NACA RM 
E53F17, Aug. 1953, 15 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abs. 116, p. 16, 7/5/57.) 


Application of Stream-filament Tech- 
niques to Design of Diffuser between Com- 
pressor in a Gas-Turbine Engine, by Nor- 
bert O. Stockman, NACA RM E55F06, 
Aug. 1955, 15 pp. (Declassified from Con- 
fidential by authority of NACA Res. Abs. 
116, p. 20, 7/5/57.) 

An Investigation of the High Alternat- 
ing Stresses in the Blades of an Axial 
Flow Compressor, by J. R. Forshaw, 
Gt. Brit. Aeron. Res. Council, Rep. & 
Mem. 2988 (formerly ARC Tech. Rep. 
17132; Nat. Gas Turbine Estab. Rep. 
R157), 1957, 19 pp. 

Satellite Rocket Power Plant, by J. 0. 
Crum and §. L. Gendler, Rand Corp. 
RA15027, Feb. 1957, 110 pp. (Declassi- 
fied from Secret by authority of ASTIA 
TAB U157, p. 57, 6/24/57.) 

High-Temperature Lubricants and 
Bearings for Aircraft Turbine Engines, 
NACA Subcommittee on Lubrication and 
Wear. Appendix A: High-speed Air- 
craft Missions, by C. M. Michaels, Wright 
Air Dev. Center; Appendix B: Engine 
Design Trends Affecting Lubricants and 
Bearings, by C. C. Singleterry, Bur. 
Aeron., Dept. of the Navy; Appendix C: 
Problems Encountered at High Tempera- 
tures in Lubrication Systems of Turbine 
Engines, by G. P. Townsend, Westing- 
house Electric Corp; Appendix D: 
Turboprop Gear Lubrication Problems, 
by C. J. McDowall, Gen. Motors Corp; 
Appendix E: Notes on High-Tempera- 
ture Fluids and Lubricants, by E. EF. 
Klaus and M. R. Penske, Penn. State 
College; Appendix F: NACA Research 
on Lubricants, Bearings, and Lubrication 
for High-Temperature Turbine Engines, 
by R. L. Johnson and E. E. Bisson; 
Appendix G: High-Temperature Bearing 
Problems, by F. W. Wellons. NACA 
RM E54D27, July 1954, 101 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abs. 116, p. 18, 7/5/57.) 

An Experimental Investigation of a 
Flat Ram-jet Engine on a Helicopter 
Rotor, by Robert D. Powell, Jr., and 
James P. Shivers, NACA RM L55F28, 
Jan. 1956, 27 pp. _(Declassified from 
Confidential by authority of NACA 
Res. Abs. 116, p. 23, 7/5/57.) 
Experimental Determination of Aero- 
dynamic Forces Normal to the Chord Due 
to Rotating Stall Acting on Compressor 
Blading, by Donald F. Johnson and 
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Costilow, NACA RM E57F 14, 
(Declassified from Con- 
NACA Res. 


Eleanor L. 
Aug. 1954, 27 pp. 
fidential by authority of 
Abs. 116, p. 18, 7/5/57.) 

Noise Survey of a Full-Scale Super- 
sonic Turbine-Driven Propeller Under 
Static Conditions, by Max C. Kurbjun, 
VACA TN 4059, July 1957, 20 pp. 

Performance of Inconel 550 Turbine 
Blades in a Turbojet Engine and Effects 
of Different Forging Temperatures and 
Heat Treatments, by C. A. Gyorgak, J. 
R. Johnston, and J. W. Weeton, NACA 
RM E55F08, Aug. 1955, 55 pp. (De- 
classified from Confidential by authority 
by NACA Res. Abs. 116, p. 20, 7/5/57.) 

Equilibrium Operating Conditions of 
Turbine-Compressor Gas Generators, by 
\. Capetti, L’ Aerotechnica, vol. 37, Feb. 
1957, pp. 3-6. 

Static Thrust of Turbojet with and with- 
out Afterburner; Limitations of Increas- 
ing the Compression Ratio of the Com- 
pressor, by G. Santangelo, L’ Aerotecnica, 
vol. 37, Feb. 1957, pp. 24-28. 

Inflow Conditions for Supersonic Com- 
pressor with Curved Blades, 
Levine, J. Applied Mech., vol. June 
1957, pp. 165-169. 

Design Progress: Napier Gazell—Power 
Unit for Helicopters, by Randolph Haw- 
thorne, Aviation Age, vol. 28, July 1957, 
pp. 32 33. 

Reds’ Type K Outperforms Western 
Turboprops, by Walter P. Moser, Avia- 
tion Age, vol. 28, July 1957, pp. 50-55. 

Space Dynamics Techniques Hinge on 
Propulsion System, by John de Nike, 
Jorgen Jensen and Michael Stoiko, Avia- 
lion Age, vol. 28, July 1957, pp. 26-31. 

Research and Development Trends, 
Aviation Age, Research and Development 
Technical Handbook, 1957-1958, Section 
C. Propulsion, pp. C-3-—C-9;(Chemical 
Rockets-Photonic & Ionic Rockets-Gas 
Turbines & Ramjets-Fuels & Lubricants). 

Gas Turbines for Aircraft (Table), 
Aviation Age, Research and Development 
Technical Handbook, 1957-1958, pp. C-10- 
C-17; Data Review, pp. C-18, 22-41. 

O-Ring Valve Used in Rockets, Avia- 
lion Week, vol. 67, July 15, 1957, p. 63. 

Injection by Concurrent Jets in Liquid 
Rocket Combustion Chambers, by R. 
Kling, Mme. G. Chevalerias and A. 
Maman, Recherche Aeron., no. 57, 1957, 
pp. 19-22 (in French). 

Turbo-Prop Performance Analysis and 
Flight Planning—A Viscount Operator’s 
Approach, by J. R. Baxter, J. Royal 
Aeron. Soc., vol. 61, June 1957, pp. 391- 
404, 

Ram-jets, by R. R. Jamison, J. Royal 
Aeron. Soc., vol. 61, June 1957, pp. 407 
421. 

Ramjet-Turbojet Combinations for 
Supersonic Flight, by R. T. DeVault, 
SAE J., vol. 65, July 1957, pp. 35-38. 

Rockets-on-Rotor System, by Donald 
S. Chatfield, SAE J., vol. 65, July 1957, 
pp. 51-52. 

Small Turbojets Cover Mach 2-3 
Range at 50,000-70,000 Ft. Altitudes, 
by Reece V. Hensley, Stanley R. Shapiro 
and Raymond Capiaux, SAE J., vol. 65, 
July 1957, pp. 90-93. 

Ramjets Come to Fore as Component 
Design Improves, by F. R. Garbarine, 
Aviation Age, vol. 28, Aug. 1957, pp. 66— 

Turboprop v. Turbojet, by Frank Robert- 
son, Flight, vol. 72, Aug. 2, 1957, p. 157. 


Gyron, Flight, vol. 72, Aug. 2, 1957, 
yp. 163-168. 
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The Martin USAF B-57 tactical bomb- 
er carries wing-tip fuel tanks to in- 
crease its range. When the fuel in 
them is exhausted, the tanks are jet- 
tisoned by actuating explosive bolts 
with Du Pont electric detonators (left). 


How Du Pont Detonators help the 
Martin B-57 get into battle trim 


When the Air Force pilot of Martin’s 
B-57 trims his ship for peak fighting 
efficiency, he drops his wing-tip fuel 
tanks. 

By simply pressing the control but- 
ton, these tanks are easily and quickly 
jettisoned by an explosive bolt using 
a Du Pont Electric Detonator. 

One model of the B-57 is equipped 
with a flip-top “life saver.” In the event 
the crew has to bail out, Du Pont Deto- 
nators blow the canopy off a split-sec- 
ond before the crew is catapulted into 
space. 

On still another model, a mechanical 
“hand” actuated by a Du Pont Detona- 
tor “snatches” the control column for- 
ward under the instrument panel. This 
protects the pilot’s legs from injury as 


he is shot out of the plane to safety. 

These three applications illustrate 
the unusual problems being solved by 
Du Pont’s complete line of explosive 
specialties . . . including high-tempera- 
ture aircraft expansion rivets; MS de- 
lay connectors for timing and firing 
explosive charges; explosive releases 
and a wide variety of electric squibs 
for igniting powder charges under con- 
trolled conditions. 

Maybe you have an unusual prob- 
lem that can be solved by one of these 
Du Pont Explosive Specialties. We'll be 
happy to send you more information. 
Just write to E. I. du Pont de Nemours 
& Co. (Inc.) , 2543-Nemours Bldg., Wil- 
mington 98, Delaware. 


DU PONT EXPLOSIVE SPECIALTIES — 


PRODUCTS OF DU PONT RESEARCH 


Better Things for Better Living . . . 


REG. U.S. PAT. OFF. 


through Chemistry 


3 SEE THE “DU PONT SHOW OF THE MONTH” ON CBS 
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“The might 
of 


voodoo 


was proved 
with Statham 


Transducers 


IN THE F-101' 


Statham transducers gathered 
pressure and acceleration 
data in McDonnell Aircraft 
Corporation's program to 
prove out the aerodynamic 

design and structural design. 


WHEN THE NEED 
IS TO KNOW...FOR SURE 
SPECIFY STATHAM 


Accelerometers « Load Cells 
Pressure Transducers 


Catalog, complete with prices, 
available upon request. 


INSTRUMENTS, INC. 
LOS ANGELES 64 


Index to Advertisers 


AEROJET-GENERAL CORPORATION....... 
D'Arcy Advertising Co., Los Angeles, Calif. 


AMERICAN MAcuiInE & Founpry Co 
Fletcher D. Richards, Inc., New York, N. Y. 


ApPLIED Puysics LABoRATORY, THE JoHNS Hopkins UNIVERSITY 
M. Belmont ver Standig, Inc., Washington, D. C. 


Avco MANUFACTURING CORPORATION, RESEARCH AND ADVANCED DEVELOp- 


MENT DIVISION................. 
Benton & Boles, Inc., New York, N. Y. 


Becco CuHEemicaL Division, Foop Macutnery & CHEMICAL CORPORATION. 
John Mather Lupton Co., New York, N. Y. 


Bristot Company 
James Thomas Chirurg Co., New York, N. Y. 


Century EvEctronics & IstrRuMENTs, INC 
Gibbons Advertising Agency, Inc., Tulsa, Okla. 


DIvEeRsEY ENGINEERING COMPANY 
Roark & Colby Advertising, Chicago, Il. 


pu Pont pE Nemours, E. I., anp ComPpANY 
ExpLosIvEs DEPARTMENT 
Charles L. Rumrill & Co., Inc., Rochester, N. Y. 


Puoto Propucts 
V. W. Ayer & Son, Inc., Philadelphia, Pa. 


EASTMAN Kopak ComPANY.. 
The Rumrill Company, Inc., Roe hester, 


DEVELOPMENTS, INC... 


THE GARRETT CORPORATION 
AIRESEARCH MANUFACTURING DIVISIONS 
J. Walter Thompson Co., Los Angeles, Calif. 


MACHINE WORKS... 
R. W. Sayer Co., Chicago, Ill. 


LocKHEED AIRCRAFT,COMPANY, MISSILE SysteMs Division 
Hal Stebbins, Inc., Los Angeles, Calif. 


Maraquarpt AIRCRAFT CoMPANY... 
Grant Advertising, Inc., Hollywood, Calif. 


THE Norton Company .. Third Cover 
James Thomas Chirurg Co., Boston, Mass. 


Al Paul Lefton Co., Philadelphia, Pa. 


RHEEM MANUFACTURING COMPANY 
Getz & Sandborg, Inc., Beverly Hills, Calif. 


RockETpYNE, A Division oF NorTH AMERICAN AVIATION, INC.. 
Batten, Barton, Durstine & Osborn, Inc., Los Angeles, Calif. 


SERVOMECHANISMS, INC 
Sanger-Funnell, Inc., New York, N. 


SpacE TECHNOLOGY LABORATORIES, THE RAMO-WOOLDRIDGE CORPORATION 
The McCarty Company, Los Angeles, Calif. 


StaTHAM INSTRUMENTS, INC. 
Compton Advertising, Inc., Los ‘Angeles, Calif. 


THIOKOL CHEMICAL CoORPORATION......... 
Dancer-Fitzgerald-Sample, Inc., New York, N. Y. 


VaALCOR ENGINEERING CORPORATION... . 
Keyes, Martin & Company, Newark, N. J. 


WyMANn-GorDON CoMPANY 
John W. Odlin Co., Inc., Worcester, Mass. 


JET PROPULSIO! 


>. 4 
* 
‘ 

; 
: 
: 
j 

29 
3 

~ 77 
ek 

5 


Cover 


137 


76 


136 


131 


148 


Jover 


145 


80 


